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Chapter 1 
Introduction 


In this thesis the polarized and nnpolarized photon distributions of the nucleon (proton, 
neutron), evaluated in the equivalent photon approximation, are computed theoretically 
and the possibility of their experimental determination is demonstrated. The thesis is 
based on the following publications [1-6]: 

• M. Gliick, C. Pisano, E. Reya, The polarized and unpolarized photon content of the 
nucleon, Phys. Lett. B 540, 75 (2002) [Chapter 3]. 

• M. Gliick, C. Pisano, E. Reya, I. Schienbein, Delineating the polarized and unpola¬ 
rized photon distributions of the nucleon in eN collisions, Eur. Phys. J. C 27, 427 
(2003) [Chapter 4]. 

• A. Mukherjee, C. Pisano, Manifestly covariant analysis of the QED Compton process 
in ep — > C'-yp and ep —> eyX, Eur. Phys. J. C 30, 477 (2003) [Chapter 5]. 

• A. Mukherjee, C. Pisano, Suppressing the background process to QED Compton 
scattering for delineating the photon content of the proton, Eur. Phys. J. C 35, 509 
(2004) [Chapter 6]. 

• A. Mukherjee, C. Pisano, Accessing the longitudinally polarized photon content of 
the proton, Phys. Rev. D 70 , 034029 (2004) [Chapter 7]. 

• C. Pisano, Testing the equivalent photon approximation of the proton in the process 
ep uWX, Eur. Phys. J. C 38, 79 (2004) [Chapter 8]. 



Introduction 




The results of [5] have been summarized in [7] as well. Furthermore in Section 2.5 we 
shortly recall the main hndings of the recent paper [8], also concerning the structure of 
the nucleon: 

• M. Gliick, C. Pisano, E. Reya, Probing the perturbative NLO parton evolution in 
the small-x region, Eur. Phys. J. C 40, 515 (2005). 

The equivalent photon approximation (EPA) of a charged fermion is a technical device 
which allows for a rather simple and efficient calculation of any photon-induced subpro¬ 
cess. The first explicit formulation and quantitative application of the EPA were given in 
1924 by Fermi [9], who utilized it to estimate the electro-excitation and electro-ionization 
of atoms, and also the energy loss, due to ionization, of a-particles travelling through mat¬ 
ter. In several cases, he obtained a satisfactory numerical agreement with experimental 
data. Ten years later, in order to simplify calculations of processes involving relativistic 
collisions of charged particles, Williams [10] and Weizsacker [11] further developed Fermi’s 
semi-classical treatment and extended it to high-energy electrodynamics. They observed 
that the electromagnetic field generated at a given point by a fast charged particle passing 
close to it contains predominantly transverse components. By making a Fourier analysis 
of the field, they concluded that the incident particle would produce the same effects as 
a beam of photons and computed their distribution in energy. This model assumes that 
the particle motion is not appreciably affected during the interaction, in particular that 
its scattering angle is small. 

The first field-theoretical derivations were given in the fifties by Dalitz and Yennie [12], 
Curtis [13], Kessler and Kessler [14], and later by Chen and Zerwas [15]. The equivalent 
photon method for pointlike fermions, like the electron, has been investigated and utilized 
widely, for example it has been applied to pion production in electro-nucleon collisions 
[12,13] and to two-photon processes for particle production at high energies [16-18]. A 
detailed history of the method, its various formulations and its first applications are 
contained in Kessler’s review article [19]. In more recent times, the validity of this ap¬ 
proximation has been examined by Bawa and Stirling [20] in the context of the production 
of large transverse momentum photons at the HERA collider, by comparing the exact and 
approximate cross sections. Frixione, Mangano, Nason and Ridolfi in a well-known paper 
[21] further modified the method in order to improve its accuracy and tested it in the case 
of heavy-quark electroproduction at HERA. Their results were extended by De Florian 




3 



Figure 1.1: Connection between electroproduction and photoproduction on a target a. 


and Frixione [22] to the case of a longitudinally polarized electron. 

As an example to illustrate the EPA, we shall consider the process in which an electron 
e of very high energy scatters from a target hadron a, e.g. a proton. At leading order in a, 
the hne structure constant of quantum electrodynamics (QED), the electron is connected 
to the target by one photon propagator, as depicted in the hrst diagram in Figure 1.1. 
If we denote with E and E' respectively the initial and final energies of the electron, the 
photon will carry a momentum q such that 

~ — 2 E'E '(1 — cos 6 *), ( 1 . 1 ) 

where 6 is the electron scattering angle. In the limit of forward scattering, whatever the 
energy loss, the photon momentum approaches = 0 ; therefore the reaction is highly 
peaked in the forward direction and the underlying dynamics is that of a photoproduction 
process. It can be shown that the cross section Uea, integrated over g^ with integration 
bounds 

2 2 2 ON 

W = , ?min = -/^ , (1-2) 

is given by a convolution of the probability that the electron radiates off a photon, the 
equivalent photon distribution 7 ^( 0 ;,/i^), with the corresponding real photoproduction 
cross section which is in general easier to calculate than the original one: 

dUea = p^) a^adx. (1.3) 

The connection between electroproduction and photoproduction on a target hadron a is 
shown, in terms of Feynman diagrams, in Figure 1.1. An explicit expression of 7 ^( 0 ;,/r^) 
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is given in [21], 




a 


1 + (1 — /i^(l — x) 


In 


X 


mlx'^ 


+ 2mix 


X 




mlx'^ 


(1.4) 


where x is the fraction of the electron energy carried by the photon and /i has to be 
identihed with a momentum scale of the photon-induced subprocess. While in (1.2) 
is the largest value of kinematically allowed, is not a well dehned quantity, since 
when the momentum transfer squared becomes too large (in absolute value), the EPA 
breaks down and use of (1.3) would lead to huge errors. Critical examination of the EPA 
in electron-hadron collisions, in connection with different choices for the scale fi, can be 
found in [20,21]. 

The longitudinally polarized electron-target a cross section Aaea can be obtained from 
(1.3) with the formal substitutions 


a 


Act, Ay®. 


(1.5) 


The quantities Aa are dehned in terms of cross sections crxj^x^ for incoming particles of 
dehnite helicities, 

dAa = - (d(T++ -|- da_— da+_ — da_+) 

= ^ (dc^++- da+_), (1.6) 

where the second equality follows from parity invariance of the electromagnetic interaction. 
The corresponding unpolarized cross section is obtained by taking the sum instead, namely 

= ^(dcr++ + da+_). (1.7) 

The polarized equivalent photon distribution Ay® are dehned in terms of densities for 
photons of dehnite helicity in electrons of dehnite helicity 

A7' = 7? - 7V = 7r - 7r. (1-8) 

while the sum will give the unpolarized photon distribution, 

y- = y;++yl+ = yr +y®-, (1.9) 


where the superscripts refer to the parent electron and the subscripts to the photon. One 
can show that [22] 


Ay®(a;,/i2) 


a 1 — (1 — x)^ 

27r X 


In 


/i^(l — x) 
mlx"^ 


-I- 2m\x^ 



1 

V/i2 

m‘lx‘^ ) \ 


( 1 . 10 ) 
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In the treatment of protons and neutrons, (1.3) and (1.5)-(1.9) still hold, with the 
replacement e —p or n, but a special situation arises in the calculation of their photon 
distributions, due to the fact that they are not pointlike particles. Here it is necessary 
to distinguish between elastic and inelastic scattering. In the former case the nucleon 
does not break up but is temporaly in an excited state, which can be described in terms 
of certain form factors. In the latter case, appealing to the parton model, only quarks, 
antiquarks and gluons are usually considered as the (essentially free) constituents of the 
initial nucleon, which ceases to exist and its constituents dually hadronize. The photons 
radiated off the quarks and antiquarks may be characterized by (1.4) and (1.10), and 
due to their logarithmic enhancement factor they are expected to become increasingly 
relevant at very high energies, when the momentum scale p is also large. At this level the 
“inelastic” photons can be included among the parton distributions of the nucleon. 

Therefore the total (polarized) unpolarized photon distribution of a nucleon N will be 
given by 

(A) 7 (a;, p^) = (A) 7 ei(a;) + (A) 7 inei(a;, p^), (1.11) 

where the elastic component (A) 7 ei is due to A —> 'jN and the inelastic component 
(A) 7 inei is due to A — 7 A with A 7 ^ A. As pointed out by Drees and Zeppenfeld [23] 
and by Kniehl [24], 7 ^^ cannot be obtained from (1.4) by just replacing the electron mass 
with the proton mass: this would strongly overestimate it. One has to take into account 
the effects of the form factors, which determine the scale independence of (A) 7 ei. The 
derivations of and performed in [24] and by Gliick, Stratmann, Vogelsang in [25], 
will be generalized to the neutron and to the polarized sector in Chapter 3. 

The photon content of the nucleon 7 ( 0 ;, p^) can be utilized, instead of the more common 
form factors and parton distributions, to calculate photon-induced subprocesses in elastic 
and deep inelastic ep reactions, leading to great simplihcations. For example, as shown 
by Bliimlein [26] and by De Riijula and Vogelsang [27], the analysis of the deep inelastic 
QED Compton scattering process ep eyA reduces to the calculation of the 2 —2 
subprocess ey^ —ey instead of having to calculate the full 2 —3 subprocess eg —eyg. 
Analogously one can consider just the simple 2 —2 subprocess y®y^ —(instead of 
y®g —> fj,~^p,~q) for the analysis of deep inelastic ep or ey^ —uW (instead of 

eg — uWq) for associated uW production in ep uWX. 

Similarly, yy fusion processes like y^’yP —> cc, for (heavy) lepton 

(£), heavy quark (c), charged Higgs (A^) and slepton (£) production can be easily ana- 
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lyzed in purely hadronic pp reactions, providing also an interesting possibility of producing 
charged particles which do not have strong interactions. Carlson and Lassila [28], Drees, 
Godbole, Nowakowski, Rindani [29], and Ohnemus, Walsh, Zerwas [30] initiated these 
studies; in particular, in [29] it is shown that the cross section for the pair production of 
heavy charged scalars or fermions via 77 fusion amounts to a few percent of the corre¬ 
sponding Drell-Yan qq annihilation cross sections, at energies reached at the CERN Large 
Hadron Collider (LHC). However, the disadvantage of the low production rates is com¬ 
pensated by the simple and clean experimental situation encountered when the photons 
are emitted from protons which do not break up (purely elastic processes) [30]. 

It still remains to extend the above-mentioned analysis to pd and dd reactions. More¬ 
over, analogous remarks hold for the longitudinally polarized eN and pp, pd, dd reactions 
where the polarized photon content of the nucleon A 7 (a;, p^) enters. 

In this thesis we shall concentrate on lepton-nucleon scattering, with special attention 
to the aforementioned QED Compton process, which is one of the most important reac¬ 
tions for directly measuring the photon content of the nucleon and testing the reliability 
of the EPA. Furthermore, we shall discuss how such measurements can provide additional, 
independent informations concerning the usual polarized and unpolarized structure func¬ 
tions, along the lines of [26,31-33]. 

In detail, the outline of the thesis will be as follows: 

• Chapter 2 serves as an introduction into the concepts of electromagnetic form fac¬ 
tors, structure functions and parton distributions of the nucleon, in terms of which 
elastic and inelastic eN cross sections are commonly expressed. They are presented 
here because they enter in the dehnition of the photon content of the nucleon. No¬ 
tations and conventions used in this chapter as well as in the rest of the thesis are 
summarized in Appendix A. 

• In Chapter 3 a new expression for the polarized equivalent photon distribution 
is explicitly derived. The calculation of the unpolarized one, already performed 
in [24,25,27], is also shown in detail for completeness and the resulting photon 
asymmetries are presented for some typical relevant momentum scales. 

• The production rates of lepton-photon and dimuon pairs at the HERA collider and 
HERMES experiment are evaluated in Chapter 4, utilizing the photon distributions 
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previously derived, convoluted with the cross sections relative to the subprocesses 
ey —> ey and yy —given in Appendix B. It is shown that the production 
rates are sufficient to measure the polarized and unpolarized photon content of the 
nucleon. 

• Chapter 5 is devoted to the unpolarized QED Compton scattering in ep —> eyp and 
ep eyX, with the photon emitted from the lepton. The full 2 —3 process is 
calculated in a manifestly covariant way by employing appropriate parametrizations 
of the proton’s structure functions. These results are compared with the ones based 
on the EPA, as well as with the experimental data and theoretical estimates for the 
HERA collider given in [31]. It is shown that the cross section is reasonably well 
described by the EPA of the proton, also in the inelastic channel. In addition it 
turns out that the results obtained in [31], based on an iterative approximation pro¬ 
cedure proposed by Conran and Kessler [34], deviate appreciably from our analysis 
in certain kinematical regions. Details about the kinematics of the process are given 
in Appendix C. 

• In Chapter 6 the virtual Compton scattering process in ep eyp and ep eyX, 
where the photon is emitted from the hadronic vertex, is investigated. It represents 
the major background process to QED Compton scattering. New kinematical cuts 
are suggested in order to suppress the virtual Compton scattering background and 
facilitate the extraction of the equivalent photon distribution of the proton at the 
HERA collider. The analytic expressions of the matrix elements of unpolarized 
QED and virtual Compton scattering can be found in Appendix D. 

• Chapter 7 is devoted to the QED Compton process in longitudinally polarized 
lepton-proton scattering. The kinematical cuts necessary to measure the polari¬ 
zed photon content of the proton and to suppress the major background process 
coming from Virtual Compton scattering are provided for HERMES, COMPASS 
and future eRHIC experiments. We point out that these measurements will also 
give access to the spin dependent structure function gi in a kinematical region not 
well covered by inclusive measurements. The analytic expressions of the matrix ele¬ 
ments of polarized QED and virtual Compton scattering are relegated to Appendix 
E. 

• The accuracy of the EPA of the proton in describing the inelastic process ep uWX 
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is investigated in Chapter 8. In particular, the scale dependence of the correspon¬ 
ding inelastic photon distribution is discussed. Furthermore, an estimate of the 
total number of events, including the ones coming from the elastic and quasi-elastic 
channels of the reaction, is given for the HERA collider. 

• Finally, summary and conclusions can be found in Chapter 9. 




Chapter 2 

The Structure of the Nucleon 


In physics, one of the most common ways of getting information on the structnre of 
extended objects like hadrons is to use structureless particles as projectiles which scatter 
off the hadron in question. To probe the inside of a nucleon one naturally uses charged 
lepton (electron or muon) beams: such reactions dominantly take place by the exchange 
of a single virtual photon, therefore they can be studied with relative ease and clarity. 
Effects of higher order virtual photon exchange are believed to be less than a few percent. 
According to Heisenberg’s uncertainty principle, one can resolve the target structure down 
to the scale A ~ h/Q, where Q is the momentum transfer from the lepton to the nucleon. 
Hence, the higher the energy loss of the lepton, the hner the structure that can be resolved. 
The study of such reactions represent the main subject of the present chapter. 

In Section 2.1 we examine electron-nucleon elastic scattering and discuss the signih- 
cance of the electromagnetic form factors. In Section 2.2 inelastic electron-nucleon scat¬ 
tering is described in terms of the structure functions Fi and F 2 . The spin dependent 
structure functions gi and g 2 are introduced in Section 2.3. The parton model and quan¬ 
tum chromodynamics description of inelastic electron-nucleon collisions can be found in 
Section 2.4 and in Section 2.5 respectively. 


2.1 Elastic Form Factors 

We consider the elastic electron-nucleon scattering process: 


e{l) + N{P) ^ e{l') + N{P'), 


( 2 , 1 ) 
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Figure 2.1: 


The one-photon exchange diagram for elastic electron-nucleon scattering. 


where the four-momenta of the particles are given in the brackets. At lowest order in 
perturbation theory of QED the reaction is described by a one-photon exchange diagram, 
depicted in Figure 2.1. In terms of the Dirac 7 matrices and spinors, the electron transition 
current is 

(e(f) I J“„(0) I e(l)) = -e \n(l')r n(l)], (2.2) 

with e denoting the electric charge of the proton. The vertex linking the photon with the 
nucleon is not point-like, therefore the nucleon transition current is different from ( 2 . 2 ); 
it is the most general Lorentz four-vector that can be constructed from P, P' and the 7 
matrices, sandwiched between u{P) and u(P): 

( N{P') I J"(0) I N{P)) = e [u{P') F“ u{P) ]. (2.3) 

In F“ terms involving 7 ^ are ruled out by the conservation of parity. Furhermore, using 
the Dirac equations fui^P) = mu{P) and u{P')f = mu{P'), where m is the nucleon 
mass, one can show that there are only three independent terms, 7 “, ia^^qg and g", with 

\ |7". 7'’! (2.4) 

and q being the momentum transfer. 


q = l-l'. 


( 2 . 5 ) 




2.1 Elastic Form Factors 


11 


Therefore, quite generally, 

r“ = + F,{q^)q\ (2.6) 

The coefficients Fk{q^) {k = 1,2,3) are the electromagnetic elastic form factors of the 
nucleon, which are functions of the momentum transfer squared the only independent 
scalar variable at the nucleon vertex (being P ■ q = —q‘^/2). For the electromagnetic case 
we are interested in, is a conserved current and therefore the matrix element (2.3) 
has to satisfy the condition 


( N(P') I I N{P) ) = 0, (2,7) 

which implies F^{q^) = 0 in (2.6). The Dirac equation can be used again in order to 
replace u{P')i qf^u{P) by -u(P')[2m7" — (P + P')"]m(P) in (2.6) so that one can write 

) 7” - ^ F^(q^) (P + PT . (2.8) 

At = 0, which physically corresponds to the nucleon interacting with a static electro¬ 
magnetic field, the form factors are related to the electric charge Q and the magnetic 
dipole moment /r of the nucleon; 

ePi(0)=g, ^[Pi(0) + P 2 ( 0 )]=/x. (2.9) 

For an electrically neutral particle, like the neutron, one has Pi(0) =0. If a particle has 
no anomalous magnetic moment k, dehned such that [i = (1 -l- k) e/(2m), then P 2 ( 0 ) = 0. 
Experimentally, 

Pf (0) = 1, P|(0) = Kp ~ 1.79 (2.10) 


for the proton, and 


Ff{0) = 0, F^{0) = -1.91 


( 2 . 11 ) 


for the neutron. In the following, we will always express ja in units of the nucleon magneton 
eh/{2mc), that is /ip = 2.79, fin = —1.91. 

The amplitude M relative to the process (2.1) has the form 
M = -e" [u{l'hau(l)] i |u(F')r“„ti(/’)], 


( 2 , 12 ) 
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Taking the modulus squared of the amplitude and multiplying by the appropriate phase 
space and flux factors, one flnds that the differential cross section can be written as 


da 


1 

4a/ (/ ■ Py — m?ml 


|M|2 + P-1' -P') 


dH' d?P' 
(27r)3 2/^ (27r)3 2P^ ’ 


(2.13) 


which holds with the normalization of the spinors given in (A. 16) and where me denotes 
the electron mass. After integrating over the phase space of the scattered nucleon, one 
can use the condition P' = P + qto rewrite the energy conserving ^-function as 

i(/„ + Po - /i - p;,) = ti(Pa - (Po + 9o)t = i(P'" - (P + qf) 

= 5{q^ + 2P-q). (2.14) 


Therefore (2.13) reduces to 


da 


1 

4a/(/ ■ P)2 — m^TTlg 


|M|2 


d?l' 


6{q^ + 2P-q), 


(2.15) 


with 

W?=^tLcf(l\l')Hf(P-,P% (2,16) 

where the tensors La/ 3 {l;l') and H^f{P]P') come from averaging over initial spins and 
summing over flnal spins in the products of the electron and nucleon matrix elements 
(2.2) and (2.3), when (2.12) is squared. The completeness relation (A.17) is commonly 
used to perform the spin sums. More explicitly: 

^ E <'(''> I I '('))*<"('') I -^1(0) I e(i)) 

spins 

= ^e^Tr[(/'+me)7"(/ + ^^e)7^] 

= 2e^{ri'^+ -ml)] (2.17) 


and 


//.f(P;P') 


i {N{P') I J“.(0) I N{P))'{N(P') I 7„(0) I N{P)) 

spins 

^ Tr[{f+ m)r“(f + m)r/ 




{p + p')'^{p + pY 
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+ (Fi + F2)2(gV'-gV) 


(2.18) 


In the kinematical region of high energies, the electron mass can be neglected and then 
the contraction of (2.17) and (2.18) gives 


= 2e^|4(^F2 


p2 

4m2 2 


[2(P-/')(P-/) - 

2{F^ + F2fq\Pl')\. (2.19) 


In a frame in which the nucleon is at rest and the electron moves along the axis with 
energy E and is scattered into a solid angle fl = (0, ip) with hnal energy E', i.e. 


P = (m, 0,0,0,), / = P(l, 0,0,1), /' 


E'{1, sin 6 cos99, sin 0sin99, cosO), 


( 2 . 20 ) 


then 


g2 = {l-l'f = -dPP'sin^^ 


and one can dehne the energy transfer from the electron to the target 

q ■ P 


p = E-E' = 


m 


The differential cross section (2.15) can be rewritten as 


da 


a 


dP'dfi 4P'2 sin4 I 


Fl 


2 g , 
4m2 ^ 


cos^-^(Pi + P2)^sm^ - 

2 2mF ^ 2 


5 Z/ + 


( 2 . 21 ) 


( 2 . 22 ) 


2m 


, (2.23) 


with a = e^/(47r). The form factors Pi and P2, usually referred to as the Dirac and Pauli 
form factors respectively, parametrize our ignorance of the complicated structure of the 
nucleon. In practice, however, it is better to use linear combinations of them, the Sachs 
electric and magnetic form factors 


= Fi(«2) + F2tf). (2.24) 


dehned so that no interference terms, GeGmi occur in (2.23). At q^ = 0, one has 


(0) = 1, G^(0) = 1 + ~ 2.79, 

G^(0) = 0, Gl,{Q) = -1.91. (2.25) 




14 


The Structure of the Nucleon 


Integration of (2.23) over E', taking into account that depends on E' when 6 is held 
hxed, see (2.21), gives the Rosenbluth formula [35] 

.,2 


da 


a 


4E2sin^| E 


Gl + rGl, 2^ , o 0^2 .2^ 
cos^ - + 2tGI, sin^ - 


(2.26) 


with r = —q^/Aim? and the factor 


E' 

¥ 


2E , 

1 H-sin^ 

m 



(2.27) 


arises from the recoil of the target. The Rosenbluth formula is the basis of all experimental 
studies of the electromagnetic structure of the nucleon and allows to determine the form 
factors by measuring da/dfl as a function of 6 and Experimentally it turns out that 
the form factors drop rapidly as —g^ increases: 


GVjf) ^ _ GlM) 

flp fifl 

0<G^(g') <0.1 . 



0.71 GeV^ 


-2 


(2.28) 

(2.29) 


The function given in (2.28) is merely empirical and is often called a dipole £t; a monopole 
fit would be (1 — cg^)“^. Also there is no fundamental theoretical reason for G% 

to have the same g^ behaviour: this likeness is expressed by saying that the three form 
factors scale together in g^. 


In the non-relativistic limit, E ^ m, one can see from (2.27) that E = E'] therefore 
—g^ = |qp m4. In this limit the form factors Ge and Gm are the Fourier transforms 
of the nucleon’s charge and magnetic moment density distributions, respectively [36]. 
Assuming, for example, that the charge density distribution of the proton is spherically 
symmetric, i.e. a function of r = |r| alone, and that it is normalized such that 



(2.30) 


then the exponential in its Fourier transform can be expanded for small \q\ as follows 


= I d^rp(r)e*^'^ = / rp{r)i 1 + iq ■ r 


= 1 - |g| 2 (r 2 ) + ... , 

6 


[q-r) 


+ ... 


(2.31) 


where the mean square charge radius of the proton is dehned by 


d^r p(r)r‘^. 


( 2 . 32 ) 
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Figure 2.2: The one-photon exchange diagram for inelastic electron-nucleon scattering. 


Identifying (2.31) with the expansion of (2.28), 


GW) 


Ge^^) + 5 '^ 


( dGW) 

\ 



1 - 



( dGW) 

\ dg^ 


+ ... , (2.33) 

g2=0 


one gets 




(2.34) 


The same radius of about 0.8 fm is obtained for the magnetic moment distribution. Fi¬ 
nally, the charge distribution of the nucleon has an exponential shape in configuration 
space: the Fourier transform of p{r) = with m = 0.84 GeV, gives the result (2.28) 
for 


2.2 Unpolarized Structure Functions 

Probing the nucleon with a large wavelength photon (small momentum transfer squared 
—g^) can only provide information abont its dimension. It is possible to have a better 
spatial resolution by increasing —g^, but already for —g^ > the elastic process (2.1) 
is not dominant any more and the nucleon often breaks up into hadronic debris. The 
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inelastic reaction 

e(/) + iV(P)^e(/')+X(Px), (2.35) 

where X is the undetected hadronic system, represents the most direct way to explore 
the internal structure of the nucleon; it is fully inclusive with respect to the hadronic 
state and, when the momentum transfer is not ultra-high, is dominated by one-photon 
exchange, as described by the diagram in Figure 2.2. The reaction (2.35) is described by 
three kinematic variables. One of them, the incoming lepton energy P, or alternatively the 
center-of-mass energy squared s = (/ -|- P)^, is hxed by the experimental conditions. The 
other two independent variables can be chosen among the following invariants: = —q^ 

and z/, already dehned in (2.21) and (2.22), the center-of-mass energy squared of the 'y*N 
system (that is the invariant mass squared of the hadronic system X) 


= (P + g)^ 


(2.36) 


the Bjorken variable 


Xb 


and the “inelasticity” 


2P-q Q^ + W^-m‘^ ' 


P ■ q W‘^ + — m? 

P ■ I s — rn^ 


(2.37) 


(2.38) 


In the target rest frame, where v is the transferred energy from the electron to the target, y 
is the fraction of the incoming electron energy carried by the exchanged photon, y = u/E. 
The relation connecting xb, y and is given by 


xbV 


_ _ 

s — rn? s 


(2.39) 


Since > mf the Bjorken variable xb takes values between 0 and 1, and so does y. The 
measurement of the cross section corresponding to the process (2.35) shows a peak when 
the nucleon does not break up iW ~ m) and broader peaks when the target is excited to 
resonant barion states, most of them concentrated in the range m < hF <1.8 GeV. As 
W increases one reaches a region where a smooth behaviour is set in. This is the deep 
inelastic region, where both Q and W are large compared to the typical hadron masses. 
In this kinematical domain, the reaction (2.35) is known as deep inelastic scattering. 
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If we neglect the electron mass, the cross section for deep inelastic scattering (DIS) 
can be written as 


da 


1 1 
4P- / ^ 




d?l' 

(27r)32/' ’ 


(2.40) 


where, as before, we have averaged over the initial electron and nucleon spins, and summed 
over the hnal electron spin. The leptonic tensor /') was calculated in (2.17) and the 

hadronic tensor, corresponding to the electromagnetic transitions of the target nucleon to 
all possible hnal states, is dehned as 



Px){27r)V(l + P-1'- Fx), 


(2.41) 


with 


P£,(P-. Px) = I E I I I Ji(0) I N{P)} (2.42) 

spins 


and Px denoting the total four-momentum of the state |X). The dehnition (2.41) holds 
with states normalized as in (A.13). The cross section (2.40) reduces to (2.13), with 
me = 0, when X is restricted to be also a nucleon. 

In general, if we do not average over the initial nucleon spin S, the electromagnetic 
hadronic tensor Wafs consists of a symmetric and an antisymmetric (spin dependent) part 
under a ^ P, 


= ivy (g, P) + n%(q; P,S). (2.43) 

Since is hermitian, W*^ = hTsa, and (2.43) corresponds also to break Wafs into its 
real and immaginary parts. As is symmetric under a <-> /?, when contracted with 
Wap only the symmetric piece of Wap will contribute. Furthermore, the condition 

q'^Wap = Wapq^ = 0 (2.44) 


must hold for both real and imaginary parts of Wap due to current conservation, see (2.7). 
The most general form of compatible with parity conservation and with (2.44) is 


IF7^ = 




9ap + ^]wpP-q,q^) 




Pa - 


P ■ q 


Pp-^qp]W2{P-q,q^), (2.45) 
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where Wi^ 2 {P ■ are known as the structure functions of the nucleon. They are the 
generalization to the inelastic case of the elastic form factors. If we substitute (2.45) and 
(2.17) in (2.40), then 

= 4 (/ ■ l')Wi + 2[2{P ■1){P- /') - (/ ■ l')m^ ] ^ (2.46) 

47re^m 

and, in the nucleon rest frame, where P ■ I = mE, P ■ V = mE', one gets for the cross 
section (2.40) 


da 

dfidE' 


E'^^ ( 2Wi sin^ - + IT 2 cos^ 


(2.47) 


The cross section has again the characteristic angular dependence that was found for 
ep —> ep. One usually introduces the longitudinal and transverse structure functions 


Ft{xb-,Q‘^) — 2,xbEi{xb,Q‘^), (2.48) 

El{xb,Q'^) = E2{xb,Q‘^) — 2xbEi{xb,Q‘^), (2.49) 


which correspond to the absorption of transversely and longitudinally polarized virtual 
photons respectively, and in (2.48), (2.49) are expressed in terms of the following dimen¬ 
sionless structure functions 


E,{xb,Q^) = mW,{P-q,q^) (2.50) 

F2{xb,Q^) = vW2{F-q,q^). (2.51) 


Bjorken [37] argued that in the limit 


cx), xb =- hxed, 

2mu 


(2.52) 


now referred to as Bjorken limit, Ei and E 2 approximately scale, namely depend on xb 
only. This behaviour is already present for > 1 GeV^ and is surprisingly in contrast to 
the strong dependence, roughly as of the elastic form factors of the nucleon. On 
the other hand the elastic form factors of a pointlike particle like the muon are constants 
independent of Q^\ the ep ep scattering cross section is given by (2.23) with Fi = 1 
and E 2 = 0. Hence scaling seems to be an indication of scattering from charged pointlike 
constituents of the nucleon, the partons, and historically its observation [38-40] inspired 
the so-called parton model [37,41]. 


Moreover, structure function measurements show that El ^ F 2 , suggesting the spin- 
1/2 property of partons, since a (massless) spin-1/2 particle cannot absorb a longitudinally 
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polarized photon [42], In contrast, spin-0 (scalar) partons could not absorb transversely 
polarized photons and so we would have Fi = 0, i.e. Fl = -F2, in the Bjorken limit. 
Partons are nowadays identihed with the quarks of quantum chromodynamics (QCD). 


To conclude, the hadronic tensor Wap describes the unknown coupling of the virtual 
photon to the nucleon in terms of the structure functions, which can be extracted from 
experiments. However, the hadronic tensor can also be computed from models; in this 
case it is useful to develope a technique which allows one to extract the structure functions 
from a knowledge of Wafs- Dehning the projection operators [43] 


Vf = 


Vf = 


3P ■ q 
S'lre'^a 


_p«p/3 _ g 


af3 


pa pfS 


3^ 


afS 


with 


P ' Q 2 
a = — -h m , 


2x 


B 


and using (2.45), (2.50), (2.51), one can see that 


(2.53) 


(2.54) 


= Fl 


(2.55) 


and 


= F 2 . 


(2.56) 


2.3 Polarized Structure Functions 


Polarized DIS, involving the collision of a longitudinally polarized electron on a polarized 
(either longitudinally or transversely) nucleon, provides a different, but equally important 
insight into the structure of the nucleon. As mentioned in the previous section, if we do 
not average over the nucleon spin, the hadronic tensor will consist also of an antisymmetric 
part, Imposing (2.44), can be expressed in terms of the polarized structure 

functions Gi and G 2 as follows [44,45] 

H'4 = ie„g,,q>\^mS^G,(P ■ q, q'‘) + |(P ■ q)S<’ - (S ■ q)P^] 



(2.57) 
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where S is the covariant spin vector of the nucleon, whose essential properties are 


^.p = 0, 


(2.58) 


Clearly hh^ changes sign under reversal of the nucleon’s polarization. From the cross 
section formula (2.40), one notices that Gi and G 2 cannot be obtained from an experi¬ 
ment with just a polarized target. Both the electron and the nucleon must be polarized, 
otherwise the term drops out. Analogously to (2.43), the leptonic tensor has to 

be generalized to 

= + (2.59) 

where s is the spin four-vector of the electron, dehned such that s ■ / = 0, = —1, and 

is given in (2.17). The additional, antisymmetric, term can be calculated using the 
spin projector operator (A.20). If in (2.2) we make the replacement 

u{l )—^ ^ (1+ 7 V)m(/,s), (2.60) 

then we can again utilize (2.17), without the factor 1/2 alredy included in (2.60), to 
compute the leptonic tensor and perform the sum over the initial electron spins with help 
of the completeness relation. Having inserted the projection operator, only one of the two 
possible polarizations will contribute. We get 

Lap = ^ e^Tr[(l -h 7^) (/ -f me)'ya{J'+ me)^p] (2.61) 

and the term proportional to 7® will give 

= 2ie^meeappaS^{l - I'Y . (2.62) 


For a high energy [E 3> mg), longitudinally polarized electron, the spin vector is 

2 A 1 

5“ = —r, Ae = ±-, (2.63) 

TTlg 2 

and (2.62) becomes 


= 2ie^eappJ^q'', 


(2.64) 


where the helicity of the electron Ag has been fixed to be -1-1/2. The amplitude squared 
in the cross section (2.40) will have the form 


(2.65) 
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with 

= 4|l(/.9)(S.,)-,"(S./)]mGi 

+ h[ (S . ,)(P . /) - (P. ,)(S . 1 ) IG^I (2,66) 

and given in (2.46). The difference of cross sections with nucleons of opposite 

polarizations will single out only the antisymmetric part of the leptonic and hadronic 
tensors, namely the second term in (2.65). For a longitudinally polarized nucleon (that is 
polarized along the incoming electron direction), with the kinematics specihed in (2.20), 
the spin vector reads 

^=(0,0,0,!) (2.67) 

and the polarized cross section is given by 

dAa 1 d(T+ d(T_ 
dE'dfi “ 2 [ dE'dfi ~ dE'dfi _ 

= -‘^^[{E + E'cose)mG^-Q^G2l ( 2 - 68 ) 

with the subscripts ± meaning ±5*. Taking the sum instead of the difference in the hrst 
line of (2.68), one recovers the result (2.47) for the unpolarized cross section. 

Similarly to the unpolarized case, one introduces the structure functions 

= mVGi(P ■ g,g^), (2.69) 

92ixB,Q'^) = mu^G2{P ■ q,q‘^), (2.70) 


which are observed to approximately scale in the deep inelastic region, and the projectors 
[43] 


vf = [■^^^ + 2(P-Q)a;J^A + (Q- 


27re‘^b \ im? 


+ 2(P ■ q)xB S^ + {q- S)qx^Pr,e<^^^\ (2.71) 


b = —Am 


+ 2(P. - («■ S)T 


(2.72) 
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such that 


and 


VfW^p = g2 (2.73) 

VfW^fS = 91+92. (2.74) 


2.4 Parton Model 


In the parton model the nucleon is considered to be made of collinear, free constituents, 
each carrying a fraction ^ of the nucleon four-momentum: the quarks and antiquarks. Here 
we limit the discussion only to quarks, the extension to antiquarks being straightforward. 
The cross section of deep inelastic scattering is then described as the incoherent sum of 
all the electron-quark cross sections dd: 


da 


[ d^(l{^, + S)da 

q,s "'0 


(2.75) 


where s] S) is the number density of quarks q, with charge Cg in units of e, four- 
momentum fraction ^ and spin s inside a nucleon with spin S and four-momentum P. 
The cross section dd refers to the electron-quark scattering subprocess 


eil) + qm ^ ein + qik') 


(2.76) 


and, similarly to (2.15), (2.16), after integration over the struck quark phase space, reads 


dd 




d^Z' 

(2^^’ 


(2.77) 


where q = I —I', Lag is given in (2.17) and the quark tensor Wap = q, s) is the same 

as the leptonic tensor Lap, with the replacements /" —-|- q°‘. That is 


Wapi^, q, s) = q) + q, s), 


(2.78) 


with 


wlp{i,q) = 2e^[2i^PaPp + iPaqp + iqcPp-i{P ■q)9ap] 

wtp{i,q,s) = -2ie^mgeappaS^q'' 


(2.79) 




Parton Model 


and the quark mass for consistency is taken to be rriq = before and after the interaction 
with the virtual photon. By comparison of (2.75) with (2.40), and using the relation 

mp + qf) ^ + 2^^ ■ q) = 777^ - xb)^ (2.80) 

2F ■ q 

one can express the hadronic tensor Waf^ in terms of the quark tensor Wa/s as follows 

^Wa,(3{q-,P,S) = - XB)q{^,s]S)Wc,/3{^,q,s). (2.81) 

27r ^ ^ ZP ■ q Jq 4 

From (2.53)-(2.56) and (2.79)-(2.81), one obtains the parton model predictions for the 
unpolarized structure functions: 

PiixB) = (2.82) 

Q 

and 

P2{xb) = XB^^e^qi^XB) = 2 xbFi{xb)- 
<1 

where the unpolarized quark number densities q{xB) are dehned as 

q{xB) = ^q{xB,s-,S). 

S 

From (2.82) and (2.83), the Callan-Gross relation [46] follows 

Fl{xb) = F-2{xb) - 2xbFi{xb) = 0, (2.85) 

and it turns out that a mesurement of the structure function F 2 {xb) allows us to determine 
the momentum distributions of partons in the nucleon. 

From (2.71)-(2.74) and (2.79)-(2.81), the polarized nucleon structure functions are 
obtained: 

giixB) = ^^e^AqixB), ( 2 . 86 ) 

92ixB) = 0, (2.87) 

where 


(2.83) 


(2.84) 


Aq{xB) = q{xB,S] S) - q{xB, -S] S) 
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is the difference between the number densities of quarks with spin parallel to the nucleon 
(s = S) and those with spin anti-parallel (s = —S). Fixing the nucleon to be longitudi¬ 
nally polarized with positive helicity, (2.84) and (2.88) can be rewritten in terms of parton 
densities with definite helicity, with notation analogous to (1.6) and (1.7): 

(1{xb) = gtixs) + Aq{xB) = qt{xB) - q-{xB), (2.89) 

where Aq{xB) measures how much the parton q “remembers” of its parent nucleon po¬ 
larization. 

The structure function QiIxb) yields information on how the helicity of the nucleon is 
distributed among its parton constituents, while q-z^xb) has not a simple interpretation 
in the parton model. It can be shown [43] that, if we allow the partons to have some 
transverse momentum fcj_ inside the nucleon, then g 2 {xB) is non-zero. However, it cannot 
be calculated without making some model of the k± distribution. 

2.5 Structure Functions in QCD 

The parton model is only the zero-th order approximation to the real world; quarks 
and antiquarks are not free particles, they interact by emitting and absorbing gluons. A 
detailed discussion of QCD, the theory which describes the strong intractions of quarks 
and gluons, can be found in [42,47,48]. 

From an empirical point of view, one observes that the scaling predicted by the parton 
model is violated. Structure functions appear to depend on Q^, although in a relatively 
mild way, logarithmically. This behaviour arises from perturbative QCD and represents 
the original, and still one of the most powerful, quantitative test of the theory. The 
radiation of gluons produces the Q^-evolution of the quark (and antiquark) distributions 
in (2.89), furthermore it determines the appearence of the unpolarized and polarized gluon 
distributions, defined in a way similar to (2.89), 

9{xb,Q^) = gl{xB,Q‘^) + g-{xB,Q^), 

Ag{xB,Q^) = gXixB.Q"^) - gt{xB,Q'^)] (2.90) 


g^{xBiQ‘^) and g^{xB,Q^) being the densities associated to the positive and negative 
circular polarization states of the massless, spin-1 gluon. Moreover, Similarly to QED, it 




2.5 Structure Functions in QCD 


25 


is possible to define an effective “fine-structure constant” for QCD, 


a, = 


47r ’ 


(2.91) 


with Qs being the strong coupling. Furthermore it is convenient to introduce the dimen¬ 
sional parameter A, because it provides a description of the dependence of on the 
renormalization scale (in DIS usually identified with the scale of the probe Q). The defi¬ 
nition of A is arbitrary; one possibility is to write as as an expansion in inverse powers of 
InQVA^, 


4% 


(3i InlnQVA^ 


+ 0 


(2.92) 


/3olnQVA2 PI (lngVA2)2 Vln'(gVA2) 

where Pq = 11 —2n//3, Pi = 102 —38n//3 and Uf is the number of quarks with mass less 
than the momentum scale Q. Equation (2.92) illustrates the asymptotic freedom pTopeitj: 
CKs —0 as cx) and shows that QCD becomes strongly coupled at g ~ A. Therefore 

perturbative calculations (and the parton model) are reliable only for large momentum 
transfer. The value of A depends on the renormalization scheme adopted and must be 
determined from experiment. 

Once the parton distributions are fixed at a specific input scale = Qq, mainly by 
experiment, their evolution to any > Ql is predicted by perturbative QCD. If we 
define, in the unpolarized sector, the flavor nonsinglet distributions 


(Jns- = u — u, d — d, 

(3 'ns+ = (u + u) - {d + d), {u + u) + {d +d) - 2{s + s), 


(2.93) 


and the singlet combination 


s = (g + g), 

q=u,d,s 


at NLO QCD the evolution equations take the form 


(2.94) 


and 


^ g^) = Tns± C) gxsit , 


^ / S(a;B,gQ \p / S 

V 9{xb,Q^) 


(2.95) 


(2.96) 
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Figure 2.3: Measurements of the proton structure function F 2 {xb, Q^) from HERA colla¬ 
borations and fixed target experiment. The data are shown as a function of for various 
fixed values of xb (denoted as x in the plot). Curves represent a NLO QCD fit to the 
data. 


where t = IuQ^/Qq and the convolution (®) is dehned by 

{P ^q){xB,Q‘^) = [ — p(—'] q{y,Q‘^). 

Jxs y \y J 


The splitting functions are given by 


-Pns± — 


asm 

2ti 






and 


P 


asm 

27r 


P^°\xb) + 



(2.97) 


(2.98) 


(2.99) 
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Xb 

Figure 2.4: Comparison of two NLO QCD fits, corresponding to different input distribu¬ 
tions [8], with the measurements of the proton structure function F 2 at small xb [62]. The 
results and data for the bins in = 1.5 CeV^ and 2 CeV^ have been multiplied by 0.75 
and 0.85, respectively, as indicated. 

where 

with, = 0,1; (2.100) 

V Pm P 9 9 J 

being given in (2.92). The LO expressions are entailed in (2.92) and in (2.98)- 
(2.99); they can be obtained by simply dropping all higher order terms {fii, 

At LO (2.95), (2.96) reduce to the well-known DGLAP evolution equations [49-52]. The 
complete set of NLO splitting functions in the commonly used MS factorization scheme 
has been calculated [53-55] and can be found, for example, in [48]. 

Equations (2.95)-(2.99) hold also in the polarized sector, with the formal replacements 
q Aq and P AP. The polarized splitting functions are known up to NLO in the 
MS scheme [56-58] and are listed in [59]. 
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Fl extraction from HI data (for fixed W=276 GeV) 



Figure 2.5: Measurements of the proton structure function Fl by the HI collaboration at 
HERA at a fixed value of W = 276 GeV [65]. 


The resulting NLO parton distributions are directly related to physical quantities, 
such as structure functions, by a convolution with calculable, process dependent, coeffi¬ 
cient functions. For consistency, the choice of the factorization convention must be the 
same for both the coefficient functions and the splitting functions underlying the parton 
distributions. Within the MS scheme, F 2 {xb,Q^) is given by 


1 

Xb 


F2{xb,Q^) 


^ I q{xB, + qixB, Q^) 

n ^ 


as{Q^) 

271 


Cg {q + q) + 2Cg ® g 


( 2 . 101 ) 


where the unpolarized coefficients functions Cq^g are given, for example, in [60] and in 
[61]. Measurements of the proton structure function F 2 together with a NLO QCD fit, 
presented in [62], are shown in Figure 2.3. As one can see, QCD appears to predict 
correctly the dependence of structure function over four orders of magnitude. 

Recently, a dedicated test of the validity of (2.101) at very low xb in the perturbative 
regime, > 1 GeV^, has been performed [8] and the results are shown in Figure 2.4. A 
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X 


Figure 2.6: Measurements of the nucleon structure function gi as a function of xb (denoted 
as X in the plot) at = 5 GeV^ [67]. The data are from the following experiments: El55 
(solid circles), E143 (open circles), SMC (squares), HERMES (stars) and E154 (crosses). 
The solid curves correspond to a NLO QCD fit, while the dashed curves are from a purely 
phenomenological fit. 


good agreement with recent precision data for F 2 [63], restricted to 

3 X 10“^ < < 3 X 10-^ 1.5 GeV^ < < 12 GeV^ (2.102) 

has been found, as well as with the present experimental determination of the curvature 
of F 2 [64]. 

The structure function Fl{xb, Q^) is only non-zero at order Ug in perturbation theory, 
i.e. Fb^xbjQ"^) = 0{as)', deviations from the Gallan-Gross relation (2.85) are evident in 
Figure 2.5, taken from [65]. The data are much less precise than the ones for F 2 {xb, Q‘^), 
but QGD seems to work well also in this case. 
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Similarly to F 2 {xb,Q‘^), gi{xB,Q^) can be written as 


gi{xB,Q^) = ellAq{xB,Q^) + M{xb,Q^) 


+ 


2tt 


ACg ® (Ag + Aq) + 2ACg ® Ag 


(2.103) 


where ACq^g can be found in [66]. In Figure 2.6 results on gi presented in [67] are shown 
as a function of at = 5 GeV^. 


At LO the coefficients {A)Cq^g vanish; from (2.101) and (2.103) it turns out that 
the gluon does not directly contribute to the structure functions, but only indirectly via 
the Q^-evolution equations. The sums in (2.101) and in (2.103) usually run over the light 
quark-flavors q = u, d, s, since the heavy quark contributions (c, b, ...) have preferrably to 
be calculated perturbatively from the intrisic light quarks {u, d, s) and gluon (g) partonic 
constituents of the nucleon. This treatment of the heavy quarks underlies the unpolarized 
GRV98 [68] and polarized GRSVOl [69] parton distributions, which will be used in our 
numerical estimates in the perturbative regime. For our studies in the low-Q^ region, 
where perturbative QGD is not applicable, we resort to the purely phenomenological 
parametrization ALLM97 of F 2 [70] and to the parametrization BKZ of gi [71], shortly 
described in Ghapters 5 and 7 respectively. 


To conclude, it has been shown that the inclusion of 0(a) QED corrections to the 
parton evolution modihes only very slightly equations (2.95)-(2.96) [72-74], so we will 
not consider such effect. We will concentrate on the other consequence [73,74] of the 
emission of photons from quarks, the appearence of the (inelastic) photon distributions 
of the proton and the neutron. 




Chapter 3 

The Equivalent Photon Distributions 
of the Nucleon 


In this chapter the polarized and nnpolarized photon content of protons and neutrons, 
evaluated in the equivalent photon approximation, are presented. In particular, the uni¬ 
versal and process independent elastic photon components turn out to be uniquely deter¬ 
mined by the well-known electromagnetic form factors of the nucleon and their derivation 
is shown in detail. The inelastic photon components are obtained from the corresponding 
momentum evolution equations subject to the boundary conditions of their vanishing at 
some low momentum scale. The resulting photon asymmetries, important for estimating 
cross section asymmetries in photon-induced subprocesses are also presented for some 
typical relevant momentum scales. 

3.1 Unpolarized Photon Distributions 

As already mentioned in the Introduction, the concept of the photon content of (charged) 
fermions is based on the EPA, the equivalent photon approximation. Applied to the 
nucleon N = p, n it consists of two parts, an elastic one due to A^ —yA^ and an inelastic 
part due to N ^ yX with X ^ N. Accordingly the total photon distribution of the 
nucleon is given by 

y(a;, /i^) = yei(a;) yinei(a:, /i^), (3.1) 

where x is the fraction of the nucleon energy carried by the photon and /x is a momentum 
scale of the photon-induced subprocess. The two components are discussed separately in 
the following. 
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3.1.1 Elastic Component 

The elastic photon distribution of the proton, 7 g[, has been presented in [24] and can be 
generally written as 

2 2 

_ m X i/'r/ \ on \ 

where 

H,(t) ^ 7(0 + r7(0 = + rGi,(«) ,3 3 , 

1 + r 

with r = —t/Arn?, m being the nucleon mass, and where Ge and Gm are the Sachs elastic 
form factors discussed in Section 2.1. 

The result (3.2) can be obtained extending to an unpolarized nucleon N the analysis 
[21] for a photon emitting unpolarized electron. We consider the process 

N{P) + a{l)^ N{P')+X, (3.4) 

where the target a is a massless parton {P = 0), P^ = P'^ = m?, and X is a generic 
hadronic system. The corresponding cross section can be written as 

P')W<.dk, l) . ( 3 . 5 ) 

where 

q = -k = P' -P, (3.6) 

k being the four-momentum of photon emitted by the nucleon. The tensor H^f{P-, P') is 
dehned in (2.18) and, using (2.24), (3.3) and (3.6), can be rewritten in a more compact 
form as 

Hf{P- P') = [H^{q^){2P - A;)“(2P - kf + Gl,{q^){q^- rA;^)]. (3.7) 

The partonic tensor Wag{k,l), as shown in (2.45), can be decomposed as 

Waq{k,l) = (^-gap + ^ k^kg^ Fi{q^, k ■ 1) 

~Jk^ ^ ^ ^^{q ,k-1) , (3.8) 
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with the structure functions Fi and F 2 dehned in (2.50) and (2.51). As pointed out in 
[21], in the limit 0, Wap{k,l) must be an analytic function of therefore q^Waq 

has to vanish for = 0. This implies 

{k-l) 


r 


■F2{q\k-l)=F^{t),k-l)+0{q^), 


(3.9) 


and the terms 0{q^) will not be considered in the following. Furthermore, one can intro¬ 
duce the variable 

a; = ^, (3.10) 

which represents the fraction of longitudinal momentum carried by the emitted photon, 
assumed to be real and collinear with the parent nucleon, that is 


k = xP. 


(3.11) 


From (3.7)-(3.10) one hnds 


Hf{P- P')W^p{k, 1) = 47re^<' AH,{q 


m H-^ 

X x^ 


2q^GUqn'>F,{0,Fk). 


If we write the four-momenta of the incoming and outgoing nucleons as 
P = (E, 0,0, E/3), P' = (P',O,P'/3'sin0,P'/3'cos0), 
with ^_ 

r m3 


p = \n- 


/ 5 ' = \ 1 - 


P2 ’ ^ V P'2 ’ 

then the phase space for the scattered nucleon will be given by 




(3.12) 


(3.13) 


(3.14) 


(3.15) 


= 2nP E dP d cos 9, 

where the azimuth integration has already been carried out. The integration variables 
(P', cos 6*) can be replaced by (g^, x)] using the following relations 


g2 = 2m^ -2EE\1- f3f3'cos 9), 
EUl+/3^cos9) 

X = 1 — 


(3.16) 

it can be proved that the Jacobian of this change of variables is 2P'/3'2 and (3.15) becomes 


d^P' 


— = Trdg dx. 


(3.18) 
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The cross section (3.5) can now be written as 


daNa{P, 0 = - 


where t = and 


a 1 f r 2m4x ^/1 


2ti t \ \ t 


+ 2^— — !^ Hi(t) + a^a{k,l) dtdx, (3.19) 


1 TTP^ 

<r-,a(k,l) = = —Fi{0,k-l) 


(3.20) 


is the cross section for the process 'y{k) + N{P) X for a real photon. Integrating over 
t one gets 

dajva(P, /) = 7 ei(a;) a^a{k, 1) dx, (3.21) 


where 7 ei(a;) is given by 


7 ei(a;) = -- 


2 --1 + 


2m?'x 


Hi{t) + xG\i{t) 


(3.22) 


The extrema of integration can be determined from (3.16) and (3.17); from the latter we 


1 — cos^ 6 


1- a/1- 


{A^ + w? cos^ 6*)(1 — cos^ 9) 


where 


A = E{lp(3){l-x). 


Expanding (3.23) in powers of 6*, for 6* -C 1, one gets [21] 


(A^ — 


E' = ^ + 0(9% 


and (3.16) becomes 


2 m^x^ E{1 + P){A^ — m^f 2 


q = - 


1 — X 


r + 0{9^ 


The value tmax = Q'max is obtained by taking 0 = 0, namely 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


^max 


2 2 
m X 

1 — X 


(3.27) 


The minimum value of the photon virtuality tmin can be computed in a simple way by 
imposing that the invariant mass W'^ of the produced hadronic system X be bounded 
from below [21]: 

= (k + ly > ( 3 . 28 ) 
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from which it comes, for hh^in = 

train = —‘^k ■ I + = —2xP ■ I + m? = —x{s — fu'^) + 771^, (3.29) 

where s is the nucleon-parton centre-of-mass energy squared. 

The Sachs form factors which appear in (3.22) are conveniently parametrized by the 
dipole form proportional to (1 —1/0.71 GeV^)“^ as extracted from experiment, see (2.28). 
As pointed out in [24], this implies that the support from values t —0.71 GeV^ to 
the integral in (3.22) is suppressed, hence, in the kinematical region xs 3> one may 
integrate from tmin = —oo to tmnx = —nUx‘^/{l — x) so as to obtain the universal process 
independent 7 ei in (3.2) . 

Equation (3.2) can now be analytically integrated, remembering that, from (2.28) and 
(3.3), for the proton we have 

= (1 + «^)“^ > Hl{t) = + ar)-^ (3.30) 

with /Xp = 1 + Kp ~ 2.79 and a = 4m^/0.71 GeV^ ~ 4.96, while for the neutron 

G^it) = KnT{l + aT)~^ , Gl^{t) = Kn{l + aT)~^ , i7”(t) = K^r(l + ot)”"^ , (3.31) 

with Kn — —1.79. After integration, the elastic component of the equivalent photon 
distribution reads [1], for the proton 

(3.32) 
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with a_ = a — 1. For arriving at (3.32) we have utilized the relation 


dr- 


r^fl + arF 

4(l-x) ^ ' 


= -Aal + 4- 


X 


x'^z^ ’ 


which will be also relevant for the polarized photon contents to be presented below. The 
result in (3.32) agrees with the one presented in a somewhat different form in [24]. 

If we integrate (3.21) taking = 1 (pointlike particle), with the integration 

bounds (3.27) and fmin = —we recover the result (1.4) for the unpolarized photon 
content of the electron. 


3.1.2 Inelastic Component 


As pointed out in [27], the complete function 7 ( 0 ;, /i^) in (3.1) could be built-up by adding 
to the elastic contribution all resonant [34] and non-resonant hnal hadronic states, and 
their interferences. Alternatively one could guess an inclusive or ’continuous’ 7 inei(a^, At^), 
based on the parton model, where the photon is emitted by one of the quarks in the 
nucleon [25]. In this latter picture, which we adopt, the addition of the resonant and 
continuous contributions may be double-counting. The inelastic part in (3.1) is then 
given by the leading order (LO) QED evolution equation [25] 


d7inei(a^,h^ 

dlUyU^ 




q=u,d,s 


dy 


X 


y 


where 


Pzqiy) = 


y 


(3.36) 


(3.37) 


is the quark-to-photon splitting function and q{y,iji^)i Q'(|/,/i^) are respectively the quark 

(-) (-) (-) (-) 

and antiquark distribution functions of the nucleon at LO QCD [ 68 ], with u^=dP ', (P‘=u^, 
(-) (-) 

s^=s^. Equation (3.36), which states that the probability to hnd a photon in the nucleon 
is given by the convolution of the probabilities to hnd hrst a quark inside the nucleon and 
then a photon inside the quark, is integrated subject to the ‘minimal’ boundary condition 


7inei(a;,/io) = 0 (3-38) 

at [ 68 ] Hq = 0.26 GeV^. The boundary condition (3.38) is obviously not compelling and 
affords further theoretical and experimental studies. Since for the time being there are 
no experimental measurements available, the ‘minimal’ boundary condition provides at 
present a rough estimate for the inelastic component at /i^ 3> /ig- 
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3.2 Polarized Photon Distributions 

The previous analysis can be extended to the polarized sector, i.e., to 

= A7ei(a;) + Ayine^a;, /i^). (3.39) 

As before, the elastic and inelastic parts will be stndied separately. 


3.2.1 Elastic Component 

The elastic part A 7 ei(a;) in (3.39) is determined via the antisymmetric part of the tensor 
describing the photon emitting nncleon N 

Hf{P- P') = ^e^Tr [(1 + 7^^)(f+ m)r“(f' + m)r^] (3.40) 

for the process 

N(P; S) + a(l; s) ^ N(P') + X (3.41) 

where a being a parton with four-momentnm I initially kept off-shell and S, s are the 
polarization vectors [ 22 ] satisfying the transversality condition S ■ P = 0 and s ■ I = 0. 
Eqnation tensor has been obtained in a way similar to (2.61), with 7 “ —T". In terms of 
the Dirac and Panli form factors Fi, 2 {t) the elastic vertices T" are given by 

r“ = (Fi + F2)7“-^F 2 (P + P')“- (3.42) 

2m 

The analysis has been carried ont originally in [1] and is a straightforward extension of the 
calculation [ 22 ] of the polarized eqnivalent photon distribntion resnlting from a photon 
emitting electron. The antisymmetric part of (3.40), as calculated in [1], reads 

Hf^(p-P') = 2te^mGl,e’^^P^SpK + 2tGM{F2/2m) 

X (P + P')“ -{P + SpP^P',, (3.43) 

with, as before, k = P — P'. One can show that 

(P + p'Yehp--' - (P + SpP,P'^, = {k ■ S)P^k^>e^^^^' 

+ {2m^ + P ■ k)e^P^^'Spk^>, (3.44) 

where we made use of the e-identity 

gafi^Ppaa' _ ^afS^ppaa' _|_ gap^/Spaa' _|_ ^fippa' _|_ ^aa'^jSpcrp. 


(3.45) 
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hence (3.43) can be rewritten in a more compact form [75] as 

P') = UgeGmS^ - ^M{GM -GE)k-S ^ 

-L ~1~ ^ Ttlj 

The polarized cross section relative to the process (3.41) is given by 

0 = ^ H:r(P; E) Wi^(k, 1) , ( 3 , 47 ) 

where is the antisymmetric part of the partonic tensor, describing the polarized 

target a(/; s), which is expressed in terms of the usual polarized structure functions gi 
and g 2 , see (2.57) together with (2.69) and (2.70), 

1) = e^f,p,kP g^{q\ k ■ l)s^ + g^^q^ ^ ' 0 ^ • (3-48) 

The spin vectors for the incoming nucleon and parton can be written as 

S^ = Ns{Pc.-^l}j, = (3.49) 


where the normalization factors Ns and Ng are related to S'^ and by [76] 


Nl = 




(3.50) 


(3 = Gl- 


(3.51) 


y (p-/) 2 ’ 

and are hxed in order to satisfy the condition IP^I = |s^| = 1. 

Putting the parton on-shell {p = 0) and using the dehnition (3.10), together with (3.49), 


we have 




2m?x 


1 H— G^ 

X 


^ _ X + ^ Gm{Gm — Ge) 

q^ X 1 + T 


gi{q\k-l), (3.52) 


where all the terms proportional to g 2 drop from this equation. Then (3.47) becomes 


= -IT- 


2m^x 21 2 

i H 2 

q^ X 


2 1 + 


rri^x ll Gm{Gm ~ Ge) 


q^ X 


1 + r 


2^9i{0,kG) 

^ — P~l — ’ 
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which holds in the limit 0, after changing the variables of integration as for the 

unpolarized case. This expression can be related to the polarized cross section for real 
photon-parton scattering, which can be computed by convoluting the partonic tensor with 
the antisymmetric part of the photon polarization density matrix [76] 

= - epel) = (3.54) 

where is the photon polarization vector and is its spin vector 

ta = Nt ^ , (3.55) 

with Nt chosen so that \t‘^\ = 1. We get (g^ —0) 

l) = jN ^ . ( 3 . 56 ) 

Combining the last equation with (3.53) and integrating over g^, one gets the analogous 
of (3.21) for a polarized process 

dAaAra(T’, 1) = A7ei(a;) Aa^a(/c, l)dx, (3.57) 

with [1] 

A7ei(a;) = -^ / yS 2-x + — - — G^(t) - 2 1 - x +—^ GM(t)F 2 (t)> 
^^min V L J L J J 

= -^ jGM(t)\ 2-x + —^ F,(t)+xF 2 (t)'> (3.58) 

^imin V L J J 

where the hrst term proportional to in the hrst line corresponds to the pointlike result 
of [22]. Following [24], we again approximate the integration bounds by fmin = —oo and 
^max = —m^a;^/(l — x) as in (3.2) in order to obtain an universal process independent 
polarized elastic distribution. Using, in addition to (3.30) and (3.31), 


Ff(t) 

1 + fXpT 2 

= (1 + ar) , 

F|(t) = ^(l + ar)-^ 

1 + r 

(3.59) 

F-(t) 

= 2Kn (1 + ar) 

1 + r 

F^{t) = Kn (1 + ar) ^ 

1 + r 

(3.60) 


equation (3.58) yields for the proton 

^ 7 ^ 1 ( 2 ^) = { (2-3^) + +2aa;^ I + 2Kp (^1 - x + / - 2^^! , 

(3.61) 
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-4 -3 -2 -1 

10 10 10 10 X 

Figure 3.1: The polarized and unpolarized total photon contents of the proton, A'y^ and 
7 ^, according to (3.1) and (3.39) at some typical fixed values of (in GeV^). The 
-independent elastic contributions are given by (3.61) and (3.32). 


and for the neutron 

|a^(l - x)I + 4 (^1 - X + (3.62) 

with I and I being given in (3.34) and (3.35). 

Integrating (3.58) between tmin = —/i^ and tmax given in (3.27), with the form factors 
Gm = Fi = —1 and F 2 = 0, we obtain the polarized photon distribution of the electron 
( 1 . 10 ). 


3.2.2 Inelastic Component 


The inelastic contribution [1] derives from a straightforward extension of (3.36), 
dA'jinefix, p‘^) a ^ ^2 


din p3 


a 


q=u,d,s 


■ ^^'>9 ( ~ ) [^9{y, + Ag(|/, p^)] (3.63) 
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Figure 3.2: As in Figure 3.1 but for a linear x scale. 

where 

1 — (1 — 

AP,,{y) = -^^ = 2 - y (3.64) 

y 

is the polarized quark-to-photon splitting function and Aq{y, /i^), Aq{y, jjf) are the polari¬ 
zed quark and antiquark distribution functions of the nucleon. We integrate this evolution 
equation assuming again the not necessarily compelling ‘minimal’ boundary condition 

A7i„ei(a;,Po) = 0, (3.65) 

according to |A 7 inei(x,/Xg)! < jjff) = 0, at /Xg = 0.26 GeV^ using the LO polarized 

parton densities of [69]. These latter two equations together with (3.63) yield now the 
total photon content A'y{x,y‘^) of a polarized nucleon in (3.39). 

3.3 Numerical Results 

Our results for A'yP{x,y‘^) in (3.39) are shown in Figure 3.1 for some typical values of 
/i^ up to yf = = 6467 GeV^. For comparison the expectations for the unpolarized 
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Figure 3.3: The asymmetry of the polarized to the unpolarized photon content of the 
proton as defined in (3.66) at various fixed values of pf (in GeV^) according to the results 
in Figure 3.1. The -dependence of the elastic contribution to is caused by the p^- 
dependent total unpolarized photon content in the denominator of (3.66). For illustration 
the p'^-independent elastic ratio is shown as well. 


in (3.1) are depicted as well. The /i^-independent polarized and unpolarized 
elastic contributions in (3.61) and (3.32), respectively, are also shown separately. Due to 

the singular small-y behavior of the unpolarized parton distributions y q {y,p^) in (3.36) 
as well as of the singular in (3.32) as a; —^ 0, the total x'yP{x,p‘^) in Figure 3.1 

increases as x —> 0, whereas the polarized xA'y^{x, p^) —>• 0 as x —^ 0 because of the 
vanishing of the polarized parton distributions yA q {y, p^) in (3.63) at small y and of the 
vanishing xA'j^^^x) in (3.61) at small x. In fact, xA'jP^x, p"^) is negligibly small for x ^ 
10 “^ as compared to xjP{x,p'^). 

For larger values of x, x > 10“^, xA 7 ^(x,/i^) becomes sizeable and in particular is 
dominated by the /^^-independent elastic contribution xA'j^fx) at moderate values of p‘^, 
p^ A 100 GeV^ (with a similar behavior in the unpolarized sector). This is evident from 
Figure 3.2 where the results of Figure 3.1 are plotted versus a linear x scale. 
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Figure 3.4: As Figure 3.1 hut for the neutron, with elastic polarized and unpolarized 
contributions being given by (3.62) and (3.33). 


The asymmetry Ap{x,^‘^) is shown in Figure 3.3 where 

A^{x, pf) = [A7ei(a:) + Ayine^x, pf)] /'^{x, p^) (3.66) 

with the total unpolarized photon content of the nucleon being given by (3.1). To illustrate 
the size of Ay^j relative to the unpolarized we also show the /i^-independent ratio 
Ayei(3^)/7ei(^) Figure 3.3 which approaches 1 as x —>• 1. 

The polarized photon distributions Ay^(x,/r^) shown thus far always refer to the so 
called ‘valence’ scenario [69] where the polarized parton distributions in (3.63) have flavor- 
broken light sea components Au ^ Ad ^ As, as is the case (as well as experimentally 
required) for the unpolarized ones in (3.36) where u ^ d ^ s. Using instead the somehow 
unrealistic ‘standard’ scenario [69] for the polarized parton distributions with a flavor- 
unbroken sea component Au = Ad = As, all results shown in Figures 3.1-3.3 remain 
practically almost undistinguishable. The same holds true for the photon content of a 
polarized neutron to which we now turn. 
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Figure 3.5: As in Figure 3.4 but for a linear x scale. 


The results for A'y‘^{x, gif) are shown in Figure 3.4 which are sizeably smaller than the 
ones for the photon in Figure 3.1 and, furthermore, the elastic contribution is dominant 
while the inelastic ones become marginal at x 0.2. For comparison the unpolarized 
7 "(x, jjf) in (3.1) is shown in Figure 3.4 as well. Here, 7 ”) in (3.33) is marginal and xj'ffx) 
is non-singular as x —>• 0 with a limiting value X72(x)/a = {fna) ~ 0.078. Thus the 
increase of X 7 "'(x,/x^) at small x is entirely caused by inelastic component xjlf^fx, in 
(3.36), due to the singular small-?/ behavior of ?/ <? (?/, /i^), which is in contrast to X 7 ^(x, /x^) 
in Figure 3.1. 

These facts are more clearly displayed in Figure 3.5 where the results of Figure 3.4 are 
presented for a linear x scale. Notice that again the polarized xA 7 "'(x,/x^) —> 0 as x —> 0 
because of the vanishing of the polarized parton distributions yA q (?/,/x^) in (3.63) at 
small y and of the vanishing of xA 7 gj(x) in (3.62) at small x. 

Finally, the asymmetry H!)(x,/x^) defined in (3.66) is shown in Figure 3.6 which is 
entirely dominated by the elastic contribution for x ^ 0 . 2 . As in Figure 3.3 we illustrate 
the size of the elastic A'y^fx) relative to the unpolarized 7 ei(x) by showing the ratio 
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Figure 3.6: As Figure 3.3 but for the neutron asymmetry according to the results in Figure 

3.4. 


^7d/7ei Figure 3.6 as well. Notice that A 7 "[/ 7 gi —| as a; ^ 1 in contrast to the case 
of the proton. 

Clearly, the nucleon’s photon content 7 ( 0 ;,/i^) is not such a fundamental quantity 
as are its underlying parton distributions /(x, = q, q, g or the parton distributions 

/'>'(x, of the photon, since 7 ^(x, is being derived from these more fundamental 
quantities. Moreover, its reliability remains to be studied. We shall try to carry out this 
task starting on Chapter 5. 

A FORTRAN package (grids) containing our results for as well as those for 

7 (x, can be obtained by electronic mail. 











Chapter 4 

Measurement of the Equivalent 
Photon Distributions 


In the previous chapter we estimated the polarized and unpolarized equivalent photon 
distributions of the nucleon (A) 7 (a;, consisting of two components, 

(A)7(a;,/i^) = (A)7ei(a;) + (A)7inei(a;, (4.1) 


where the elastic parts (A) 7 ei are uniquely determined by the well-known electromagnetic 
form factors ^ 1 ^ 2 ( 5 '^) of the nucleon. The inelastic components are hxed via the boundary 
conditions 


(A)7inel(3^, Po) 0 


(4.2) 


at pg = 0.26 GeV^, evolved for > pg, according to the LO equations 
d(A)7inei(a;,p2) a 2 T dy 


din ij? 


It Z f A)9(!/.y) + (A)9(!/.7)]. (4.3) 


with the unpolarized and polarized parton distributions in LO taken from [68,69]. As 
already stated, the boundary conditions are not compelling but should be tested expe¬ 
rimentally. However at large scales the results become rather insensitive to details at 
the input scale pg and thus the vanishing boundary conditions yield reasonable results for 
(A)7inei which are essentially determined by the quark and antiquark (sea) distributions 
of the nucleon in (4.3). At low scales however, (A) 7 inei(a;, p^) depend obviously on the 
assumed details at the input scale pg. Such a situation is encountered at a hxed target 
experiment, typically HERMES at DESY. At present it would be too speculative and 
arbitrary to study the effects due to a non-vanishing boundary 7inei(a^, ho) 7^ 0- Rather this 
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should be examined experimentally if our expectations based on the vanishing boundary 
turn out to be in disagreement with observations. 

In the present chapter we consider two processes which offer a clear opportunity to 
gain information on the photonic structure of the nucleon: muon pair production in 
electron-nucleon collisions eN efj,~^fj,~X via the subprocess and the 

QED Compton process eN —eyX via the subprocess ey for both the HERA 

collider experiments and the polarized and unpolarized fixed target HERMES experiment 
at DESY. The production rates of lepton-photon and dimuon pairs are evaluated in the 
leading order equivalent photon approximation and it is shown that they are sufficient 
to facilitate the extraction of the polarized and unpolarized photon distributions of the 
nucleon in the available kinematical regions [2], On the other hand, it should be noted that 
a study of NN via y^y^ —in hadron-hadron collisions is impossible 

[29, 30] due to the dominance of the Drell-Yan subprocess . The logarithmic 

enhancement of the photon densities is not enough to overcome completely the extra factor 

in the yy fusion process. 

Measurements of (A)y(a;,/i^) are not only interesting on their own, but may provide 
additional and independent informations concerning (A)^q^ in (4.3), in particular about 
the polarized parton distributions which are not well determined at present. 


4.1 Theoretical Framework 

In this section we present the kinematics and cross section formulae of the reactions 
under study, following the lines of Appendix D in [77]. The unpolarized and polarized 
subprocess cross sections for yy —and ey —> ey can be found, for example, in [78]; 
their derivation is shown in detail in Appendix B. 


4.1.1 Dimuon Production 

We consider first deep inelastic dimuon production ep —> ep,'^fi~X via the subprocess 




(4.4) 
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e 


/i+ 


Figure 4.1: Lowest-order Feynman diagram for dimuon production in ep collisions. The 
crossed u-channel diagram is not shown. 

as depicted in Figure 4.1. With the four-momenta of the particles given in the brackets 
in (4.4), the Mandelstam variables are dehned as 

s = {ki + /c 2 )^, i = {ki - /i)^, M = (^2 - ^i)^- (4.5) 

Suppose that the photon 7 ^^ carries a fraction f of the electron’s momentum and that a 
similar dehnition for x exists for the photon 7 ^. Then in the e — p center-of-mass system 
the four-momenta ki and ^2 of the colliding photons, assumed to be collinear with the 
parent particles, can be written as 

h = ^ (1, 0 , 0 ,1), h = ^ (1, 0 , 0 , -1), (4.6) 

where the positive z axis is taken to be along the direction of the incident electron and s 
is the squared center-of-mass energy, which satishes the condition 

s = fxs. (4.7) 

The four-momenta of the outgoing muons can be written in terms of their rapidities yi ^2 
and momentum components Iti 2 ^ transverse with respect to the 2 ; axis, 

= ^Ti (cosh ?/i, 1,0, sinhyi), ^2 = ^T 2 (coshj/ 2 ,-1, 0, sinh j/ 2 ). (4.8) 
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In general, an outgoing particle with energy E has component of the velocity along the 2 : 
axis given by 


/9 = 


4 

E 


and its rapidity y can be dehned so that 


E = /7’COsh|/, lz = lTsaih.y, 


(4.9) 


(4.10) 


where Iz and It are respectively the longitudinal and transverse components of its mo¬ 
mentum 1. Therefore, the relation between y and (3 is given by 


l/ = arctan/3 = ^In^^ 

^ -L 

and, substituting (4.9) in (4.11), one has also 

1 , E + lz 

!/ = X In 


(4.11) 


(4.12) 


2 E -I, 

For a massless particle, Iz = EcosO, 6 being the center-of-mass scattering angle, and 
(4.12) assumes the much simpler form 


, e 

y = — In tan - , 
y 2 


(4.13) 


called pseudorapidity and very convenient experimentally, since one needs to measure 
only 6 in order to determine it. 

The cross section relative to the inclusive production of the two muons, in the most 
differential form, is given by 

da 


d.^da; 




2s 


s) s) |M|2 {2TTE6‘^{ki + k2-li-l2 


d% 


dH. 


(27r)32Ei (27r)32E2 ’ 

(4.14) 


where d is the cross section of the subprocess 7 ^ 7 ?’ — yEyr. In the spirit of the leading 
order equivalent photon approximation underlying (4.14), we shall adopt the LO photon 
distribution of the proton 7 ^’(a;, s) given by (3.1), together with (3.2) and (3.36), as well 
as the LO equivalent photon distribution of the electron 7®(4'5), 




27r 




mt 


(4.15) 
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where me is the electron mass. Equation (4.15) is obtained from (1.4), retaining only the 
leading logarithmic term and identifying the scale /i^ with s. The phase space elements 
of the two muons can be written as 


2E, 


- d^ZTi <5yi = TT d/|, d?/i, 


(4.16) 


where Ixi = |ZrJ, * = 1,2, and the azimuth integrations have been carried out in the last 
equality. Furthermore, the original four-dimensional ^-function in (4.14) can be split into 
its energy, transverse momentum and longitudinal momentum parts: 

+ k 2 -li-l 2 ) = h - % cosh 1/1 - cosh 1 / 2 ^ - It^) 

^ ^ ~ ~ sinh|/i - It^ sinh|/ 2 ^ , (4.17) 


therefore, at lowest order, the transverse momentum components of the ^-function ensure 
that the muons are produced with equal and opposite transverse momenta. We dehne 


Pt = Iti — I 


T2 


(4.18) 


and the integrations over ^ and x in (4.14) can be carried out using the two remaining 
(5-functions in (4.17). Finally, making use of (4.16), (4.18) and the dehnition 

d(j 1 


di IhTTS^ 


|M|2, 


we get [77] 


——y = ^7^(C, s) s) ^ , 


(4.19) 


(4.20) 


dyidy2dpT ' ' ' ' dt 

where the dependence of momentum fractions ^ and x on the variables yi, y 2 , pr is given 
by 


and 


{ = (e" + e»"), 


X = (e -t- e ^^) 


(4.21) 


(4.22) 


The dimuon invariant mass squared s in (4.5) can also be expressed in terms of the 
rapidities of the two muons and the transverse momentum of one of them as 


s = 2p^[l -f cosh(//i - 1 / 2 )]; 


(4.23) 
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using this relation one can derive, from (4.20), the cross section differential in yi, y 2 and 
s: 

. = i^rT—- h -^ s) xY{x, s) ^ , (4.24) 

dyidy 2 ds 2[1 + cosh(|/i - y 2 )\ dt 

with 

i = —^PTVse~^^, u = —^pTVse~'^'^. (4-25) 


At HERA (s = AEgEp) rapidities are commonly measured along the proton beam 
direction, hence one should replace yi with —yi (or, equivalently, exchange ^ with x) 
in (4.21), (4.22) and (4.25), since the e — p center-of-mass rapidities yi were deffned to 
be positive in the electron forward direction. Being rapidities additive quantities under 
successive boosts, the laboratory-frame rapidities of /i+ and p.“, pi and 772 , are related to 
yi and 7/2 by 


r]i = yi + In 



(4.26) 


where the last term in (4.26) is the rapidity relative to the boost along the axis from 
the laboratory to the center-of-mass frame, calculated according to (4.11) with velocity 
f3 = {Ep — Ee)/{Ee + Ep), Ep and E^ being the colliding proton and electron energies. In 
terms of pi and 772 , (4.21) and (4.22) are given by 


e = 


X = 


Vs 1 



(4.27) 



1 ’ 


/ e»?i + e»?2 ^ 


(4.28) 

2Ep 

4- e-^% 

1 ’ 


and can be used to estimate the dimuon production process in the laboratory-frame, to¬ 
gether with (4.20) and (4.24)-(4.26). Alternatively, we can choose as independent variables 
771 , 772 and using the relation 


5 = AeEj 


aVl 


+ e 


m 


(4.29) 

obtained from (4.27), we are able to calculate the Jacobian of this change of variables and 
finally we get [ 2 ] 

da d^El 


1 cosh( 77 i - 772 ) e *?! + e dt 


d77id772d,^ 

where the cross section for the subprocess 7 ^ 7 ^ 
laboratory-frame, 

dd 2^0^ ft u\ 47ra^ 

di In i 


(4.30) 

given in (B. 8 ) reads, in the 


cosh (771 - 772 ). 


(4.31) 
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The last equality in (4.31) follows from (4.25) with Ui —> —yi and (4.26), that is 

i = —2^pTEee^^, u = —2^pTEee^^. (4.32) 

Furthermore, from (4.23) and (4.29), remembering that pi — yj = pi — pj, one hnds 


Pt = 


‘2^E, 


e-^i + e-’?2 


(4.33) 


At the hxed-target experiment HERMES (s = 2mEe), where the axis is chosen to 
be along the electron beam, (4.6)-(4.25) still hold and (4.26) has to be replaced by 


Pi = yi + In 


2E, 


m 


(4.34) 


as now Eg/[Eg + m) is the velocity of the boost from the laboratory to the center-of-mass 
frame. Therefore, in (4.27)-(4.32) one has to make the following replacements 

m 


Ej, —>■ — , 
P o ’ 


Vi 


-Vi 


(4.35) 


with Pi now corresponding to the rapidities of the observed particles with respect to the 
electron beam direction. 

Furthermore, at HERMES one may study also p^{x,s) as well as the polarized 
given in (3.39), (3.58) and (3.63), by utilizing, from (1.10), 


Ar(e,5) = 


a l-(l-0' 


In — 


(4.36) 


2t[ ^ ml 

in the spin dependent counterpart of (4.30), while the relevant LO cross section for the 
polarized subprocess is given by (B.13), namely 

dA(j dd 

d£ di 


(4.37) 


4.1.2 Electron-Photon Production 

For the Compton process ep —> epX proceeding via the subprocess 

e(/) + Yik) e{l') + p{k'), (4.38) 

as depicted in Figure 4.2, we dehne the variables 


s = (/ + kY, i = {I — I'Y, u = {I — k'Y- 


(4.39) 
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e 


P 




Figure 4.2: Lowest-order Feynman diagram for Compton scattering in ep collisions. The 
crossed u-channel contribution is not shown. 


The kinematics of the process is quite similar to the one of the reaction ep pPprX, 

discussed above. In particular if one fixes ^ = 1 and drops the terms d^ and 
(4.6)-(4.8), (4.14), (4.16)-(4.33) are still valid, with the obvious replacements 

ki I, k2 —>■ k, h k, I 2 —>■ k', 

2/2 ^ 2 / 7 , Ve, (4-40) 


where |/e ,7 and ry,^^ are respectively the center-of-mass and laboratory rapidities of the 
produced (outgoing) electron and photon. Hence at HERA (4.30) is substituted by [2] 


da 




pne _j_ g»?7 


da 


- -—--- xY(x, sj— , 

1 + cosh(?7e — V-y) 6“'^= + c-'^'y dt 


dr]edrj.y 

where x is fixed by (4.28), that is s = AxEeEp. According to (B.34), 

dd 2TTa^ 

dt P 

with 


s u 

~ H—r 
u s 


- - = 1 + ^77, 

u 


(4.41) 


(4.42) 


(4.43) 


which can be derived from (4.29), (4.32) and (4.33). 

At the HERMES experiment s = 2xmEe, and (4.41)-(4.43) still hold, but with r 7 e ,7 —^ 
—rje^^. The equivalent of (4.41) for longitudinally polarized incoming particles is obtained 
by replacing the photon distribution plx, s) and subprocess cross section da/dt with 
their spin dependent counterparts Ap{x,s) and dAd/dt. From (B.37), 

dAd 27ra^ 


dt 


(4.44) 
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Figure 4.3: Event rates for QED Compton (ej ej) and dimuon production (77 — 
) processes at the HERA collider. The cuts applied are as described in the text. 

These expressions, as well as the ones relative to dimuon production at the HERMES 
experiment presented in the previous section, apply obviously also to the COMPASS pp 
experiment at CERN whose higher incoming lepton energies {E^ = 50 — 200 GeV) enable 
the determination of at lower values of x as compared to the corresponding 

measurements at HERMES. (Notice that for a muon beam one has obviously to replace 
me by in (4.15) and (4.36)). 

4.2 Numerical Results 

We shall present here the expected number of events for the accessible x-bins at HERA 
collider experiments and at the hxed-target HERMES experiment subject to some rep¬ 
resentative kinematical cuts which, of course, may be slightly modihed in the actual 
experiments. These cuts entail s > Smin, 7min A Vi < Vinux and Ei > Emm where Ei are 
the energies of the observed outgoing particles. 
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Figure 4.4: Expected statistical accuracy of the determination of'y^{{x), {§)) via the (a) 
QED Compton process and (b) the dimuon production process at the HERA collider. The 
numbers indicate the average scale (s) (in GeV^ units) for each x-bin. 


The relevant integration ranges at HERA are fixed via 

^mi' 




< a; < 1, 


Af^EgEp 

with s given by s = AxfEeEp while rji are constrained by 

, xEr, 

Vi + V2 = In 


which follows from (4.27) and (4.28). The relation 


E,, 


(l + e^’^O - 1 


(4.45) 


(4.46) 


T]i - r]j = In 


(4.47) 
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Figure 4.5: Event rates for the QED Compton process at HERMES using an (un)polarized 
proton target. The upper (solid and dashed) curves refer to an unpolarized proton, whereas 
the lower ones refer to a polarized proton target. The cuts applied are as described in the 
text. 


as obtained from the outgoing particle energy Ei = p^ cosh rg and its transverse momen¬ 
tum (4.33), further restricts the integration range of rji^ as dictated by Ei > For 

the QED Compton scattering process (4.41), .^ = 1, ? 7 i = pe, P 2 = V'y (4.45)-(4.47). At 
HERMES Ep ^ m/2 and pi —> —pi in the above expressions with pi the outgoing particle 
rapidity with respect to the ingoing lepton direction. 

In the following we shall consider Emin = 4 GeV. For the QED Compton scattering 
process we further employ Smin = 1 GeV^ so as to guarantee the applicability of perturba¬ 
tive QCD, i.e., the relevance of the utilized s), see (4.3) with p? = s. For the dimuon 
production process we shall impose Smin = m^[\['(2S')] = (3.7 GeV)^ so as to evade the 
dimuon background induced by charmonium decays at HERMES (higher charmonium 
states have negligible branching ratios into dimuons); for HERA we impose in addition 
-Smax = m^[T(lS')] = (9.4 GeV)^ in order to avoid the dimuon events induced by bot- 
tomium decays. Finally, at HERA we consider ?7min = —3.8, ?7max = 3.8 and at HERMES 
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Figure 4.6: Expected statistical accuracy of the determination of (a) 7 ^((a;), (s)) and (b) 
A'yP{{x), {§)) via the QED Compton process at HERMES using an (un)polarized proton 
target. The numbers indicate the average scale (s) (in Ge]R units) for each bin. 


Vmin = 2.3, 7max = 3.9. The integrated luminosities considered are /^hera = 100 pb ^ and 
-^HERMES = 1 fb ^ 

In Figure 4.3 the histograms depict the expected number of dimuon and QED Compton 
events at HERA found by integrating (4.30) and (4.41) applying the aforementioned 
cuts and constraints. The important inelastic contribution due to being calculated 
according to (4.3) using the minimal boundary condition, is shown separately by the 
dashed curves. 

To illustrate the experimental extraction of 7 ^(x, s) we translate the information in 
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Figure 4.7: As in Figure 4-5 but for a neutron target. The negative signs at some lower-x 
bins indicate that the polarized total cross section and/or inelastic contribution is negative. 


Figure 4.3 into a statement on the accuracy of a possible measurement by evaluating 
(s)) at the averages {x), {§) determined from the event sample in Figure 4.3. 
Assuming that in each bin the error is only statistical, i.e. S'j^ = ± 7 ^/ \/Nun, the results 
for x'y^/a are shown in Figure 4.4. It should be noticed that the statistical accuracy 
shown will increase if 7 ^ei(^>/^o) 7 ^ 0 in contrast to our vanishing boundary condition 
used in all our present calculations. Our results for the QED Compton process in Figures 
4.3 and 4.4 are, apart from our somewhat different cut requirements, similar to the ones 
presented in [27]. 

Apart from testing 'y^{x, s) at larger values of x, the fixed-target HERMES experiment 
can measure the polarized A'yP{x, s) as well. In Eigure 4.5 we show the expected number 
of QED Compton events for an (un)polarized proton target. 

The accuracy of a possible measurement of ■y^/x), {§)) and A'y^/x), {§)) is illustrated 
in Eigure 4.6 where the averages (x), (s) are determined from the event sample in Eigure 
4.5 by assuming that the error is only statistical also for the polarized photon distribution. 
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Figure 4.8: As in Figure 4-6 but for a neutron target. 


i.e. 5(A7^’) = ±{VN^/\ANun\)AY. 

The analogous expectations for an (un)polarized neutron target are shown in Figures 
4.7 and 4.8. It should be pointed out that, according to Figures 4.6(b) and 4.8(b), 
HERMES measurements will be sufficiently accurate to delineate even the polarized 
distributions in the medium- to small-x region, in particular the theoretically more spe¬ 
culative inelastic contributions. 

For completeness, in Figures 4.9 and 4.10 we also show the results for dimuon produc¬ 
tion at HERMES for (un)polarized proton and neutron targets despite the fact that the 
statistics will be far inferior to the Compton process. 

The dimuon production can obviously proceed also via the genuine Drell- Yan subpro- 
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Figure 4.9: As in Figure 4-5 but for dimuon production at HERMES using (un)polarized 
proton and neutron targets. The lower solid and dashed curves refer to a polarized nucleon 
target and the negative signs indicate that the polarized cross sections are negative. 

cess qq —> pAp,~ where one of the (anti)quarks resides in the resolved component of the 
photon emitted by the electron. However, as already noted in [79], this contribution is 
negligible as compared to the one due to the Bethe- Heitler subprocess 77 —. The 
unpolarized dimuon production rates at HERA where also studied in [79,80] utilizing, 
however, different prescriptions for the photon content of the nucleon. 

Exact expressions for the Bethe-Heitler contribution to the longitudinally polarized 
-yN — p~^p~X process are presented in [81] but no estimates for the expected production 
rates at, say, HERMES or COMPASS are given. 

4.3 Summary 

The analysis of the production rates of lepton-photon and muon pairs at the colliding 
beam experiments at HERA and the fixed-target HERMES facility, as evaluated in the 
leading order equivalent photon approximation, demonstrates the feasibility of determi¬ 
ning the polarized and unpolarized equivalent photon distributions of the nucleon in the 
available kinematical regions. The above mentioned production rates can obviously be de¬ 
termined in a more accurate calculation along the lines of [34], involving the polarized and 
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Figure 4.10: As in Figure 4-6 but for dimuon production at HERMES for (un)polarized 
proton and neutron targets. The statistical accuracy for the inelastic contributions is 
similar to those shown for the total result, except for the almost vanishing A 7 ;”g[. 

unpolarized structure functions gi ^2 and Fi, 2 , respectively, of the nucleon. The expected 
production rates are similar to those obtained in our equivalent photon approximation, as 
discussed in detail in the next chapters. It thus turns out that lepton-photon and muon 
pair production at HERA and HERMES may provide an additional and independent 
source of information concerning these structure functions. 

























Chapter 5 

The Unpolarized QED Compton 
Scattering Process 


The QED Compton scattering in high energy electron-proton collisions ep eyX is one of 
the most important processes for an understanding of the photon content of the proton. 
In addition, it can also shed some light on the proton structure functions 
[31,33,34] in the low-Q^ region, where they are presently poorly known [26]. 

The QED Compton scattering has been recently analyzed in [31], where the above 
mentioned alternative descriptions were confronted with the experimental data and it 
was found that the description in terms of Fi^ 2 , he. for X ^ p, is superior to the one in 
terms of the inelastic photon distribution Henceforth we shall refer to the 

description in terms of Fi ^2 as exact to distinguish it from the approximations involved in 
the EPA. 

It should be noted, however, that the analysis in [31, 34] utilized the Callan-Gross 
relation [46] Fl^xb^Q"^) = F2{xb,Q‘^)—‘2^xbFi{xb,Q‘^) = 0. This relation is contaminated 
by higher order (NLO) QCD corrections as well as by higher twist contributions relevant 
in the low-Q^ region which may invalidate the assumptions underlying the exact analysis. 
Furthermore, the analysis in [31,34] was carried out within the framework of the helicity 
amplitude formalism [34]. The implementation of experimental cuts within this formalism 
is nontrivial and affords therefore an iterative numerical approximation procedure [34,82] 
whose first step corresponds to —k"^ = = 0, where k is the momentum of the virtual 

photon. 


It is this second issue that we intend to study here. We shall replace the noncovariant 
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helicity amplitude analysis of [34] by a standard covariant tensor analysis whose main 
advantage, besides compactness and transparency, is the possibility to implement the 
experimental cnts directly and thus avoid the necessity of employing an iterative approx¬ 
imation of limited accnracy. The hrst issne concerning the Fl contribntions affords some 
estimates of this poorly known strnctnre fnnction and we refrain from its study here. 

In Section 5.1 it is specihed what QED Compton scattering is and how it can be 
selected from the reaction ep —eyX. In Sections 5.2 and 5.3, we calculate its exact 
cross section for the elastic and inelastic channels. The nnmerical resnlts are discussed in 
Section 5.4. The snmmary is given in Section 5.5 and the kinematics in Appendix C. 


5.1 Radiative Corrections to Electron-Proton Colli¬ 
sions 

The lowest order Feynman diagrams describing the process ep —> eyX, with a real photon 
emission from the electron side, are shown in Figure 5.1. The corresponding amplitudes 
contain the denominators -|- /')^ and k‘^{l — k'Y, therefore the main contribntions to 
the cross section come from conhgurations where one or both denominators tend to zero. 
These conhgnrations have different experimental signatnres and they are described in the 
following [31,34]. 

• The dominant contribntion stems from the so-called bremsstrahlung process, which 
corresponds to a conhgnration where k"^, (/' -|- k'Y, (/ — k')"^ stay all close to zero. It 
involves qnite large connting rates bnt in order to be measnred, reqnires a specihc 
small-angle detector, because the outgoing electron and photon have small polar 
angles and escape throngh the beam pipe. 

• Either (/' -|- k')"^ ~ 0 or (/ — k'Y — 0, bnt the momentnm sqnared of the exchanged 
virtnal photon k"^ is hnite: the ontgoing photon is emitted either along the hnal 
or the incoming electron line and this conhgnration corresponds to the so called 
radiative corrections to electron-proton scattering. 

In the hrst case, (/' -|- k'Y — 0, the cross section is dominated by the contribntion 
given by the hrst diagram, and this kind of events is called Final State Radiation 
(FSR). It is nsnally very difficult to distinguish this process experimentally from 
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1 



Figure 5.1: Feynman diagrams considered for ep —> e'yX, with a real final state photon 
(k'^ = 0). 

non-radiative deep inelastic scattering ep eX, since the outgoing electron and 
photon are almost collinear. 

In the second case, {I — k'fi ~ 0, the main contribution to the cross section is given 
by the second Feynman diagram in Figure 5.1 and such events are classified as Initial 
State Radiation (ISR). In the detector one observes only the outgoing electron, the 
final photon being emitted along the incident electron line. 

• The virtuality of the exchanged photon is small, ~ 0, but both (/' + k'fi and 
(/ — k'Y are finite: the produced hadronic system goes straightforwardly along the 
incident proton line, the outgoing electron and photon are detected under large polar 
angles and almost back-to-back in azimuth, so that their total transverse momentum 
is close to zero. This configuration is referred to as QED Compton scattering, since it 
involves the scattering of a quasi-real photon on an electron. This process will thus 
be selected by performing a cut on the total transverse momentum of the outgoing 
electron and photon or on the acoplanarity (5.47) of the electron-photon system. 

As pointed out in [26], the corresponding cross section is large, despite the fact that 
it contains an additional factor a = 1/137 compared to the tree-level cross section 
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for ep eX. This follows since the emission of a large transverse momentum 
photon can lead to a reduction of the true momentum squared transferred to the 
proton: (/ — I'Y in ep eX is shifted to {I — I' — k'Y in ep —eqX, which can 
become of the order of the proton mass squared or even smaller, as will be discussed 
below, see (5.45). The reduction in and the corresponding increase of the cross 
section in (5.40), compensates for the smallness of the additional factor a. Hence 
QED Compton scattering can provide a tool to investigate the small-Q^ behaviour of 
the proton structure functions. The first measurements of F 2 using QED Compton 
scattering have been published by the HI collaboration at HERA [33]. 


5.2 Elastic QED Compton Scattering 

We consider elastic QED Compton scattering: 

e{l)+p{P)^e{l')+^{k')+p{P'), (5.1) 

where the four-momenta of the particles are given in the brackets. We introduce the 
invariants 


s = (P + /)2, t = k‘^, (5.2) 

where k = P — P' is the four-momentum of the exchanged virtual photon. Moreover, we 
will make use of the Mandelstam variables (4.39) relative to the subprocess e{l)'y{k) —> 
e{l')'y{k'). The photon in the final state is real, = 0. We neglect the electron mass 
everywhere except when it is necessary to avoid divergences in the formulae and take the 
proton to be massive, = P'^ = m?. The relevant Feynman diagrams for this process 
are shown in Figure 5.1, with X being a proton and Px = P'■ The squared matrix element 
can be written as 

VPP = k'), (5.3) 

where ,k') is the leptonic tensor (B.31), given also in [75,83] and H^f{P,P') 

is the hadronic tensor defined in the first line of (2.18), with = p, in terms of the 
electromagnetic current J^. If we use the notation 

dPS„(p;pu = (2x)‘‘5(p - J]pi) n 


(5.4) 
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for the Lorentz invariant iV-particle phase-space element, the total cross section will be 

^ f dPS2+i{l + P; k', P')\ Mel I' . (5.5) 

Equation (5.5) can be rewritten following the technique of [24], which we slightly modify 
to implement the experimental cuts and constraints; in particular all the integrations will 
be performed numerically. Rearranging the (2 -|- l)-particle phase space into a sequence 
of two 2-particle ones, (5.5) becomes: 

<T,i(s) = j iPS^(l + P-,l‘ + k', P')j^Hf(P, P')X„4l, k) . (5.6) 

The tensor Xap contains all the informations about the leptonic part of the process and 
is dehned as 

Xag{l, k) = J dPS2{l + k- k')T^f3{l; k') (5.7) 

and Ta/s can be written as [3] 

Tafsih I', k') = / - Qap {s^ + + 2tt) -|- 2s I alp + T.U 

SU I ^ 

+ (^ + + h^a) ~ ~ {^Q-k'p -f Igk'^ 

+ + (5-8) 

It can be shown that 

dPS^(l + k; I', k') = /‘y , (5.9) 

lo TT^yS — t) 

with (p* denoting the azimuthal angle of the outgoing e — 7 system in the e — 7 center-of- 
mass frame. For unpolarized scattering, Xap is symmetric in the indices a, (3 and can be 
expressed in terms of two Lorentz scalars, Xi and X 2 : 

Ws(i,*:) = 2||3y(S,«)+X2(S,«)]([j^i-J;]) 

+ [7(9,f) + X2{s,t)]{t(Ja:fi - ko.k/3)\. (5.10) 

Mi,t) ^ jCt^l‘-rx„i,{l,k), 

(s - t)‘ 


with 


(5.11) 
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Using the leptonic tensor (5.8) and also the relations 


we obtain 


-tt 


47r2(s —1)3 2Ti‘^{s — tY 


= e 




_ 4 - 2ts + 25^ + 2st + , A 

— e -r—-—-^-=A2(^S,t,tJ. 


167r^(.s — t) At[‘^s{s — t){t — s — i) 

The invariants Xi{s,t,i), with i = 1,2, are related to Xi{s,t) by 


Xi{s,t) = 27r / diXi{s,t,i). 


The integration limits of i are: 


^max 0) ^min <S T t T 


S — t 




(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 

(5.17) 


where rrie is the mass of the electron. We point out that the kinematical cuts employed 
by us prevent the electron propagators to become too small and thus the divergences 
are avoided, so we can safely neglect the electron mass in the numerial calculation. The 
hadronic tensor in the case of elastic scattering can be expressed in terms of the common 
proton form factors as in (3.7), namely 


Hf{P, P') = [i7i(t)(2P - kY{2P - kf + Glj{t){tg^^ - k^k% (5.18) 


where Pi(t), already introduced in (3.3), is given by 

GW) - (tUrWGlW) 


H,(t) = 


1 — t/Aw? 


Using 


hnally we get [3] 


o-eds) = 


dPS2{l + P]l' + k\P') = 


dt 


8n{s — m?) 


(5.19) 


(5.20) 


a 


87r(s — w?Y 


pin 

ds 


dt 


max ^ 4 - / ‘'max 


r*27r 


dt / d(fi^ 


‘'min ^ ‘'min 


s — m 
s — t 
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s — t 


o — 777^ \ ^ /s 2777'^ ^ ^ 

X I —-1 ] [3Xi(s, t, t) + X2{s, t, t)] H-—[Xi(s, t, t) + X2{s, t, t)] 


t 


+ t, i) 


H^{t)+X2{sC,i)Gl,{t)\, 


(5.21) 


where Smin denotes the minimum of s and tmin.max are given by 


^min, max 


= 2m^ — 


2s L 


(s + wS){s — s + mS) ± (s — mS) a/(s — s + — Asm? 


(5.22) 


It is to be noted that in (5.21) we have shown the integration over (p* explicitly, because 
of the cuts that we shall impose on the integration variables for the numerical calculation 
of the cross section. The cuts are discussed in Section 5.4. 

The EPA consists of considering the exchanged photon as real, so it is particularly 
good for the elastic process in which the virtuality of the photon |t| is constrained to be 
small (< 1 GeV^) by the form factors. It is possible to get the approximated cross section 
^EPA exact one in a straightforward way, following again [24]. If the invariant 

mass of the system e —7 is large compared to the proton mass, Smin 3> mS, one can neglect 
\t\ versus s, w? versus s, then 


Xi(s,t,t) ^ Xi(s,0,t) = 0, 


(5.23) 


and 


-N 2s d(j(s, i) 
X 2 (s, t, t) ^ X 2 (s, 0 , t) = - 


(5.24) 


TT dt 

where the differential cross section for the real photoproduction process ey —> ey is given 
in (B.34). We get: 


^el(.) ^ = 


yi-m/y3)2 


dx 




d(j(a;s, i) 
dt 


(5.25) 


where x = s/s and is the elastic contribution to the equivalent photon distribution 

of the proton (3.2): 


^ ‘'min 


1 fl 


~ I — ^~r 

X \x / t 


m 


21 


Hi(t) + G^(t) >, 


with 


2 2 
m X 

1 — X 


(5.26) 


(5.27) 
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To clarify the physical meaning of x, let us introduce the variable 

l-k , , 

I, = — , (5.28) 

It is possible to show that represents the fraction of the longitudinal momentum of the 
proton carried by the virtual photon, so that one can write 


k = x^P + k, 


with k ■ I = 0. Using (5.2) one gets 


x^ 


s — t 

9 ’ 

s — 


(5.29) 


(5.30) 


which reduces to x in the EPA limit, see also (3.10) and (3.11). One can also dehne the 
leptonic variable xp 


Xl 


Q1 

2P-{1- /') ’ 


(5.31) 


where Qf = —t. When t ~ 0, it turns out that also xi ~ x. 


5.3 Inelastic QED Compton Scattering 

To calculate the inelastic QED Compton scattering cross section, we extend the approach 
discussed in the previous section. In this case, an electron and a photon are produced in 
the hnal state with a general hadronic system X. In other words, we consider the process 

e(/) + p{P) ^ e(r) + 7(A:') + X{Px), (5.32) 

where Px = Xlx momenta of the produced hadrons. Let the 

invariant mass of the produced hadronic state X to be lU; (5.2) still holds with = —t. 
The cross section for inelastic scattering will be 

O'inel(’S) = 7 -T-yr / dPS' 2 +v(^ +-P; -Pvi, •• •,-Pvjv) I ^inel |^, ( 5 . 33 ) 

2[S — m^) J 

where 

M|„el P = Cv)r„;,(i, k; I', k') 


(5.34) 
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is the squared matrix element and the tensor Px) has already been introduced in 

(2.42). If we rearrange the (2 + iV)-particle space phase into a sequence of a 2-particle 
and a iV-particle one, we get 

(5.35) 

where is given by (5.10) and is the hadronic tensor for inelastic scattering 

dehned in (2.41), that is, using the notation (5.4), 


= I dPSxiP - k- Px,,Px^) //“fi • 


(5.36) 


The hadronic tensor is parametrized in terms of Fi and F 2 , the usual structure functions 
of the proton, as in (3.8), i.e. 

-(QY^ + k‘^k^)Fi(xB,QA 

Q V 


+ {2xbP‘^ - k^) {2 xbP‘^ - k^)^ ( 5 . 37 ) 

2xb 


where is the Bjorken variable (2.37), 


Xb = 


2P ■ i-k) Q^ + W^-m'^ ' 


Using 


as before, we get [3] 


dPS 2 {l + P;k + k\Px) = 


8n{s — nP) 


(5.38) 


(5.39) 




A7t{s - W?Y Jy^2 


M^max rip-Wf /‘Qmax ^Umax /‘2vr 

L f i 

'^min >^Smin Vmin ^ tmin U 


\\fc^s-nP{^ s-nP\ . 2 \( 2 is-nP) 1 

""uyjwi « + ““ {Q^(» + QV~P 

+ C-^Jjl3X,(S,QM) + W(I,Q^^)] + 

+ l■n(s. ot t) + ( 3 ^ £)] - .Y.(S, i) 


X F2(xb,Q^)y - 


(5.40) 
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Here Xj(s, with i = 1,2, are given by (5.14)-(5.15) with t replaced by —Q^. The 

limits of the integration over are: 


Vmi 


mm,max 


2s 


(s + m?){s — s + W'^) 

T (s - m^) y/{s-s + - 4 swA , (5.41) 


while the extrema of t are the same as (5.17). The limits hh^nmax given by: 


^min (m + TTljr) , ^max (v^ '\/'Smin ) , 


(5.42) 


where is the mass of the pion. 

In the EPA, we neglect compared to s and compared to s as before. Using 
(5.23) and (5.24), we can write 


^inel('S) 


_EPA 

^inel 




dx 


dt 7rnel(a^> XS. 


dd(a;s, i) 
dt 


(5.43) 


where again x = s/s and if^^ex^x^xs) is the inelastic contribution to the equivalent photon 
distribution of the proton [84]: 


7fnel(a^> XS) 



( i + ii-yf 
V 


2m^x‘^\ 

(5.44) 


The limits of the integration can be approximated as 


Vmm 


2 2 
X m 

1 — a; ’ 


(5.45) 


and we choose the scale Q^ax to be s. Our expression of 7 ?g[(a;,a;s) differs from [26] by a 
(negligible) term proportional to m?. Following [31, 34] we shall use the LO Callan-Gross 
relation 


Fl{xb,Q^) = F2{xb,Q^) -‘2xbFi{xb,Q‘^) = 0 (5.46) 


in our numerical calculations. 
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5.4 Numerical Results 

In this section, we present an estimate of the cross section, calculated both exactly and 
using the equivalent photon approximation of the proton. We have used the same kine- 
matical cuts as used in [31] for the HERA collider, which are slightly different from the 
ones in [34], They are imposed on the following laboratory frame variables: energy of the 
final electron E', energy of the final photon polar angles of the outgoing electron and 
photon, 6e and 6^ respectively, and acoplanarity angle 0, which is defined as 

0 = I TT - 10^ - 0e| I, (5.47) 

where (j)^ and 0e are the azimuthal angles of the outgoing photon and electron respectively 
(0 < (j)^, 0e < 2 7r). The cuts are given by: 

E', E; > 4 GeV, E'^ + E'^> 20 GeV, (5.48) 

0.06 < 0e,^7 < TT-O.Oe, (5.49) 

0 < 0 < ^ . (5.50) 

The energies of the incoming particles are: E^ = 27.5 GeV (electron) and Ep = 820 GeV 
(proton). In our conventions, we fix the laboratory frame such that 0e = 0, so the 
acoplanarity will be 0 = \'k — These cuts reflect experimental acceptance constraints 
as well as the reduction of the background events due to emitted photons with {V+k'Y ~ 0 
and/or (/ — k'Y ~ 0, which are unrelated to the QED Gompton scattering process (for 
which —k"^ = 0 but with both (/' + k'Y and (/ — k'Y finite), i.e. photons emitted 

parallel to the ingoing (outgoing) electron [26,34], as discussed in Section 5.1. 

In contrast to [34], we find that the cuts (5.48)-(5.50) are not sufficient to suppress the 
background due to photon emission at the hadron vertex, discussed in the next chapter. 
Such a background has been subtracted from the measurements in [31] and from the 
numerical estimates in [31,32], which will be shown for comparison with our results in 
Figure 5.2 and in the Tables 5.1 and 5.2. 

We numerically integrate the elastic and inelastic cross sections given by (5.21) and 
(5.40). To implement the cuts in (5.48)-(5.50), we express E', E/, cos6*e, cos6*^ and cos0 
in terms of our integration variables s, f, f, (p* (and in the inelastic channel), as 
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explained in Appendix C. More explicitly, we use (C.27)-(C.31), (C.12), (C.13), together 
with (C.16)-(C.18) for the elastic channel and (C.36)-(C.38) for the inelastic one. The 
cuts imposed on the laboratory frame variables restrict the range of our integrations nu¬ 
merically. In this way, we are able to remove the contributions from outside the considered 
kinematical region. 

In the calculation of the elastic cross section, the electric and magnetic form factors 
are empirically parametrized as dipoles 

GEit) = 7-—^^ , GmW = 2.79 GEit). (5.51) 

Following [31], in the evaluation of (5.40) we have used the ALLM97 parametrization of 
the proton structure function F 2 {xb,Q‘^) [70], which provides a purely phenomenological, 
Regge model inspired, description of F2{x b, Q^), including its vanishing in the = 0 limit 
as well as its scaling behaviour at large Q^. The ALLM97 parametrization is supposed 
to hold over the entire range of xb and studied so far, namely 3 x 10“® < xb < 0.85 
and f) < < 5000 GeV^, above the quasi-elastic region (hF^ > 3 GeV^) dominated 

by resonances. We do not consider the resonance contribution separately but, using the 
so-called local duality [85], we extend the ALLM97 parametrization from the continuous 
(IV^ > 3 GeV^) down to the resonance domain ((m^ -|- m)^ < < 3 GeV^): in this 

way it is possible to agree with the experimental data averaged over each resonance, as 
pointed out in [31]. 

The elastic contribution to the ERA was calculated using (5.25) subject to the addi¬ 
tional kinematical restrictions given by (5.48)-(5.49). For the inelastic channel we have 
used (5.43) together with (5.44), the cuts being the same as in the elastic case. We have 
taken Fe = 0 and used the ALLM97 parametrization of F 2 , in order to compare consis¬ 
tently with the exact cross section. We point out that in (3.2), as well as in [1,2,27], 
F 2 {xb, Q^) in 7((iei(^) expressed in terms of parton distributions for which the LO 

GRV98 parametrization [68] was used, together with Q^in = 0-26 GeV^ so as to guaran¬ 
tee the applicability of perturbative QGD [25]. The new ll^^ex^x^xs) gives slightly higher 
results than the ones obtained with the photon distribution in (3.2). 

The Gompton process turns out to be dominated by the elastic channel, in fact after 
Monte Garlo integration, we hud that dei = 1.7346 nb, while dinei = 1.1191 nb. The 
approximated calculation gives the results: = 1.7296 nb and = 1.5183 nb. This 

means that in the kinematical region under consideration, the total (elastic -|- inelastic) 
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Figure 5.2: Cross section for Compton process at HERA-Hl. The cuts applied are given 
in (5.48)-(5.50). 


Exact: 


total 

elastic 


EPA: - total 


cross section calculated using the EPA agrees with the exact one within 14% and that 
the approximation turns out to be particularly good in describing the elastic process, 
for which the agreement is within 0.3%. This is not surprising since in the EPA one 
assumes = 0, which is not true especially in the inelastic channel and the inelastic 
cross section receives substantial contribution from the non-zero region. In terms of 
the kinematical cuts, the EPA corresponds to the situation when the outgoing electron 
and the hnal photon are observed under large polar angles and almost opposite to each 
other in azimuth, so that the acoplanarity is approximately zero. For elastic scattering 
there is a sharp peak of the exact cross section for 0 = 0, contributions from non-zero 0 
are very small in this case. But the inelastic cross section receives contribution even from 
non-zero 0, so that in this case the discrepancy from the approximated result is higher. 
The discrepancy of the total cross section with the approximate one is thus entirely due 
to the inelastic part. 

In Figure 5.2 we have compared the total cross sections (exact and EPA) in different 
Xi bins, in the region 1.78 x 10“^ < xi < 1.78 x 10“^. Figure 5.3 shows that the agreement 
improves slightly for bins in the variable x^. Since x^ ~ xi for — 0, the elastic process 
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Figure 5.3: Cross section for Compton process at HERA-Hl. The bins are in x^. The 
cuts applied are given in (5.48)- (5.50). 


is not sensitive to this change of variables. We point out again that in the EPA limit 

(Q^ = C) xi = = X. 

In Figure 5.4 we show the exact and the EPA cross section in xi and bins together 
with the experimental results and the estimates of the Compton event generator, already 
presented in [31]. Except for two bins, our exact result agrees with the experiment within 
the error bars. The slight difference of our exact result and the one of [31] may be due 
to the fact that in [31] the cross section is calculated using a Monte Carlo generator in a 
step by step iteration [34, 82] which starts by assuming Q'^ = 0, while we did not use any 
approximation. Our exact result is closer to the EPA in most of the kinematical bins as 
compared to [31]. The total cross section in the EPA lies above the exact one in most of 
the bins. 

For completeness, we have shown the numerical values of the exact and EPA double 
differential cross sections, both for the elastic (Table 5.1) and inelastic (Table 5.2) contri¬ 
butions. The kinematical bins are the same as in [31]. The exact results when the bins 
are in x^ instead of xi are also shown. 
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Qi^ (GeV^) (c) Qj^ (GeV^) (‘I) 


Figure 5.4: Double differential cross section for QED Compton scattering at HERA-Hl. 
The data are taken from [31]. The kinematical bins correspond to Table 5.1. The contin¬ 
uous line corresponds to our exact calculation, the dotted line to the calculation in [31] 
and the dashed line to the ERA. 
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X bin 

Q[ bin 

^el 

<'7” 

<1 


1.78 X 10-5 - 5.62 X 10-5 

1.5 -2.5 

2.428 X 102 

2.342 X 10^ 

2.446 X 10^ 

2.461 X 10^ 

1.78 X 10-5 - 5.62 X 10-5 

2.5 -3.5 

5.099 X 10^ 

4.71 X lOi 

5.201 X lOi 

5.051 X IQi 

5.62 X 10-5 - 1.78 X 10-^ 

1.5 -5.0 

5.279 X 102 

5.319 X 10^ 

5.259 X 10^ 

5.247 X 10^ 

5.62 X 10-5 - 1.78 X 10-^ 

5.0-8.5 

2.396 X 102 

2.327 X 10^ 

2.404 X 10^ 

2.395 X 10^ 

5.62 X 10-5 - 1.78 X 10-^ 

8.5 - 12.0 

8.496 X 10^ 

8.32 X lOi 

8.559 X lOi 

8.571 X IQi 

1.78 X 10-^ - 5.62 X 10-^ 

3.0 - 14.67 

2.080 X 102 

2.036 X 10^ 

2.056 X 10^ 

2.061 X 10^ 

1.78 X 10-^ - 5.62 X 10-^ 

14.67 - 26.33 

1.373 X 102 

1.388 X 10^ 

1.373 X 102 

1.372 X 10^ 

1.78 X 10-^ - 5.62 X 10-^ 

26.33 - 38.0 

3.712 X 10^ 

3.86 X lOi 

3.720 X lOi 

3.695 X IQi 

5.62 X 10-^ - 1.78 X 10-3 

10.0 - 48.33 

5.947 X 10^ 

5.71 X lOi 

5.918 X lOi 

5.921 X IQi 

5.62 X 10-^ - 1.78 X 10-3 

48.33 - 86.67 

3.714 X 103 

3.85 X lOi 

3.715 X lOi 

3.704 X IQi 

5.62 X 10-^ - 1.78 X 10-3 

86.67 - 125.0 

1.056 X 10^ 

1.028 X IQi 

1.057 X lOi 

1.054 X IQi 

1.78 X 10-3 - 5.62 X 10-3 

22 - 168 

1.913 X 10^ 

1.877 X IQi 

1.909 X lOi 

1.909 X IQi 

1.78 X 10-3 - 5.62 X 10-3 

168 - 314 

1.239 X 10^ 

1.229 X IQi 

1.239 X lOi 

1.238 X IQi 

1.78 X 10-3 - 5.62 X 10-3 

314 - 460 

5.917 

6.02 

5.915 

5.914 

5.62 X 10-3 - 1.78 X 10-2 

0 - 500 

4.811 

5.76 

4.890 

4.890 

5.62 X 10-3 - 1.78 X 10-2 

500 - 1000 

9.271 

9.22 

9.264 

9.271 

5.62 X 10-3 - 1.78 X 10-2 

1000 - 1500 

2.572 

2.65 

2.571 

2.573 

1.78 X 10-2 - 5.62 X 10-2 

0 - 1500 

8.238 X 10-1 

6.8 X 10-1 

9.085 X 10-1 

9.086 X 10-1 

1.78 X 10-2 - 5.62 X 10-2 

1500 - 3000 

2.431 

2.69 

2.430 

2.434 

1.78 X 10-2 - 5.62 X 10-2 

3000 - 4500 

6.336 X 10-1 

7.7 X 10-1 

6.328 X 10-1 

6.345 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

10 - 6005 

3.120 X 10-1 

4.27 X 10-1 

3.120 X 10-1 

3.117 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

6005 - 12000 

2.437 X 10-1 

2.13 X 10-1 

2.438 X 10-1 

2.436 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

12000 - 17995 

0.000 

0.000 

0.000 

2.461 X 10-2 


Table 5.1: Double differential (elastic) QED Compton scattering cross section, dei is the 
exact result in (5.21), corresponds to the results in [31]. The x-bins refer to xi in 
(5.31) except for a*) where they refer to x^ in (5.30). is given in (5.25) where 

X = x^. Q[ is expressed in GeV^ and the cross sections are in pb. 
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X bin 

(52 bin 

^inel 

Len 

^inel 

^inel 

EPA 

^inel 

1.78 X 10-5 - 5.62 X 10-5 

1.5 -2.5 

7.802 X 10^ 

6.170 X lOi 

7.367 X lOi 

1.483 X 102 

1.78 X 10-5 - 5.62 X 10-5 

2.5 -3.5 

2.799 X 10^ 

2.050 X lOi 

4.029 X lOi 

3.255 X lOi 

5.62 X 10-5 - 1.78 X 10-^ 

1.5 -5.0 

2.299 X 102 

2.296 X 102 

2.298 X 102 

3.967 X 102 

5.62 X 10-5 - 1.78 X 10-^ 

5.0-8.5 

1.108 X 102 

1.062 X 102 

1.450 X 102 

2.016 X 102 

5.62 X 10-5 - 1.78 X 10-^ 

8.5 - 12.0 

3.340 X 10^ 

2.890 X lOi 

6.048 X lOi 

7.751 X lOi 

1.78 X 10-^ - 5.62 X 10-^ 

3.0 - 14.67 

2.029 X 102 

2.287 X 102 

1.228 X 102 

2.104 X 102 

1.78 X 10-^ - 5.62 X 10-^ 

14.67 - 26.33 

9.644 X 10^ 

9.230 X lOi 

1.164 X 102 

1.476 X 102 

1.78 X 10-^ - 5.62 X 10-^ 

26.33 - 38.0 

2.742 X 10^ 

2.570 X lOi 

4.431 X lOi 

4.298 X lOi 

5.62 X 10-^ - 1.78 X 10-3 

10.0 - 48.33 

1.011 X 102 

1.064 X 102 

5.077 X lOi 

7.555 X lOi 

5.62 X 10-^ - 1.78 X 10-3 

48.33 - 86.67 

4.485 X 10^ 

4.590 X lOi 

5.304 X lOi 

4.897 X lOi 

5.62 X 10-^ - 1.78 X 10-3 

86.67 - 125.0 

1.228 X 10^ 

1.132 X lOi 

1.887 X lOi 

1.462 X lOi 

1.78 X 10-3 - 5.62 X 10-3 

22 - 168 

4.320 X 10^ 

4.917 X lOi 

2.225 X lOi 

2.791 X lOi 

1.78 X 10-3 - 5.62 X 10-3 

168 - 314 

1.831 X 10^ 

1.735 X lOi 

2.117 X lOi 

1.849 X lOi 

1.78 X 10-3 - 5.62 X 10-3 

314 - 460 

6.314 

5.760 

8.303 

9.046 

5.62 X 10-3 - 1.78 X 10-2 

0 - 500 

1.277 X 10^ 

1.432 X lOi 

6.438 

7.627 

5.62 X 10-3 - 1.78 X 10-2 

500 - 1000 

1.086 X 10^ 

9.890 

1.152 X lOi 

1.450 X lOi 

5.62 X 10-3 - 1.78 X 10-2 

1000 - 1500 

2.734 

2.600 

3.201 

4.067 

1.78 X 10-2 - 5.62 X 10-2 

0 - 1500 

2.787 

2.500 

1.321 

1.439 

1.78 X 10-2 - 5.62 X 10-2 

1500 - 3000 

3.118 

2.150 

3.149 

3.855 

1.78 X 10-2 - 5.62 X 10-2 

3000 - 4500 

7.718 X 10-1 

6.600 X 10-1 

8.421 X 10-1 

1.004 

5.62 X 10-2 - 1.78 X 10-^ 

10 - 6005 

7.203 X 10-1 

1.460 X 10-1 

4.677 X 10-1 

4.924 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

6005 - 12000 

3.739 X 10-1 

2.110 X 10-1 

3.830 X 10-1 

3.845 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

12000 - 17995 

3.738 X 10-2 

4.300 X 10-2 

4.182 X 10-2 

3.849 X 10-2 


Table 5.2: Double differential (inelastic) QED Compton scattering cross section, ainei is 
the exact result in (S.fO), <7^®° corresponds to the results in [31], The x-bins are as in 
Table 5.1, i.e. refer to xi in (5.31) except for where they refer to x.y in (5.30). (Tjnei^ 
is given in (5.43) where x = x.y. Qf is expressed in GeV^ and the cross sections are in 
pb. 
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The EPA elastic cross section agrees within 1% with the exact one for all the xi bins. 
The agreement becomes slightly better if we consider bins. For the inelastic channel, 
the discrepancies from the EPA results are considerably higher. Our ’exact’ results he 
closer to the EPA compared to [31] in almost all the bins. The result in x^ bins shows 
better agreement with the EPA compared to the xi bins, especially for higher x^. The 
discrepancy with the EPA is about 20 — 30% in most of the bins, higher in some cases. 

5.5 Summary 

In this chapter, we have estimated both the elastic and the inelastic QED Compton 
scattering cross sections for unpolarized incoming electron and proton. Our approach 
for the calculation of the total cross section is manifestly covariant and we have used 
the same cuts as in [31]. The numerical estimates of the exact cross section for different 
kinematical bins are presented and compared with the EPA and the experimental results. 
The exact cross section in the elastic channel agrees within 1% with the approximate one. 
The discrepancy is thus due to the inelastic channel and, in the next chapter, new cuts 
will be suggested in order to reduce it. A comparison with the Monte Carlo results of [31] 
is also shown. For both elastic and inelastic cross sections, our exact results are closer to 
the EPA as compared to [31]. The agreement is even better if the bins are in x^ instead 
of Xi. 




Chapter 6 


Suppression of the Background to 
QED Compton scattering 


In this chapter we perform a detailed study of the QED Compton scattering process 
(QEDCS) in ep e'yp and ep eyX, depicted in Figure 6.1, together with the ma¬ 
jor background coming from the virtual Compton scattering (VCS), where the photon is 
emitted from the hadronic vertex, shown in Figure 6.2. The two processes can be distin¬ 
guished experimentally because they differ in the kinematic distributions of the outgoing 
electron and photon. We suggest new kinematical cuts to suppress the VCS background, 
which turns out to be important in the phase space domain of the HERA experiment. We 
also study the impact of these constraints on the QEDCS cross section. We show that 
in the phase space region suggested and accessible at HERA, the photon content of the 
proton provides a reasonably good description of the QEDCS cross section, also in the 
inelastic channel. 

The VCS cross section in the inelastic channel is estimated utilizing effective parton 
distributions of the proton. In the elastic channel, to make a relative estimate of the VCS, 
we take the proton to be pointlike and replace the vertex by an effective vertex [4]. 

In Sections 6.1 and 6.2 we present the cross section in the elastic and inelastic channels 
respectively. Numerical results are given in Section 6.3. The summary is presented in 
Section 6.4. The matrix elements are explicitly shown in Appendix D. 
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+ Crossed. 


Figure 6.1: Feynman diagrams for the QED Compton proeess (QEDCS). X = p (and 
Px = P') corresponds to elastic scattering. 



Figure 6.2: As in Figure 6.1 but for the virtual Compton scattering (VCS) background 
process. 


6.1 Elastic Channel 

The elastic channel of the process under study, 

e{l)+p{P)^e{n+p{P')+^{P), (6.1) 

is described by the Feynman diagrams in Figures 6.1 and 6.2, with X being a proton 
and Px = P' ■ As in the previous chapter, we neglect the electron mass me everywhere 
except when it is necessary to avoid divergences in the formulae and take the proton to 
be massive, P^ = P'^ = m?. The corresponding cross section, using the notation (5.4) for 
the Lorentz invariant iV-particle phase-space element and the variables dehned in (4.39) 
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and in (5.2), is given by (5.5): 

^el(5) = ^ [ dPS2+l{l + P-, k', P')\ Mel I' . (6.2) 

2(s — m^) J 

Rearranging the (2 + l)-particle phase space into a sequence of two 2-particle ones as we 
did in Section 5.2, see (5.6), (5.7) together with (5.3), we get 

aei(s) = —^-^[^dPS2{l + P-,l' + k\P')dPS2{l + k-k,k')iM^, (6.3) 
2(s — m^) J 2n 

where, as before, k = P — P' is the momentum transfer between the initial and the hnal 
proton. Substituting (5.9) and (5.20) in (6.3) we obtain the hnal formula 




r-27r 


O-ellS) = 


2(47r) 


m 


2\2 


ds 


dt 


dt 


dcp* 


t) 


M. 


(6.4) 


where if* is the azimuthal angle of the outgoing e — 7 system in the e — 7 center-of-mass 
frame. The limits of integrations in (6.4) follow from kinematics and are given explicitly 
by (5.17) and (5.22). However, as it will be discussed in Section 6.4, we will impose 
additional kinematical cuts relevant to the experiment at HERA. 


As already shown in (5.3), the amplitude squared of the QED Compton scattering can 
be written as 

I j,,^QEDOS |2 _ ^ p,) _ (g 5) 

where Tafs is the leptonic tensor given by (5.8) and Pl^f is the hadronic tensor (5.18), 
expressed in terms of the electromagnetic form factors of the proton. The full cross section 
for the process given by (6.4) also receives a contribution from the VCS in Figure 6 . 2 . 
The cross section for this process can be expressed in terms of off-forward or generalized 
parton distributions [75,86]. In addition, there are contributions due to the interference 
between the QEDCS and VCS. In order to make a numerical estimate of these effects, 
one needs some realistic parametrization of the off-forward distributions. Our aim is to 
estimate the VCS background so as to hnd the kinematical cuts necessary to suppress 
it. We make a simplihed approximation to calculate the VCS cross section. We take the 
proton to be a massive pointlike fermion, with the equivalent 7 *p vertex described by a 
factor —i 7 "Fi(t). Incorporating the background effects, the cross section of the process 
in (6.1) is given by (6.4), where | M^x now becomes 

I Mel 1^ = I f + I |2 _ 2 gfJe (6.6) 
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The interference term will have opposite sign if we consider a positron instead of an 
electron. 

The explicit expressions of | f, \ f and 2 4ie are given 

in Appendix D.l. The effect of the proton mass is small in the kinematical range of 
HERA. 


6.2 Inelastic Channel 


We next consider the corresponding inelastic process, where an electron and a photon are 
produced in the final state together with a general hadronic system X: 

e{l) + p{P) ^ e(/') + j{k') + X{Px). (6.7) 

The exact calculation of the QEDCS rates follows our treatment in Chapter 5 based 
on the ALLM97 parametrization [70] of the proton structure function F 2 {xb,Q‘^)- For 
the purpose of evaluating the relative importance of the VCS background we resort to a 
unified parton model estimate of the VCS and QEDCS rates. The cross section within 
the parton model is given by 


dxB dQ^ ds dt d(p* 


q{xB,Q^" 


ds dQ"^ di d(p* 


and dd‘^ is the 


where q^XB^Q"^) are the quark and antiquark distributions of the initial proton, q = 
M, d, s, h, d, s. Furthermore, = —k"^ = —(/' + k' — /)^, xb = 2 pq-k) 
differential cross section of the subprocess 

e(/) + q(p) e(l') + -f(k') + q(p'), (6.9) 


which, similarly to (6.4), can be written as 


Ml |^ 


dsdQ^dtdv?* 2(47r)^(s - m2)2 (s - t) 

The relevant integrated cross section is obtained inserting (6.10) into (6.8), 


1 r'^max rQl, 

^mel(s) = TUl-TTf -/ 

2(47r)4(5 - m2)2 


dQ2 C'n 

^ L 


( 6 . 10 ) 


dt / dip* 

Jo 

( 6 . 11 ) 
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where we have traded the integration variable xb with 

= {p — kY = - — — . ( 6 . 12 ) 

Xb 

The limits of integration are given explicitly by (5.17), (5.41) and (5.42) with Smin = ''^e- 
Further constraints, related to the HERA kinematics, will be discussed in the numerical 
section. Similar to the elastic channel, we have 

I M<? 1^ = I M<?QEDCS 1^ I 7Vf'?vcs 1^ - 2 3fJeM'?QEDCS]pfgVCs*_ (6.13) 

Again, the interference term will have opposite sign for a positron. Furthermore, we 
introduce the auxiliary invariants S = {p' + k'Y and U = {p' — kY, which can be written 
in terms of measurable quantities, 

S = U = i-S + QY (6.14) 

Xl 

with Xl already dehned in (5.31). The explicit expression of (6.13) is deferred to Appendix 
D.2. It is also given in [87,88] for a massless proton. 

6.3 Numerical Results 

In this section the numerical results are presented. In order to select the QEDCS events, 
certain kinematical constraints are imposed in the Monte Carlo studies in [31,32]. As in 
Section 5.4 the following laboratory frame variables are used: energy of the hnal electron 
Eg, energy of the hnal photon polar angles of the outgoing electron and photon, 9e 
and 9.y respectively, and acoplanarity angle 0, which is dehned as 0 = | tt — |0.y — 0e| |, 
where 0.y and 0e are the azimuthal angles of the outgoing photon and electron respectively 
(0 < (j)y, 0e < 2 7r). The cuts are given in column A of Table 6.1 (from hereafter, they will 
be referred to as the set A). The energies of the incoming particles are: Eg = 27.5 GeV 
(electron) and Ep = 820 GeV (proton). 

So far the photon and the electron in the hnal state have been identihed only in the 
backward part of the HI detector at HERA. To select signals where there are no hadronic 
activities near the two electromagnetic clusters, the hnal hadronic state must not be found 
above the polar angle = 7r/2 [32]. Motivated by this experimental arrangement, we 
have identihed 9h with the polar angle of the hnal quark q' in the subprocess eg eyg'. 




86 


Suppression of the Background to QED Compton scattering 


A 

B 

E;, E; > 4GeV 

E' + E; > 20 GeV 

0.06 < 9e, 9.y < 71 — 0.06 

0 < 7r/4 

9h < T^l‘^ 

E', E; >4GeV 

E' + E; > 20 GeV 

0.06 < 9e, 9.y < 71 — 0.06 

s>Q^ 

S> s 


Table 6.1: A: cuts to simulate HERA-Hl detector. B: cuts introduced in this chapter. 

It can be shown that 9h is given by 

cos dh = cos 6qi = —— {xsEp — Ee — E^ cos 9e — E' cos 6*^) (6.15) 

Eg/ 

and Eqi = xsEp + Ee — E'^ — E'^ being the energy of the hnal parton. Here we have assumed 
that the hnal hadrons are emitted collinearly with the struck quark q'. For the elastic 
process 9h = 9pi, the polar angle of the scattered proton, can be obtained by substituting 
= 1 in the above expression. Thus we impose the additional condition [32] 

9h < 7r/2 (6.16) 

on the cross section. However, no constraint on the hadronic hnal state was used in the 
cross section calculation presented in [32]. Inclusion of (6.16) reduces the QEDCS cross 
section by about 10%. 

In the kinematical region dehned by the constraints mentioned above, the contributions 
from the initial and hnal state radiations, unrelated to QED Compton scattering, are 
suppressed, see Section 5.4 and [26,27,31,34]. Furthermore, we checked that the event 
rates related to the elastic VCS process and its interference with elastic QEDCS are 
negligible compared to the ones corresponding to pure elastic QEDCS. This is expected 
because the elastic QEDCS cross section is very much dominated by the small values of 
the variable —t, compared to —t, see (D.9) and (D.IO). Such an observation is similar to 
that of [31], where the elastic DVCS background was calculated using a Regge model in 
diherent kinematical bins. Our estimate was done taking the proton to be pointlike with 
an ehective vertex, as discussed in Section 6.1. We hnd that, in this approximation, the 
elastic QEDCS cross section dihers from the actual one calculated in Chapter 5 by about 
3% within the range dehned by the kinematical constraints. 
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Figure 6.3: Double differential cross section for QED Compton scattering at HERA-Hl. 
The kinematical bins correspond to Table 5.1. The continuous line describes the total 
(elastic + inelastic) cross section subject to the set of cuts A in Table 6.1. The dotted 
line shows the same results when the constraint on 9h is removed. 
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Figure 6.4: Cross section for the QEDCS and VCS processes (inelastic) at HERA-Hl. 
The bins are in s — S, expressed in GeV^. The cuts applied are listed in Table 6.1, set 
B (except S > s). The continuous line corresponds to the QEDCS cross section with 
ALLM97 parametrization of F 2 {xb,Q^), the dashed line corresponds to the QEDCS cross 
section using the effective GRV98 parton distributions in (6.17) and the dashed dotted line 
corresponds to the VCS cross section using the same effective distributions. 


Figure 6.3 shows the total (elastic + inelastic) QEDCS cross section in xi — bins 
with Qf = —t, subject to the cuts of set A. For comparison we have also plotted the 
cross section without the cut on 6h, similar to the analysis in Chapter 5. This additional 
constraint affects the result only in the inelastic channel. 

We checked that the upper limit in (6.16) reduces the contribution from the inelastic 
VCS reaction. In order to calculate it, one needs a model for the parton distributions 
q{xB,Q‘^)- However, in the relevant kinematical region, can be very small and may 
become close to zero, where the parton picture is not applicable. Hence, in our estimate, 
we replace the parton distribution q{xB,Q^) by an effective parton distribution [4] 

=-7z^^—7^q{xB,Q‘^+ Ql), (6.17) 

“r Cl 


where a = 1/4 and Ql = 0.4 GeV^ are two parameters and q{xB,Q‘^) is the NLO GRV98 
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Figure 6.5: Cross section for QED Compton scattering in bins of x.y as calculated with 
the ALLM97 (full line) and the (GRV98)gjj (dashed line) parametrization of F 2 {xb,Q^), 
respectively, as compared to the VCS background cross section (dot-dashed line). The cuts 
employed are: a) as in set A, b) as in set B of Table 6.1. The dotted line in Figure 6.5 
(a) shows the VCS cross section subject to the set of cuts A without the constraint on 6^. 


[68] parton distribution. Ql prevents the scale in the distribution to become too low. 
Equation (6.17) is motivated by a similar form used in [70,71] for the parametrization 
of the structure function F 2 {xb,Q'^) in the low region. It is clear that at high 
q{xB,Q‘^) q{xB,Q‘^)- 

In this chapter, we introduce a new set of cuts, which are given in the column B of 
Table 6.1 (and will be referred to as the set B) for a better extraction of the equivalent 
photon distribution of the proton as well as to suppress the VCS background. These 
cuts will be compared to the set A in the following. Instead of the constraint on the 
acoplanarity, namely 0 < 7r/4, where the upper limit is actually ambiguous, we impose 
s > Q^. The relevance of the cut S > s can be seen from Figure 6.4. This shows the 
cross sections of the QEDCS and VCS processes in the inelastic channel, calculated using 
(6.11) and subject to the kinematical limits of set B (except A > s), in bins of s — 
Figure 6.4 shows that the VCS cross section is higher than QEDCS for bins with s > S 
but falls sharply in bins for which s is close to S and becomes much suppressed for > s . 
This is expected because S corresponds to the quark propagator in the VCS cross section. 
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Xi 


Figure 6.6: Cross section for QED Compton scattering at HERA-Hl subject to the cuts 
of set B in Table 6.1, in (a) xi bins, (b) bins. The continuous line corresponds to 
our exact calculation using ALLM97 parametrization of F 2 {xb,Q‘^), the dot-dashed line 
corresponds to the same in the EPA, the dashed line shows the elastic contribution. 

see (D.29), and a lower value enhances this contribution. In fact the sharp drop of the 
VCS cross section in bins where ^ > s is due to the fact that both the propagators s, it 
in the QEDCS cross section are constrained to be smaller than S, U for VCS in these 
bins, see (D.28), (D.29). The QEDCS cross section is always enhanced by the factor 
in the denominator of (D.28) coming from the virtual photon, which can be very small in 
the kinematical region of interest here. This plot shows that imposing a cut on S can be 
very effective in reducing the background contribution from VCS. The interference between 
inelastic QEDCS and VCS gives negligible contribution. We have also shown the QEDCS 
cross section using the ALLM97 parametrization of F 2 {xb, Q^) as calculated in Chapter 5. 
The discrepancy between this and the one calculated using the parametrization in (6.17) 
is less than 5% in almost all the bins, and maximally 7% in two bins. 

In Figure 6.5 (a) we have shown the inelastic QEDCS and VCS cross sections in bins 
of x-f, dehned in (5.28), subject to the cuts of set A. The VCS cross section is much 
suppressed in the smaller x-^ bins but becomes enhanced as x-^ increases, which indicates 
that such a set of cuts is not suitable to remove the background at higher x.y. The situation 
will be the same in xi bins, with xi dehned in (5.31). 
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Figure 6.5 (b) shows the cross sections but with the set B. The background in this case 
is suppressed for all x.y bins, which means that such a cut is more effective in extracting 
QEDCS events also for higher x^. In addition, we have plotted the QEDCS cross section 
in terms of the structure function F 2 {xb, Q'^), using the ALLM97 parametrization. Figure 
6.5 shows that our parametrizaton gives a reasonably good description of the proton, at 
least for the QEDCS process, in all the bins. However, this parametrization has been 
used only to make a relative estimate of the background events. In fact, a quantitative 
estimate of the inelastic VCS events has not been presented in [31,32]. 

Eigures 6.6 (a) and 6.6 (b) show the QEDCS cross section in bins of xi and x.y, respec¬ 
tively, subject to the constraints of set B. The elastic cross section has been calculated 
using (6.4)-(6.5), as in Chapter 5. The inelastic cross section is given by (5.40) in terms 
of the structure functions Fi{xb, Q^) and F 2 {xb, Q^)- We have assumed the Callan-Gross 
relation and used the ALLM97 parametrization [70] for F 2 {xb, Q^)- In this way the results 
presented in Figure 6.6, labelled as exact, are free from the parton model approximations 
in Figures 6.4 and 6.5. In the same plot, we have also shown the total cross section cal¬ 
culated in terms of the EPA, according to (5.25) and (5.43). Figure 6.6 (b) shows much 
better agreement between the approximate cross section based on the EPA and the exact 
one. For Figure 6.6 (a), the discrepancy is about 3 — 8% in the hrst three bins, between 
10 — 20% in three other bins and 15% in the last bin. In Figure 6.6 (b) it is 1 — 6% in 
Eve bins, 12 — 13% in two bins and about 8% in the last bin. The discrepancy of the 
exact cross section, integrated over x.y, with the approximate one, when subject to the 
constraints of set B is 0.38% in the elastic channel and 4.5% in the inelastic one. The 
total (elastic -|- inelastic) discrepancy turns out to be 2.26%, which should be compared 
to the values 14%, already observed in Chapter 5 when subject to the set A, except the 
one on 9h, and 24% when this one is imposed too. 

As shown in Section 5.4, the elastic QEDCS cross section is described very accurately 
by the EPA. It is thus more interesting to investigate the inelastic channel in this context. 
The elastic QEDCS events can be separated from the inelastic ones by applying a cut on 
9h. We have found that, with the restriction 9h > 0.1°, the elastic events are rejected and 
all the inelastic events are retained in the cross section. A lower limit on 9^ higher than 
1° removes a substantial part (more than 30%) of the inelastic events. 

Table 6.2 shows the exact inelastic QEDCS cross section in xi and bins, subject to 
the cuts A. We have also shown the cross section in the EPA with the same constraints 
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xi bin 

Qf bin 

<7inel 

EPA 

^inel 

^inel 

^EPA* 

^inel 

1.78 X 10-5 - 5.62 X 10-5 

1.5 -2.5 

5.511 X 10^ 

1.483 X 102 

1.062 X 102 

1.298 X 102 

1.78 X 10-5 - 5.62 X 10-5 

2.5 -3.5 

1.992 X 10^ 

3.257 X lOi 

3.925 X lOi 

2.888 X lOi 

5.62 X 10-5 - 1.78 X 10-^ 

1.5-5.0 

1.720 X 102 

3.967 X 102 

2.937 X 102 

3.369 X 102 

5.62 X 10-5 - 1.78 X 10-^ 

5.0-8.5 

8.355 X 10^ 

2.015 X 102 

1.407 X 102 

1.764 X 102 

5.62 X 10-5 - 1.78 X 10-^ 

8.5 - 12.0 

2.609 X 10^ 

7.752 X lOi 

4.334 X lOi 

6.880 X lOi 

1.78 X 10-^ - 5.62 X 10-^ 

3.0 - 14.67 

1.613 X 102 

2.103 X 102 

2.330 X 102 

1.720 X 102 

1.78 X 10-^ - 5.62 X 10-^ 

14.67 - 26.33 

7.639 X 10^ 

1.477 X 102 

1.194 X 102 

1.283 X 102 

1.78 X 10-^ - 5.62 X 10-^ 

26.33 - 38.0 

2.269 X 10^ 

4.229 X lOi 

3.554 X lOi 

3.759 X lOi 

5.62 X 10-^ - 1.78 X 10-3 

10.0 - 48.33 

8.425 X 10^ 

7.555 X lOi 

9.953 X lOi 

5.980 X lOi 

5.62 X 10-^ - 1.78 X 10-3 

48.33 - 86.67 

3.745 X 10^ 

4.897 X lOi 

5.638 X lOi 

4.168 X lOi 

5.62 X 10-^ - 1.78 X 10-3 

86.67 - 125.0 

1.066 X 10^ 

1.462 X lOi 

1.644 X lOi 

1.253 X lOi 

1.78 X 10-3 - 5.62 X 10-3 

22 - 168 

3.846 X 10^ 

2.791 X lOi 

3.773 X lOi 

2.104 X lOi 

1.78 X 10-3 - 5.62 X 10-3 

168 - 314 

1.622 X 10^ 

1.849 X lOi 

2.289 X lOi 

1.543 X lOi 

1.78 X 10-3 - 5.62 X 10-3 

314 - 460 

5.836 

9.043 

7.827 

7.641 

5.62 X 10-3 - 1.78 X 10-2 

0 - 500 

1.202 X 10^ 

7.624 

9.281 

4.923 

5.62 X 10-3 - 1.78 X 10-2 

500 - 1000 

1.010 X 10^ 

1.450 X lOi 

1.242 X lOi 

1.190 X lOi 

5.62 X 10-3 - 1.78 X 10-2 

1000 - 1500 

2.584 

4.067 

3.033 

3.281 

1.78 X 10-2 - 5.62 X 10-2 

0 - 1500 

2.712 

1.439 

1.536 

7.261 X 10-1 

1.78 X 10-2 - 5.62 X 10-2 

1500 - 3000 

2.967 

3.855 

3.091 

2.922 

1.78 X 10-2 - 5.62 X 10-2 

3000 - 4500 

7.423 X 10-^ 

1.004 

7.210 X 10-1 

7.204 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

10 - 6005 

7.083 X 10-1 

4.923 X 10-1 

3.684 X 10-1 

2.354 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

6005 - 12000 

3.637 X 10-1 

3.845 X 10-1 

2.611 X 10-1 

2.124 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

12000 - 17995 

3.638 X 10-2 

3.847 X 10-2 

1.884 X 10-2 

1.537 X 10-2 


Table 6.2: Double differential QED Compton scattering cross section (inelastic) in xi and 
Qf bins. (Jinei and correspond to the exact cross section subject to the cuts A and B 
of Table 6.1 respectively, and correspond to the cross sections in the EPA and 

subject to the cuts A and B respectively. Qf is expressed in GeV^ and the cross sections 
are in pb. 
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bin 

bin 

^inel 

EPA 

^inel 

^inel 

^EPA* 

^inel 

1.78 X 10-5 - 5.62 X 10-5 

1.5 -2.5 

5.191 X 10^ 

1.483 X 102 

9.932 X lOi 

1.298 X 102 

1.78 X 10-5 - 5.62 X 10-5 

2.5 -3.5 

2.839 X 10^ 

3.257 X lOi 

5.176 X lOi 

2.888 X lOi 

5.62 X 10-5 - 1.78 X 10-^ 

1.5 -5.0 

1.761 X 102 

3.967 X 102 

2.992 X 102 

3.369 X 102 

5.62 X 10-5 - 1.78 X 10-^ 

5.0-8.5 

1.101 X 102 

2.015 X 102 

1.775 X 102 

1.764 X 102 

5.62 X 10-5 - 1.78 X 10-^ 

8.5 - 12.0 

4.573 X IQi 

7.752 X lOi 

7.189 X lOi 

6.880 X lOi 

1.78 X 10-^ - 5.62 X 10-^ 

3.0 - 14.67 

1.006 X 102 

2.103 X 102 

1.442 X 102 

1.720 X 102 

1.78 X 10-^ - 5.62 X 10-^ 

14.67 - 26.33 

9.299 X 10^ 

1.477 X 102 

1.347 X 102 

1.283 X 102 

1.78 X 10-^ - 5.62 X 10-^ 

26.33 - 38.0 

3.564 X 10^ 

4.299 X lOi 

5.166 X lOi 

3.759 X lOi 

5.62 X 10-^ - 1.78 X 10-3 

10.0 - 48.33 

4.408 X 10^ 

7.555 X lOi 

5.386 X lOi 

5.980 X lOi 

5.62 X 10-^ - 1.78 X 10-3 

48.33 - 86.67 

4.476 X 10^ 

4.897 X lOi 

6.082 X lOi 

4.168 X lOi 

5.62 X 10-^ - 1.78 X 10-3 

86.67 - 125.0 

1.611 X 10^ 

1.462 X lOi 

2.222 X lOi 

1.253 X lOi 

1.78 X 10-3 - 5.62 X 10-3 

22 - 168 

2.019 X 10^ 

2.791 X lOi 

2.174 X lOi 

2.104 X lOi 

1.78 X 10-3 - 5.62 X 10-3 

168 - 314 

1.896 X 10^ 

1.849 X lOi 

2.418 X lOi 

1.543 X lOi 

1.78 X 10-3 - 5.62 X 10-3 

314 - 460 

7.594 

9.043 

9.766 

7.641 

5.62 X 10-3 - 1.78 X 10-2 

0 - 500 

6.058 

7.624 

5.219 

4.923 

5.62 X 10-3 - 1.78 X 10-2 

500 - 1000 

1.077 X 10^ 

1.450 X lOi 

1.263 X lOi 

1.190 X lOi 

5.62 X 10-3 - 1.78 X 10-2 

1000 - 1500 

3.019 

4.067 

3.449 

3.281 

1.78 X 10-2 - 5.62 X 10-2 

0 - 1500 

1.267 

1.439 

7.874 X 10-1 

7.261 X 10-1 

1.78 X 10-2 - 5.62 X 10-2 

1500 - 3000 

3.005 

3.855 

3.018 

2.922 

1.78 X 10-2 - 5.62 X 10-2 

3000 - 4500 

8.121 X 10-1 

1.004 

7.650 X 10-1 

7.204 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

10 - 6005 

4.550 X 10-1 

4.923 X 10-1 

2.665 X 10-1 

2.354 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

6005 - 12000 

3.729 X 10-1 

3.845 X 10-1 

2.619 X 10-1 

2.124 X 10-1 

5.62 X 10-2 - 1.78 X 10-^ 

12000 - 17995 

4.083 X 10-2 

3.847 X 10-2 

2.117 X 10-2 

1.537 X 10-2 


Table 6.3: As in table 6.2 but for x.y bins. 
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(the last two cuts of set A are not relevant in this case). The discrepancy with the EPA 
is quite substantial. We have also shown the results with the cuts B, both the exact and 
the one in terms of the EPA, in the same table (the constraint s > is not relevant for 
the EPA). The discrepancy between the exact and the EPA here is much less and on the 
average it is 24%. 

Table 6.3 is almost similar, the only difference is that the bins are now in x.y. With 
the cuts of set A, the discrepancy now is on the average 44%, whereas, with the cuts B, 
the average discrepancy is 17%. 

Our results show that the extraction of the equivalent photon distribution 7 ^( 0 ;,/x^) is 
very much dependent on the kinematical constraints utilized to single out QEDCS events, 
in particular on the one on acoplanarity. The kinematical limits presented here are much 
more appropriate than those suggested in [31] for a reliable extraction of 7 ^’(a;,/i^). It is 
also clear that this discrepancy is entirely due to the inelastic channel. 

6.4 Summary 

To summarize, in this chapter we have analyzed the QED Compton process, relevant 
for the experimental determination of the equivalent photon distribution of the proton 
jj4). We have also calculated the major background process, namely virtual Compton 
scattering, assuming an effective parametrization of the parton distributions of the proton, 
both in the elastic and inelastic channels. The elastic VCS is suppressed compared to the 
QEDCS, in the phase space region accessible at HERA. We have shown that a constraint 
on the invariants S' > s is very effective in removing the inelastic VCS background. 
Furthermore, the selection of the QEDCS events in the process ep eyX is sensitive to 
the specihc kinematical limits, in particular to the upper limit of the acoplanarity angle 
0, which was used in the recent analysis [31, 32] of events as observed with the HERA-Hl 
detector. Instead of the acoplanarity, one can also directly impose cuts on the invariants, 
like s > (both of them are measurable quantities), which directly restricts one to 
the range of validity of the EPA. With these constraints, the total (elastic + inelastic) 
cross section agrees with the EPA within 3%. Thus, we conclude that by choosing the 
kinematical domain relevant for this approximation carefully, it is possible to have a more 
accurate extraction of 7 ^( 0 ;,/i^). 




Chapter 7 

The Polarized QED Compton 
Scattering Process 


This chapter, based on [5, 7], is devoted to the study the QED Compton scattering process 
in ip e'jp and ip —> eyX, where the initial lepton and proton are longitudinally 
polarized. We show that, when the virtuality of the exchanged photon is not too large, 
the cross section can be expressed in terms of the polarized equivalent photon distribution 
of the proton. We provide the necessary kinematical cuts to extract the polarized photon 
content of the proton at HERMES, COMPASS and eRHIC (the future polarized ep collider 
planned at BNL). In addition, we show that such an experiment can also access the 
polarized structure function at HERMES in the low region and at eRHIC 

over a wide range of the Bjorken scaling variable xb and Q^. 

The structure function Qi^xb-iQ"^) and its first xb moment in the low region have 
been studied in fully inclusive measurements at SLAC [89], HERMES [90,91] and JLab 
[92,93]. The most recent measurements by CLAS [94] are in the kinematical region 

= 0.15 — 1.64 GeV^. The low region is of particular interest because contributions 
due to nonperturbative dynamics dominate here and thus the transition from soft to hard 
physics can be studied. In fact the measurements in [94] clearly indicate a dominant 
contribution from the resonances and at higher they are below the perturbative QCD 
evolved scaling value of gi. This in fact illustrates the necessity of further investigation of 
gi{.XB, Q"^) in the transition region. In these fixed target experiments, low is associated 
with low values of thus the covered kinematical region is smaller compared to the 
unpolarized data. Data on gi{xB,Q‘^) for small xb and in the scaling region are missing 
due to the absence of polarized colliders so far (with the exception of RHIC, which has 
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started operating in the polarized mode for pp collisions only very recently). The small xb 
region is again interesting; it is the region of high parton densities, and measurements in 
this region will provide information about the effects of large [osln^ 1/xb]’^ resummation 
and DGLAP evolution, and also about the ’’soft” to ’’hard” scale transition [59, 95,96]. A 
better understanding of gi{xB,Q‘^) in this region is necessary in order to determine its hrst 
moment experimentally. The kinematics of QED Compton events is different from the 
one of inclusive deep inelastic scattering due to the radiated photon in the hnal state and 
thus it provides a novel way to access gi{xB, in a kinematical region not well covered 
by inclusive measurements, as already stated for F 2 {xb,Q‘^) in the preceding chapters. 

In Sections 7.1 and 7.2 we derive the analytic expressions of the cross section for 
the polarized QED Compton process in the elastic and inelastic channels, respectively. 
In Section 7.3 we discuss the background coming from virtual Compton scattering and 
also the interference between the two processes. The numerical results are presented 
in Section 7.4. A short summary is given in Section 7.5. The kinematics of the QED 
Compton process is described in Appendix D. The analytic expressions of the matrix 
elements can be found in Appendix E. 

7.1 Elastic QED Compton Scattering 

We consider QED Compton scattering in the elastic process: 

e{l)+p{P) ^ e{n+^{k')+p{P'), (7.1) 

where the incident electron and proton are longitudinally polarized and the four-momenta 
of the particles are given in brackets. Instead of the electron, one can also consider a 
muon beam (COMPASS); the analytic expressions will be the same. We make use of the 
invariants (4.39) and (5.2), k = P — P' is the four-momentum of the virtual photon. As for 
the unpolarized reaction, discussed in Chapter 5, we neglect the electron mass everywhere 
except when it is necessary to avoid divergences in the formulae and take the proton to 
be massive, P^ = P'^ = m^. The relevant Feynman diagrams for this process are shown 
in Figure 5.1, with X being a proton and Px = P'■ The matrix element squared can be 
written as 

I f = ^Hf\P, P')T^p{k k'), (7.2) 
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P') and I', k') being the antisymmetric parts of the hadronic and leptonic 

tensors respectively, which contribute to the polarized cross section. As before we use 
the notation (5.4) for the Lorentz invariant A^-particle phase-space element. The spin 
dependent counterpart of (5.6) reads 

1 r d s 1 

J — dPS,(l + P; r + k', P')A'* (i, *:) . (7.3) 

The tensor is antisymmetric in the indices a, f3 and is dehned as 

XUr *;) = / dPS^d + k-, ((, k; I', k'), (7.4) 

where k; k') is the antisymmetric part of the leptonic tensor (B.31), 

Aip^ r 

T„^^(/, k- k') = —{s - t)E + {u - t)l'P kC (7.5) 

sn L J 

contains all the informations about the leptonic part of the process and can be 
expressed in terms of the Lorentz scalar AX: 

AX(s;t), (7.6) 

with 


A.Y(s,() = -2AyP“'’*, (7.7) 

being the antisymmetric part of the photon polarization density matrix given in 
(3.54)-(3.55). We dehne the function AX(s,f,t) as 


AX(s,t) = 


di / dip* AX{s,t,i)] 


the integration limits of t are the same as in (5.17) and ip* is the azimuthal angle of the hnal 
electron-photon system in the electron-photon center-of-mass frame. Hence AX{s,t,i) 
can be obtained from the relation AX = —2 explicitly 


AX{s,t,i) = 


su{s — t) 




The hadronic tensor for polarized scattering (3.46), i.e. 

= -iehm-'X-khGEGMS, - - Ge) k ■ S 

I + T 


—P 
2 ”■ ’ 


(7.10) 




98 


The Polarized QED Compton Scattering Process 


is expressed in terms of the electric and magnetic form factors of the proton Ge and Gm', 
moreover r = —t/Aim? and S is the spin vector of the proton, which fulhls S ■ P = Q and 
= — 1. If we express S as 

S^ = Ns{Pc.-^l}\, (7.11) 


S^ = Ns\Pe.-j^Ly ( 

with Ns given in (3.50) and (3.51), then from (7.3), (7.6) and (7.8)-(7.10) we have 




max / t-max 


min ^ ‘'mm 


dt / d(p* AX{s,t,i) 

Jo 


( s — m^ 2m^ s — t \ , , 

x[(2—+ — 

^ / g - ^ s-t \ Gm{Gm - Ge) 

\ s — t t s — m? ) 1 + T 


(7.12) 


the integration bounds being the same as in (5.17) and (5.22). These bounds are modihed 
due to the experimental cuts which we impose numerically. In the EPA limit, we neglect 
I t I vs. s and w? vs. s and get 


A -N A -A 4q;^ /s u\ 2s dAd 

AX s,t,t ^ AX s,0,t = — =- 

S \U si TT dt 


(7.13) 


where dAd/dt is the polarized differential cross section for the real photoproduction pro¬ 
cess ey —> ey (B.37) and AX(s, 0) = —4sAd. The elastic cross section then becomes 


Aaei ^ Aa,f ^ = 




dt Ay^i(a; 


dAd(a;s, i) 


(7.14) 


where mg is the mass of the electron and Ayg[(a;) is the elastic contribution to the polarized 
equivalent photon distribution of the proton (3.58), 

a 2mV5^ 2 o/^i , GuiGu - Ge)' 

Ay,(x) = - -(^2- x+—jG„- 2(^i-x+— j ——j 

(7.15) 


2 — X + 


G — 2 1 — X + 


with X = s/s, tmin = —oo and tmax = —m^a;^/(l — x). 


7.2 Inelastic QED Compton Scattering 


We next consider the corresponding inelastic process 


e(;) + p(F) ^ e(r) + + X(Px). 


(7.16) 
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We take the invariant mass of the produced hadronic system to be hh; moreover = 
—k‘^ = —t and the Bjorken variable xb is related to them via (5.38). The cross section 
for inelastic scattering reads 


( 5 ) — 


327r3(s - m2)2 7,^2, 


^{VS-Wf pQL 
ds / 


(7.17) 


where is given by (7.6) and IT^ is the hadronic tensor, see (2.57) together with 
(2.69) and (2.70), 


= 2nie^—— SafSpak^ 9 i{xb, Q‘^)S„ + g 2 {xB, Q^) - -— -P, 


(7.18) 


Hence the cross section takes the form 


A(Tin»l(s) = 


47r(s - w?Y Jw2, 


’‘(CS-mf /-Qmax (^Q‘2 

ds / —-r- 


2. g2 + 


s — m? W'^ + Q‘^ — m? 2w? f s + 

s + g 2 Q 2 g 2 _ .^2 


giixB.Q"^ 


W2 + g2 _ ^2 


g2ixB:Q ) fAX(s,g ), 


(7.19) 


and AX^SjQ"^) is given by (7.8) and (7.9). The limits of the Q^, and t integrations 
are given in (5.41), (5.42) and (5.17) respectively. In the region of validity of the EPA, 
s 3> wS and s ^ the cross section becomes 


AcTinel ~ A(t)!!'i^ 


’(l-m/pS)2 


dt A7f^gi(a;,a;s) 


Ad(j(a;s, i) 


(7.20) 


where again x = s/s and A 7 ?g[(a;, xs) is the inelastic contribution to the polarized equi¬ 
valent photon distribution of the proton; 


^7fnel(a^i XS) 


a dy 7*5max (Iq 2 / 2wSx‘^ 

27rX y Jq^. V ^ 

^ min 


2(71 



(7.21) 


where we take the scale Q^ax be s and g^in = x‘^m‘^/{l — x). Here we have neglected 
the contribution from g 2 {xB,Q‘^)- Expressing gi^XB^Q"^) in terms of the polarized quark 
and antiquark distributions, one can conhrm that the above expression reduces to that 
given in (3.63). However, in this case, one chooses the minimal (but not compelling) 
boundary condition A 7 ?^gj(a;, gg) = 0 at a scale Q‘1 = 0.26 GeV^. The expression (7.21) 
is free from this particular boundary condition. 
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7.3 Background from Virtual Compton Scattering 


The processes ep e'yp and ep e'yX receive contributions from the virtual Compton 
Scattering (VCS), when the photon is emitted from the proton side as well as from the 
interference between the QED Compton scattering (QEDCS) and VCS, see Figures 6.1 
and 6.2. The polarized cross section for the elastic process, analogously to (6.4), reads 


A(Tel(s) 


1 

2(47r)"^(s 



(7.22) 


where 

I AMei 1^ = I f + I f -2?R:e (7.23) 


is the matrix element squared of the process. The integration bounds are the same as in 
(7.12). The interference term will have opposite sign if we consider a positron instead of 
an electron. The cross section of the VCS process is expressed in terms of generalized 
parton distributions and one needs a realistic model for a quantitative estimate of this 
background [86]. Here, as in the preceding chapter, in order to find the cuts to suppress 
the VCS, we make a simplihed assumption; we take the proton to be a massive pointlike 
particle with an effective 'y*p vertex, —i'y°‘Fi{t). The explicit expressions for the matrix 
elements are given in Appendix E. 


Particularly interesting for our purpose of extracting the polarized photon distribution 
of the proton is the inelastic channel. Here we use a unihed parton model to estimate the 
VCS and QEDCS rates, similar to (6.8). The cross section within the parton model is 
given by 


dAcTinel 

dxB dQ'^ ds didif* 


E 


Aq{xB,Q'^) 


dAd^ 

ds dQ2 di dip* 


(7.24) 


where Aq{xB, Q^) are the polarized quark and antiquark distributions of the initial proton, 
q = u, d, s, u, d, s and dAd"^ is the differential cross section of the subprocess 


e(/) + i[{p) e(t) + pik') + g{p'), 


(7.25) 


g being a longitudinally polarized quark in a longitudinally polarized proton and q a quark 
in the hnal state. The integrated cross section reads 
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HERMES 

COMPASS 

eRHIG 

Ee = 27.5 GeV 

Ef, = 160 GeV 

Ep = 250 GeV 

Ee = 10 GeV 

0.04 < Oe, 9-f < 0.2 

0.04 < Of,, 9-^ < 0.18 

0.06 < 9e, 9-f < 71 — 0.06 

E', E; >4GeV 

E;, E;>4GeV 

E;, E; >4GeV 

s > 1 GeV^ 

s > 1 GeV=^ 

s > 1 GeV^ 

s > Q‘^ 

s > Q‘^ 

s > Q‘^ 


Table 7.1: Energies, angular acceptance and kinematical cuts for the HERMES, COM¬ 
PASS and eRHIC experiments. 


^<^inel('S) — 






2(47r)^(s - ^ Jw2^ 


dW^ 


rQ 

Iq^. 

^ mi 


2 

max 


r*27r 


d* / f rN I d,- 

2 Q Jo 


X 


(s + QY 




(7.26) 


where 




— I A QEDCS 


+ 




vcs 


- 2 


(7.27) 


and the limits of integrations are given in (5.17), (5.41) and (5.42). The explicit expres¬ 
sions of the matrix elements are given in Appendix E. As for the study of the unpolarized 
process, it is useful to introduce the invariants .S' = (p' -|- k'Y and U = {p' — kY, which 
can be written in terms of measurable quantities and satisfy the relation (6.14). 


7.4 Numerical Results 

In this section we present our numerical results. The cuts used for the HERMES, COM¬ 
PASS and eRHIC kinematics are given in Table 7.1. As for the spin independent process, 
the constraints on the energies and polar angles of the detected particles reduce the back¬ 
ground contributions coming from the radiative emissions (when the hnal state photon 
is emitted along the incident or the hnal lepton line), because they prevent the lepton 
propagators to become too small. The QED Compton events are singled out at HERA 
by imposing a maximum limit on the acoplanarity angle 0 dehned in (5.47). We have 
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observed in the preceding chapter that, instead of this limit on 0, the constraint s > Q^, 
which is applicable experimentally, is more efficient in extracting the equivalent photon 
distribution from the exact result. Here we use this constraint. 


The unpolarized cross section has been calculated using the formulae in Chapter 4; for 
the numerical estimates we have used the ALLM97 parametrization [70] of the structure 
function F 2 {xb,Q‘^)- We have taken Fl^xbjQ"^) to be zero, assuming the Callan-Gross 
relation. In the polarized cross section, we have neglected the contribution from g 2 {xB, 
and used the parametrization [71] for gi{xB,Q'^)- In this parametrization, gi^XBiQ"^) is 
described in the low-Q^ region by the generalized vector meson dominance (GVMD) model 
together with the Drell-Hearn-Gerasimov-Hosoda-Yamamoto sum rule and, for large Q^, 
giixB,Q^) is expressed in terms of the NLO GRSVOl [69] parton distributions (standard 
scenario) in terms of a suitably dehned scaling variable 


X = 


Q^ + Ql 

Q^ + Ql + W^- m? 


(7.28) 


with Ql = 1.2 GeV^. The scale Q‘^ is changed to Q'^ + Qo, so as to extrapolate to the low- 

region. It is to be noted that for QED Compton scattering, the effects of Fb{xb,Q^) 
and g 2 {xB, Q^) have to be taken into account in a more accurate study as their effect may 
become non-negligible in the low-Q^ region. However, this is beyond the scope of the 
present work. 

Before discussing the results for specihc experiments, it is interesting to investigate 
some general properties of the total cross section. Figures 7.1 (a) and (b) show the total 
QEDCS cross sections, polarized and unpolarized respectively, as functions of the incident 
lepton energy Ei. We have imposed the constraints in the second column of Table 7.1 on 
the energies and angles of the outgoing particles, as well as those on s. Both polarized and 
unpolarized cross sections increase sharply with reach a peak at around Ei = 20 GeV 
and then start to decrease. The cross sections in the inelastic channels are also shown, 
which have similar trends except that the peak in the polarized case is broader. 


7.4.1 HERMES 

Figures 7.2 (a) and (b) show the total (elastic+inelastic) polarized and unpolarized QED 
Gompton scattering cross sections respectively, in bins of the variable x^, defined in (5.28), 
for HERMES kinematics, subject to the cuts of Table 7.1. We have taken the incident 
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Ej (GeV) 


20 40 60 80 100 120 140 160 180 200 
Ej (GeV) 


Figure 7.1: QEDCS cross section vs. energy of the incident lepton; (a) polarized, (b) 
unpolarized. The continuous line is the total cross section and the dashed line is the 
cross section in the inelastic channel. The cuts imposed are given in the central column 
of Table 7.1. We have used the ALLM97 parametrization of F 2 and the Badelek et al. 
parametrization of gi. 


electron energy to be Ef. = 27.5 GeV. We also show the cross section calculated in the 
EPA. Furthermore the contributions due to the inelastic channel of the reaction are plot¬ 
ted. The cross section, integrated over Xy, agrees with the EPA within 7.1% (unpolarized) 
and 4.8% (polarized). From the hgures it is also clear that the agreement in the inelastic 
channel (2.5% in the polarized case) is much better than for HERA kinematics discussed 
in Chapter 5. This is because at HERMES can never become too large (maximum 13.7 
GeV^), subject to our kinematical cuts, which is expected in a hxed target experiment. 
The agreement is not so good without the constraint s > 1 GeV^. Figure 7.2 (c) shows 
the asymmetry, which is dehned as 

All = ~ (7.29) 

(T++ -h (T+_ 

where -|- and — denote the helicities of the incoming electron and proton. They are 
calculated with the same set of constraints. The asymmetry is quite sizable at HERMES 
and increases in higher Xy bins. The asymmetry in the EPA is also shown. It is interesting 
to note that the discrepancy between the exact cross section and the one evaluated in the 
EPA, evident in Figures 7.2 (a) and (b), actually gets canceled in the asymmetry; as a 
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Figure 7.2: Cross section for QED Compton scattering (QEDCS) at HERMES in bins 
of (a) polarized, (h) unpolarized, (c) the asymmetry; for the polarized cross section 
the Badelek et al. parametrization of gi (BKZ) and for the unpolarized cross section the 
ALLM97 parametrization of F 2 have been used; (d) polarized inelastic cross section for 
QEDCS (long dashed), VCS (dashed-dotted) and the interference (dashed) at HERMES in 
the effective parton model. The bins are in s — S, expressed in CeV^. The continuous line 
is the QEDCS cross section using the BKZ parametrization of gi{xB, Q^)- The constraints 
imposed are given in Table 7.1. 
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consequence Figure 7.2 (c) shows an excellent agreement, except in the last bin, between 
exact and approximated results. We have also calculated the expected statistical errors 
for each bin, using the following formula, valid when the asymmetry is not too large [97], 



where Ve and Vp are the polarizations of the incident lepton and proton, respectively, C is 
the integrated luminosity and aun is the unpolarized cross section in the corresponding 
bin. We have taken Ve = Vp = 0-7 and C = 1 fb“^ for HERMES. The expected statistical 
errors increase in higher bins, because the number of events become smaller. However 
the asymmetries seem to be measurable at HERMES. 


The background from virtual Compton scattering is reduced at HERA by the ex¬ 
perimental condition of no observable hadronic activity at the detectors. Basically the 
electron and photon are detected in the backward detectors and the hadronic system in 
the forward detectors. In the previous chapter, we have observed that for unpolarized 
scattering at HERA, such a constraint is insufficient to remove the VCS contribution for 
higher x^. We have proposed a new constraint S > s, where S and s can be measured 
experimentally, to be imposed on the cross section. Here, we investigate the effect of this 
constraint on the polarized cross section. To estimate the inelastic contribution coming 
from VCS, we use (7.26), together with an effective model for the parton distribution of 
the proton. The effective parton distributions are of the form 


Aq{xB, Q^) = Aq{x, + Ql), 


(7.31) 


Aqi^XBiQ"^) being the NLO GRSVOl (standard scenario) distribution functions [69]. In 
the relevant kinematical region, can be very small and may become close to zero, 
where the parton picture is not applicable. The parameter Qg = 2.3 GeV^ prevents the 
scale of the parton distribution to become too small, while h is a suitably dehned scaling 
variable, 

- ^ xb^Q"^ + Ql) 

+ XbQI 



To estimate the unpolarized background effect, we use the same expressions as in Chap¬ 
ter 6 with the effective parton distributions given in (6.17). Figure 7.2 (d) shows the 
polarized cross section in the inelastic channel at HERMES, subject to the constraints of 
Table 7.1, in bins of s — F calculated in the effective parton model. The VCS and the 
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Figure 7.3: Asymmetry in the inelastic channel in bins of xb at HERMES. We have used 
the Badelek et al. parametrization of gi. The constraints imposed are as in Table 7.1 
(except s > Q"^), together with S — s > 2 GeW. The average (in GeW) of each bin is 
also shown. 



interference contributions are also shown. The QEDCS cross section using the Badelek et. 
al parametrization of gi{xB,Q^) is also plotted. In fact, the cross section in the effective 
parton model lies close to this. Within the parton model, the VCS is suppressed when 
s < S, similar to the unpolarized case at HERA. Unlike HERA, the interference between 
QEDCS and VCS is not negligible at HERMES, although smaller than the QEDCS in 
the relevant region. Since the interference term changes sign when a positron beam is 
used instead of an electron beam, a combination of electron and positron scattering data 
can eliminate this contribution. In order to estimate the VCS in the elastic channel, one 
needs a suitable model for the polarized generalized parton distributions. However, in 
the simplihed approximation of a pointlike proton with an effective vertex as described in 
Section 7.3, the elastic VCS as well as the interference contribution is much suppressed 
at HERMES. Similar observations hold for unpolarized scattering. 

Eigure 7.3 shows the asymmetries in the inelastic channel in bins of xb- In addition to 
the cuts mentioned above and shown in table 7.1, we have also chosen S — s > 2 GeV^ to 
suppress the background. The asymmetry is small but sizable and could be a tool to access 
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gi{xB,Q‘^) at HERMES. In fact, QED Compton events can be observed at HERMES in 
the kinematical region xb = 0.02 — 0.7 and = 0.007 — 7 GeV^ (small Q^, medium xb)- 
However, from the figure it is seen that the asymmetry is very small for xb below 0.1. We 
have also shown the expected statistical error in each bin. The average value in GeV^ 
for the polarized cross section for each bin is shown, which has been calculated using the 
formula 


(Q^) 


/bin^'dAg 

/bin^Ag 


(7.33) 


7.4.2 COMPASS 

Figures 7.4 (a) and (b) show the cross sections of the polarized and unpolarized QEDCS 
process in bins of x^ for the kinematics of COMPASS. We take the energy of the incident 
muon beam to be 160 GeV, the target is a proton. The hnal muon and the photon are 
detected in the polar angle region 0.04 < 6^^, 9^ < 0.18. The cross sections in bins, subject 
to the kinematical constraints shown in Table 7.1, are much smaller than at HERMES, 
because they start to decrease with the increase of the incident lepton energy Ei as Ei 
becomes greater than about 20 GeV, as depicted in Figure 7.1. As before, the cuts 
remove the initial and hnal state radiative events. The x^ integrated cross section agrees 
with the EPA within 14.2% (unpolarized) and 15.5% (polarized). The agreement thus 
is not as good as at HERMES. From the hgures it is seen that the cross section in the 
EPA actually lies below the exact one, both for polarized and unpolarized cases. This 
discrepancy is due to the fact that the EPA is expected to be a good approximation when 
the virtuality of the exchanged photon is small. At COMPASS, with our kinematical 
cuts, can not reach a value below 0.07 GeV^ and can be as large as 144 GeV^, whereas 
for HERMES smaller values of are accessible (see the previous subsection). Figure 7.4 
(c) shows the asymmetry in bins of x.y, also in its inelastic channel. The asymmetry is 
of the same order of magnitude as in HERMES and is in good agreement with the EPA. 
We have also shown the expected statistical error in each bin, calculated using (7.30). 
We have taken Ve = Vp = 0.7 and C = 1 fb~^ for COMPASS. The statistical error is 
large in higher x^ bins. Figure 7.4 (d) shows the polarized QEDCS, VCS and interference 
contributions (inelastic) calculated in the effective parton model, in bins of s — S'. As 
in HERMES, VCS is suppressed for s < S. The interference term is not suppressed but 
using and fj,~ beams this can be eliminated. We have also plotted the QEDCS cross 
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Figure 7.4: (a), (b), (c) and (d) are the same as in Figure 1.2 but for COMPASS. The 
eonstraints imposed are given in Table 7.1. 
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Figure 7.5: Asymmetry in the inelastic channel in bins of xb at COMPASS. We have 
used the Badelek et al. parametrization of gi. The constraints imposed are as in Table 
7.1 (except s > Q'^), together with S — s > 2 GeW. The average (in GeW) of each 
bin is also shown. 

section using Badelek et al. parametrization of gi{xB, Q^)- The VCS and the interference 
contributions (elastic) are much suppressed in the pointlike approximation of the proton 
with the effective vertex. 

Figure 7.5 shows the asymmetry at COMPASS in the inelastic channel plotted in bins 
of Xb with the same set of constraints and the additional cut S — s > 2 GeV^. The 
asymmetry is sizable and can give access to gi{xB, Q'^), the kinematically allowed range is 
Xb > 0.07 . We have also shown the expected statistical errors in the bins and the average 
in each bin. Confronting Figure 7.3 and 7.5 one can see that there is no overlap in the 
kinematical region covered at HERMES and COMPASS. Higher values of are probed 
at COMPASS in the same xb range as compared to HERMES. 



7.4.3 eRHIC 

The cross sections for eRHIC kinematics, both polarized and unpolarized, are shown in 
Figure 7.6 (a) and (b) respectively, in bins of x.y. We have taken the incident electron 
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Figure 7.6: (a), (h), (c) and (d) are the same as in Figure 1.2 but for eRHIC. The 
eonstraints imposed are given in Table 7.1. 
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Figure 7.7: Asymmetry in the inelastic channel in bins ofxs at eRHIC. We have used the 
Badelek et al. parametrization of gi. The constraints imposed are as in Table 7.1 (except 
s > Q^), together with S > s. The average (in GeW) of each bin is also shown. 

energy Ee = 10 GeV and the incident proton energy Ep = 250 GeV. The cross section 
in the EPA is also shown. The kinematic constraints are given in Table 7.1. The polar 
angle acceptance of the detectors at eRHIG is not known. We have taken the range of 
9e, 0.y to be the same as at HERA. We have checked that the constraints on the energies 
and the polar angles of the outgoing electron and photon are sufficient to prevent the 
electron propagators to become too small and thus reduce the radiative contributions. 
The unpolarized total (elastic+inelastic) cross section, integrated over x.y agrees with the 
EPA within 1.6%. The agreement in the inelastic channel is about 6.3%. The polarized 
total cross section agrees with the EPA within 9.8%. The EPA in this case lies below the 
exact result in all the bins. The agreement in the inelastic channel is about 19.6%. More 
restrictive constraints instead of s > Q^, like s > 10 Q^, makes the agreement better, 
about 1.2% in the polarized case and 1.9% in the unpolarized case. Figure 7.6 (c) shows 
the asymmetry for eRHIG, in bins of x.y. The discrepancy in the cross section cancels in 
the asymmetry, as a result good agreement with the EPA is observed in all bins except 
the last one at higher x.y. The asymmetry in the inelastic channel is also shown. We have 
plotted the expected statistical error in the bins using (7.30). For eRHIG, we have taken 
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Figure 7.8: Asymmetry in bins of (GeV^) at eRHIC. We have used the Badelek et al. 
parametrization of gi. The constraints imposed are the same as in Figure 7.7. 

Re = Rp = 0.7 and C = 1 fb~^. The expected statistical error increases in higher bins. 
The asymmetry is very small for small but becomes sizable as x^ increases. Figure 
7.6 (d) shows the polarized cross section in the inelastic channel, in bins of s — F, in the 
effective parton model for eRHIC. The VCS is suppressed in all bins, especially for s < S. 
The interference contribution is negligible, similar to HERA. The effective parton model 
QEDCS cross section is also compared with the one calculated using the Badelek et. al. 
parametrization for Q^). Similar effects are observed in the unpolarized case. In 

the pointlike approximation of the proton with the effective vertex, as before, the elastic 
VCS as well as the interference contributions are very much suppressed. 

Figure 7.7 shows the asymmetry in bins oi xb in the inelastic channel, which may be 
relevant for the determination of qiIxbtQ"^) using QEDCS at eRHIC. The asymmetry is 
small but sizable, however the error bars are large and therefore good statistics is needed. 
Xb can be as low as 0.002. A wide range of can be accessed at eRHIC starting from 
0.008 to 2000 GeV^; the average value in the bins ranges from 2.4 to 315 GeV^. 

Figure 7.8 shows the total asymmetry in bins for eRHIC. The asymmetry in this 
case is bigger in each bin and the error bars are smaller than for the xb bins, except in 
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the last Q^-hin where the number of events is smaller. 

7.5 Summary 

In this chapter we have analyzed the QED Compton process in polarized Ip scattering, 
both in the elastic and inelastic channel. As for the unpolarized process, we have shown 
that the cross section can be expressed in terms of the equivalent photon distribution of 
the proton, convoluted with the real photoproduction cross section. Furthermore we have 
provided the necessary kinematical constraints for the extraction of the polarized photon 
content of the proton by measuring the QED Compton process at HERMES, COMPASS 
and eRHIC. We have shown that the cross section and, in particular, the asymmetries are 
accurately described by the EPA. We have also discussed the possibility of suppressing 
the major background process, namely the virtual Compton scattering. We point out that 
such an experiment can give access to the spin structure function gi{xB, Q^) in the region 
of low and medium in hxed target experiments and over a broad range oi xb, 
at the future polarized ep collider, eRHIC. Because of the different kinematics compared 
to the fully inclusive processes, the QED Compton process can provide information on 
gi{xB-, Q^) in a range not well-covered by inclusive measurements and thus is a valuable 
tool to have a complete understanding of the spin structure of the proton. 





Chapter 8 

uW Production in ep —^ uWX 


In [24] the unpolarized elastic photon distribution was tested in the case of uW production 
in the process ep uWp. The relative error of the cross section as calculated in the EPA 
with respect to the exact result was shown as a function of ^/s, in the range 100 < ^/s < 
1800 GeV. The agreement turned out to be very good, the approximation reproducing 
the exact cross section within less than one percent. Motivated by this results, following 
the lines of [6], our aim here is to check if the same holds in the inelastic channel. 

The process ep uWX has been widely studied by several authors [98-102]. Its 
relevance is related to the possibility of measuring the three-vector-boson coupling hPIPy, 
which is a manifestation of the non-abelian gauge symmetry upon which the Standard 
Model is based. The observation of the vector boson self interaction would be a crucial 
test of the theory. Furthermore, such a reaction is also an important background to a 
number of processes indicating the presence of new physics. The lightest Supersymmetric 
Standard Model particle has no charge and interacts very weakly with matter; it means 
that, exactly as the neutrino from the Standard Model, it escapes the detector unobserved 
and can be recognized only by missing momentum. This implies that a detailed study of 
the processes with neutrinos in the hnal states is necessary to distinguish between the new 
physics of the Supersymmetric Standard Model and the physics of the Standard Model. 
At the HERA collider energies = 318 GeV) the ep uWX cross section is much 
smaller than the one for ep eWX [99,100], also sensitive to the hPIEy coupling, due 
to the presence in the latter of an additional Feynman graph where an almost real photon 
and a massless quark are exchanged in a w-channel conhguration (w-channel pole). The 
dominance of the process ep eWX justihes the higher theoretical and experimental 
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[103] attention that it has received so far, as compared to ep uWX. One way of 
improving the problem of the low number of deep inelastic uW events at HERA would be 
to consider also the elastic and quasi-elastic channels of the reaction, as will be discussed 
in Section 8.2. 

It is worth mentioning that not all the calculations of the ep uWX event rates 
available in the literature, in which only the photon exchange is considered (see Figure 
8.1), are in agreement, as already pointed out in [102]. In particular, the numerical 
estimate of the cross section for HERA energies presented in [98,99], obtained in the EPA 
approach, is one half of the one published in [102], obtained within the framework of the 
helicity amplitude formalism without any approximation. The value given in [101] is even 
bigger than the one in [102]: all these discrepancies cannot be due to the slightly different 
kinematical cuts employed in the papers cited above and stimulate a further analysis. 
Our results agree with [98,99]. 

The plan of this chapter is as follows. In Section 8.1 we calculate the exact cross section 
for the inelastic channel in a manifestly covariant way and we show in which kinematical 
region it is supposed to be well described by the EPA. The formulae for the corresponding 
elastic cross sections, both the exact and the one evaluated in the EPA, are also given. 
The numerical results are discussed in Section 8.2. The summary is given in Section 8.3. 


8.1 Theoretical Framework 

The uW production from inelastic ep scattering, 

e(/) + p{P) ^ z/(/') + W{k^) + X{Px), (8.1) 

is described, considering only one photon exchange, by the Feynman diagrams depicted in 
Figure 8.1. The four-momenta of the particles are given in the brackets; Px = Xlx 
the sum over all momenta of the produced hadronic system. We introduce the invariants 

s={P + lY, s = {l + ky, i={l-l'Y, Q‘^ = -k‘^, (8.2) 

where k = P — Px is the four-momentum of the virtual photon. If we denote with e and 
the polarization vectors of the incoming photon and outgoing W and with u and u the 
Dirac spinors of the initial electron and hnal neutrino, the amplitude of the subprocess 
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Figure 8.1: Feynman diagrams for the process ep uWX. 


67 * —uW reads 


M = 


2v/2sin6'vK 


e^{k)e*^{k')u{r) 






(1 - 7")«.(/) 

(8.3) 


where is the mass of the W boson, 6w the weak-mixing angle and 

^ap^l = gap{‘2k - A;% + - k)^ - gpa{k' + k)p (8.4) 

describes the hFhF 7 vertex. 

The integrated cross section of the fnll reaction ep vWX has been calculated in [6] 
and can be written as (5.40), that is 
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(8.5) 


where IF^ indicates the invariant mass sqnared of the prodnced hadronic system X, ip* 
denotes the azimnthal angle of the outgoing v — W system in the v — W center-of-mass 
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frame, and 


Xb = 


+ Q2 _ ^2 


is the Bjorken variable. Furthermore, F’l, 2 ( 2 ^ 5 , are the structure functions of the 

proton and the two invariants Xi^ 2 {s,Q‘^,t), which contain all the information about the 
subprocess 67 * — uW, are given by (5.14), (5.15): 


Xl(s,g^^) = - 


X2{s,Q\i) = 


167r2(s + g^) 


The tensor VV“^ is obtained from (8.3): 


yy;a/3 _ _ ^ 


the sums over the electron and vector boson spins are performed by making use of the 
completeness relations (A. 17) and 


{k') — -gap + 


The hnal result can be expressed as 

xps.QM) = _ i-f, [(«" + + * + i) 


+ 2(g^ + s)^i + 8(g^ + s)P + si^ 


( 8 . 10 ) 


X2(s,g^^) = 


^-^{4M® (g2 + s) - 4M® [3s(g2 + s) 

2V27Tp{Q^ + s){M^-ty^ ^ ; 

+ (2g^ + s)i\ + 4iW^[s(2s + t)^ + g^(4.§^ + 2si + t^)] — M^s 
X [g^s + g^ (9.5^ + 2st — 4t^) + 4(s + t) (2.5^ + 2si + P)] 

+ g^s^[s(s + 1 ) + g^(s + 2 t)]}, (8.11) 


Gf = 


8 M^ sin^ 9i 


( 8 . 12 ) 
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In (8.5) the minimum value of s is given by the squared mass of the W boson: 


■Smin ’ 

while the limits of the integration over are: 

W"min = {rn + = (v^ - )^ 

where is the mass of the pion. The limits Q^in max given by: 
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(8.15) 


(8.16) 


Integrating Xi^ 2 (>s, t) over cp* and t, with the limits in (8.16), one recovers (4.1) and 
(4.2) of [24] respectively, times a factor of two due to a different normalization. The EPA 
consists of considering the exchanged photon as real; it is possible to get the approximated 
cross section from the exact one, (8.5), in a straightforward way, following again 
Chapter 5. We neglect compared to s and compared to s then, from (8.10)-(8.11), 


Xi{s,Q^,i) ^ Xi(s, 0,t) = 0, 


and 
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(8.18) 


TT dt 

where we have introduced the differential cross section for the real photoproduction pro¬ 
cess ey — uW: 
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1 + h/s^ 

with u = {I — k'Y = — s — t. Equation (8.19) agrees with the analytical result already 

presented in [98,99], obtained using the helicity amplitude technique. Using (8.17) and 
(8.18), we can write 
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( 8 . 20 ) 




120 


uW Production in ep —^ uWX 


where x = s/s and l/^ex^x^xs) is the inelastic component of the equivalent photon distri¬ 
bution of the proton given in (5.44). As already pointed out (see also [98]) there is some 
ambiguity in the choice of the scale of here taken to be s = xs, which is typical 
of all leading logarithmic approximations, and any other quantity of the same order of 
magnitude of s, like —t or —h, would be equally acceptable within the limits of the EPA. 
The numerical effects related to the scale dependence of the inelastic photon distribution 
are discussed in the next section. 


The cross section relative to the elastic channel, ep uWp, has been calculated in 
[24]; it can be written as (5.21), namely 
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( 8 . 21 ) 


with t = —integrated over the range already dehned by (8.15), and Smin given by 
(8.13). The limits of integration of t are the same as in (8.16) and the invariant Hi{t) is 
given in (5.19) in terms of GE{t) and GM{t), the electric and magnetic form factors of the 
proton, respectively. 


Again, in the limit s 3> and s 3> —t, the cross section factorizes and is given by 


aei(s) ~ = 
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( 8 . 22 ) 


where x = s/s and is the universal, scale independent, elastic component of the photon 
distribution of the proton introduced in (3.2). 


8.2 Numerical Results 

In this section, we present a numerical estimate of the cross sections for the reactions 
ep uWX and ep uWp, calculated both exactly and in the EPA, in the range 
100 < ^/s < 2000 GeV. We take Mw = 80.42 GeV for the mass of the W boson and 
Gp = 1.1664 X 10“^ GeV“^ for the Fermi coupling constant [104]. All the integrations are 
performed numerically. In the evaluation of (8.5) and (5.44) we assume the LO Gallan- 
Gross relation (5.46) and we use the ALLM97 parametrization of the proton structure 
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function ^ 2 ( 0 ;, Q"^) [70], which is expected to hold over the range oi xb and studied so 
far, namely 3 x 10“® < xb < 0.85 and D < < 5000 GeV^. As for the QED Compton 

scattering process discussed in the preceding chapters, we do not consider the resonance 
contribution separately but we extend the ALLM97 parametrization from the continuous 
(pp 2 ^ 3 GeV^) down to the resonance domain ((m^ + m)^ < < 3 GeV^): in this 

way it is possible to agree with the experimental data averaged over each resonance. In 
our analysis, the average value oi xb always lies within the kinematical region mentioned 
above, where the experimental data are available. On the contrary, the avarage value of 

becomes larger than 5000 GeV^ when > 1200 GeV, so we need to extrapolate 
the ALLM97 parametrization beyond the region where the data have been htted. Our 
conclusions do not change if we utilize a parametrization of F 2 {xb,Q‘^) whose behaviour 
at large is constrained by the Altarelli-Parisi evolution equations, like GRV98 [ 68 ]. 
The electric and magnetic form factors, necessary for the determination of the elastic 
cross sections in ( 8 . 21 ) and ( 8 . 22 ), are empirically parametrized as dipoles, according to 
( 5 . 51 ). 

At the HERA collider, where the electron and the proton beams have energy Eg = 27.5 
GeV and Ep = 920 GeV respectively, the cross section is dominated by the inelastic 
channel: Uei = 2.47 x 10“^ pb, while Uinei = 3.22 x 10“^ pb; therefore the expected 
integrated luminosity of 200 pb“^ would yield a total of about 11 events/year. 

Figure 8.2 shows a comparison of the inelastic cross section calculated in the EPA, 
with the exact one, Uinei, as a function of y/s, where several scales for are 

proposed, namely = E ~u, —t in (5.44). It turns out that the choice of —t does not 
provide an adequate description of Uinei, while s and —it are approximatively equivalent 
in reproducing Uinei. In particular, the choice of —it is slightly better in the range 300 < 
y/s < 1000 GeV, while s guarantees a more accurate description of the exact cross section 
for y/s > 1000 GeV. At HERA energies, = 3.64 x 10“^ pb, 3.51 x 10“^ pb and 
3.07 X 10 ^ pb for Qmax = E and —t, respectively. In the following we will £x the 
scale to be s, in analogy to our previous studies about the QED Gompton scattering 
process in ep —> eyX [3-5]. 

In [4, 5] it was suggested that the experimental selection of only those events for which 
s > restricts the kinematics of the process to the region of validity of the EPA and 
improves the extraction the equivalent photon distribution from the exact cross section. 
The effect of such a cut on the reaction ep uWX is shown in Figure 8.3 and the 
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Figure 8.2: Exact and approximated (EPA) inelastic cross sections of the process ep 
uWX as functions of ^/s. The different scales utilized in the calculation of the approxi¬ 
mated cross section are written in the brackets. 



4s (GeV) 


Figure 8.3: Exact and approximated (EPA) inelastic cross sections of the process ep 
nWX as functions of ^/s. The scale s is utilized in the calculation of the approximated 
cross section and the kinematical cut s > is imposed in the exact one. 










8.2 Numerical Results 


123 
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Figure 8.4: Exact and approximated (ERA) total ( = elastic + inelastic) cross sections 
of the process ep uWX as functions of ^/s. The exact elastic component, which is 
indistinguishable from the approximated one, is shown separately. 


reduction of the discrepancy is evident at large ^/s, but not at HERA energies, where Uei 
and (Jinei are unchanged. 

In Figure 8.4 the total (elastic + inelastic) exact cross section is depicted as a function 
of ^/s, together with the approximated one. Here the kinematical constraint s > is 
not imposed on the exact cross section. The average discrepancy is reduced to be about 
2%, due to the inclusion of the elastic channel, better described by the EPA (average 
discrepancy 0.05%). The elastic component is also shown separately, and it agrees with 
the curve presented in Figure 3 of [24]. For = 318 GeV, = 2.47 x 10“^ pb, in 
perfect agreement with the exact value Uei. 

We compare now our results with the ones already published. In [102], taking into 
account the photon exchange only (Figure 8.1) and with no further approximation, hxing 
Mw = 83.0 GeV, sin^ 9w = 0.217, Ef. = 30 GeV, Ep = 820 GeV and using the parton 
distributions [105] (Set 1), together with the cuts > 4 GeV^ and IF^ > 10 GeV^, the 
value (Tinei = 3.0 X 10“^ pb was obtained. This is in contrast to Uinei = 1-5 x 10“^ pb, 
calculated using (8.5) with the same sets of cuts, values of the energies, Mw and parton 
distributions utilized in [102]. The authors of [102] also report the value Uinei = 4.0 x 10“^ 
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pb, obtained in [101] with a similar analysis at the same energies, using = 78 GeV 
and sin^ = 0.217. The lower limit on was taken to be 0(1) GeV^, but not explicitly 
mentioned. Even with the ALLM97 parametrization, which allows us to use no cutoff on 
Q^, we get (Jinei = 3.1 x 10“^ pb, far below 4.0 x 10“^ pb. No analytical expression of the 
cross section is provided in [101,102], which makes it difficult to understand the source 
of the discrepancies. 

Finally, an estimate of the the ep uWX cross section is also given in [98,99], 
utilizing an inelastic equivalent photon distribution slightly different from the one in 
(5.44), which can be written in the form 

. (8.23) 


where 
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cut 
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Qmax — and Qcut = 1 GeV^. Equation (8.23) can be obtained from (5.44) 

neglecting the mass term and approximating the integration over Q^. In the calculation 
performed in [98,99], Qmax) is convoluted with the differential cross section for 

the real photoproduction process in (8.19). At y/s = 300 GeV, hxing Mw = 84 GeV, 
sin^ = 0.217 and using the parton distribution parametrization [105] (Set 1), we get 
(Tinei = 1-6 X 10“^ pb, very close to the value 1.5 x 10“^ pb published in [98,99]. 


8.3 Summary 

We have calculated the cross section for the inelastic process ep —> uWX, both exactly 
and using the equivalent photon approximation of the proton, in order to test its accuracy 
in the inelastic channnel and complete the study initiated in [24], limited to the elastic 
process ep uWp. The relative error of the approximated result with respect to the 
exact one is scale dependent; hxing the scale to be s, it decreases from about 10% at 
HERA energies down to 0.5% for y/s = 1500 GeV, then it slightly increases up to 3% 
for y/s = 2000 GeV. In conclusion, even if not so remarkable as for the elastic channel, 
in which the deviation is always below one percent [24], the approximation can be con¬ 
sidered quite satisfactory. We have compared our calculations with previous ones in the 
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literature and found that they are in agreement with [98,99], but disagree with [101,102], 
Furthermore, we have estimated the total number of vW events expected at the HERA 
collider, including the elastic and quasi-elastic channels of the reaction. The production 
rate turns out to be quite small, about 11 events/year, assuming a luminosity of 200 pb“\ 
but the process could still be detected. 





Chapter 9 

Summary and Conclusions 


In a photoproduction reaction initiated by a nucleon N (= p, n), the nucleon can be 
considered to be equivalent to a beam of photons, whose distribution can be computed 
theoretically. It was shown, see (3.21) and (3.57), that the cross section of the reaction 
factorizes, being given by the convolution of the universal (process independent) photon 
distribution of the nucleon, with the corresponding photon-nucleon cross section. 

The advantage of this model is that the nucleon’s photon content can be utilized to 
simplify the calculations of photon-induced subprocesses in elastic and deep inelastic eN 
and NN collisions, commonly described in terms of the electromagnetic form factors, 
structure functions and parton distributions (A)/(a;,/i^) (with f = q, q, g) introduced in 
Chapter 2. 

The polarized and unpolarized photon distributions (A) 7 (a;, p^) were evaluated in the 
equivalent photon approximation (EPA) in Chapter 3. Both of them consist of two com¬ 
ponents: the elastic ones, (A) 7 ei, due to A — 7 A, and the inelastic ones, (A) 7 inei, due to 
N — 7 A, with X ^ N. The elastic photon components turned out to be scale indepen¬ 
dent and uniquely determined by the electromagnetic form factors. The inelastic photon 
components were radiatively generated in a model where they vanish at a low resolution 
scale /iQ. At large scales the resulting distributions are rather insensitive to details at 
the input /Tq, therefore such not compelling vanishing boundary conditions are supposed 
to yield reasonable results for (A) 7 inei(a;,/x^). However at the low scales characteristic of 
hxed target experiments, like HERMES at DESY, (A) 7 inei(a;,/i^) depend obviously on the 
assumed details at /Xq. We did not investigate any effect due to 7 inei(a^, ho) 7^ should 

rather be examined experimentally if our expectations based on the vanishing boundaries 
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turn out to be in disagreement with observations. 

Chapter 4 was devoted to demonstrate the feasibility of the measurement of 
(A) 7 (a;,/i^). We studied muon pair production in electron-nucleon collisions eN —> 
via the subprocess 7 ^ 7 ^ —> and the QED Compton process eN —e 7 X 

via the subprocess > ej for the HERA collider and the polarized and unpolarized 

hxed target HERMES experiments. The muon pair production process was evaluated in 
the leading order equivalent photon approximation, hence we considered just the simple 
2 —2 subprocess 7 ^ 7 -^ —, instead of the full 2 —3 or (even more involved) 2 —4 
subprocesses 7 ®g —and eg —Similarly, the analysis of the deep inelastic 
QED Compton process was reduced to the study of the ey scattering instead of 

eg —eyg. We concluded that the production rates of lepton-photon and dimuon pairs 
would be sufficient to facilitate the extraction of the polarized and unpolarized photon 
distributions in the available kinematical regions. 

The full 2 —3 QED Compton process in ep —> eyp and ep —> eyX, with unpolarized 
incoming electron and proton, was calculated in a manifestly covariant way in Chapter 
5 by employing an appropriate parametrization of the structure function F 2 [70] and 
assuming the Callan-Gross relation (2.85) for F^. These “exact” results were compared 
with the aforementioned ones based on the EPA, as well as with the experimental data 
and theoretical estimates for the HERA collider given in [31]. Although the cross section 
in the elastic channel is accurately described by the EPA (within 1%), this is not the case 
in the inelastic channel. It turned out that the agreement with the EPA is slightly better 
in bins, where x.y, defined in (5.28), is the fraction of the longitudinal momentum of 
the proton carried by the virtual photon, compared to the bins in the leptonic variable 
xi (5.31), used in [31]. In addition the results obtained by an iterative approximation 
procedure [34, 82], commonly used [31], were found to deviate from our analysis in certain 
kinematical regions. 

The virtual Compton scattering in unpolarized ep —> eyp and ep ^'jN, where 
the photon is emitted from the hadronic vertex, was studied Chapter 6 . It represents the 
major background to the QED Compton scattering. New kinematical cuts were suggested 
in order to suppress the virtual Compton background and improve the extraction of the 
equivalent photon content of the proton at the HERA collider. The total (elastic - 1 - 
inelastic) discrepancy of the exact cross section with the approximate one was reduced 
to be about 2%, which should be compared with the value 14% relative to the cuts [31] 
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discussed in Chapter 5. 

The 2 —3 QED Compton process in longitudinally polarized lepton-proton scattering 
was analysed in Chapter 7, both in its elastic and inelastic channels. For our numerical 
estimates we utilized the BKZ parametrization [71] of the spin dependent structure func¬ 
tion Qi and we neglected g 2 . It has to be noted that the effects of g 2 , as well as of Fl in the 
unpolarized cross section, may become relevant at low scales and it would be interesting 
to take them into account in future studies. The kinematic cuts necessary to extract the 
polarized photon content of the proton and to suppress the major background coming 
from virtual Compton scattering were provided for HERMES, COMPASS and the future 
eRHIC experiments. For the HERMES and eRHIC kinematics, we found that the total 
cross section can be described by the EPA with an error estimated to be less than 10%. 

The reliability of the EPA in reproducing the cross section of the process ep vWX 
was investigated in Chapter 8. In order to examine this issue, both the subprocesses 
eyP —> vW and eg —vWq were studied. The relative error of the approximate result 
with respect to the exact one was found to be less than 10% in the inelastic channel and 
less than 1% in the elastic one, this last value being in agreement with [24]. 

To conclude, we suggested how to access the photon distributions (A) 7 (a;,/i^) in eN 
collisions, studying the accuracy of the equivalent photon approximation of the nucleon 
in describing different processes and trying to identify the kinematical regions of its va¬ 
lidity. Our Endings should now be confronted with experiments. As already mentioned, 
these measurements would not only be interesting on their own, but would provide the 
opportunity of getting additional and independent informations concerning the structure 
functions Fi ^2 and gi^ 2 , underlying the (inelastic) photon distributions. In particular, 
being the kinematics of QED Compton events difierent from the one of inclusive deep 
inelastic scattering, due to the radiated photon in the final state, it provides a novel way 
to access the structure functions in a kinematical region not well covered by inclusive 
measurements [2, 5, 26,31-33]. Hence it represents a valuable complementary tool to have 
a complete understanding of the structure of nucleons. 





Appendix A 

Notations and Conventions 


Units 

Natural units h = c = 1, with h = /i/27r, are used throughout this thesis, where h and 
c denote the Planck constant and the speed of light, respectively. 

Relativistic conventions 

The metric tensor ga/s = g°‘^, with a, /? = 0,1, 2, 3, is given by 

^00 = + 1 , ^ ^ otherwise = 0. (A.l) 

The contravariant vectors of the space-time coordinate and energy-momentum of a particle 
of mass m and energy E are given by 

= p“ = (U,p), (A.2) 

where 

E = \/+ m? . (A.3) 

The covariant vectors are 

Xa = gafiX^ = (A.4) 

Pa = gagp^ = {E,-p), (A.5) 

and hence their scalar product is dehned by 


X ■ p = Xap'^ = gagx^^p^ = tE — xp, 


(A,6) 



132 


Notations and Conventions 


where it is understood that repeated indices are summed. Furthermore the totally anti¬ 
symmetric tensor = —Sappa is dehned so that 


^0123 — +1- 


(A.7) 


7 matrices 

The Dirac 7 matrices 7 " = ( 7 °, 7 *), with i = 1,2,3, satisfy the anticommutation 
relation 




The matrix 7 ^ is dehned as 


7^ = 75 = i7°7S^7^- 


The hermitian conjugate of 7 " is taken to be 

= 7 ° 7 " 7 °, 

so that according to the dehnition (A.9) one has 

75t = 

The scalar product of the 7 matrices and any four-vector A is dehned as 


4 = = 7°A° - 7^A^ - 7^712 - 


(A.8) 

(A.9) 

(A.IO) 

(A.ll) 

(AA2) 


Spinors and normalizations 

The normalization of the one particle state \p, a) with momentum p and other quantum 
numbers a is taken to be 

{p, a I pa) = {27r)^2E6^{p — p)6aa'- (A.13) 

The Dirac spinors for fermion, u{p,s), and antifermion, v{p,s), with momentum p, spin 
s and mass m satisfy 


0 = ij) — m)u{p, s) = u{p, s){j[) — m) 
= {jj) + m)v{p, s) = v{p), s){j[) + m), 


(A.14) 

(A.15) 
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where u{p, s) = u^{p, 5 ) 7 °. The Dirac spinors are normalized in such a way that 

u{p, s)u{p, s') = 2m6ss', v{p,s)v{p,s') = -2mSss', (A.16) 

which hold also for massless particles and antiparticles. Accordingly the completeness 
relations read 

^u{p,s)u{p,s) = i^ + m), (A.17) 

S 

'^^v{p, s)v{p, s) = — m). (A.18) 


Spin projection operators 

The polarization vector s of a relativistic spin ^ particle or antiparticle with momentum 
p is a pseudovector which fulhls 

5^ = -!, s-p = 0. (A.19) 

The projection operators onto a state with polarization ±s are known to be 

V{±s) = ^ (1 =*= 3^^) for particles, (A.20) 

V{±s) = ^(1=F7^;^) for antiparticles. (A.21) 

For massless particles and antipartcles only longitudinal polarization is possible and the 
operators introduced above become helicity projectors: 

V± = ^ (f 7^) fo^ particles, (A. 22 ) 

V± = ^ (f T 7 ^) for antiparticles, (A.23) 

where the subscript + indicates positive helicity, that is spin parallel to the momentum, 
and the subscript — indicates negative helicity, that is spin antiparallel to the momentum. 





Appendix B 

Photon-induced Cross Sections 


The well-known unpolarized and polarized cross sections for dimuon production in 77 — 
and for the Compton scattering reaction ej 67 are recalculated in a manifestly 
covariant way in Appendices B.l and B.2 respectively. 


B.l Dimuon Production 


The process of photon-photon annihilation into a muon-antimuon pair 

7 (^ 1 )+ 7 (^ 2 ) ^ +AC(/2) (B.l) 

is described to lowest order in QED by the two Feynman diagrams in Figure B.l. 

We dehne the Mandelstam variables 


S=(A^l + *:2)^ 


(B.2) 


which satisfy the constraint 

s -|- t -|- h = 0, 

as the masses of the particles are neglected. The scattering amplitude reads 


(B.3) 


M = -eht(ki)e^/(k.,)u(l2) 


[ (h - h)' 




(ki-hf 


7"(fe- (' 1 ) 7 '’ 




(B.4) 


where e is the proton charge and e^*{ki), with i = 1, 2, denotes the polarization vector of 
an incoming photon with momentum ki and helicity A*. The Dirac spinors, whose spin 
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Photon-induced Cross Sections 



Figure B.l: Feynman diagrams for 77 — 


dependence is not shown explicitely here, were introduced in (A. 14) and (A. 15). 


Unpolarized cross section 

Making use of the completeness relations (A.17), (A.18) and 

^e^(/c)ej(/c)* = 

A 

one can perform the spin sums and get 

+ \M+-? 


(B.5) 




spins 


= e 


+ 


{h-hY {k^-hY, 




ki ■ k2 — h ■ k2 — h ■ ki . 

4 Tr [ l/ 2 l/l\ 


{k^-liY{k2-liY 


(B. 6 ) 


with the subscripts ± referring to the two possible values ±1 of the helicity A* of each 
photon. The traces can be easily calculated and expressed in terms of the variables (B.2), 
so that (B. 6 ) reduces to the simple form 


|Af|2 = 2e'‘(l + ?), (B.7) 

the last (interference) term in (B. 6 ) being zero due to (B.3). From the dehnition (4.19), 
that is 


da 

dt 


1 

IbTTS^ 



(B.8) 
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one gets the final result for the unpolarized cross section: 


da 2710 ^ / i it 
di s‘^ [it i 


(B.9) 


Polarized cross section 

It can be shown (see discussion below (B.13)) that the two photons annihilate into a 
muon-antimuon pair only if they have opposite helicities, that is 


\M++f = 0, 


(B.IO) 


which, substituted in the first line of (B.6), together with (B.7), gives the following result 


lAAfp = 1 | 

needed for the calculation of the polarized cross section, 


(B.ll) 


di IbTTS^ 

Substituting (B.ll) into (B.12), one obtains 


|AMp 


(B.12) 


dAd dd 27rQ;^ /1 it 


\u t 


(B.13) 


To demonstrate the validity of (B.IO), one has to inspect the helicity structure of the am¬ 
plitude (B.4). Making use of the helicity projection operators (A.22), (A.23), one can perform 
the following substitutions in (B.4) 


u±{h) = -n(/2)(l =F7^), 

v±{h) = ^ (1 =F7^)^^(h), 


(B.14) 


where the subscripts on the spinors denote the helicities of the corresponding fermions. In this 
way one notices that, when the helicities of the fermions are the same, the amplitude vanishes 
because of a mismatch in the projection operators (1 ± 7^)/2. Hence 


M++ = M++; + _ + M+ + ;_ + 


(B.15) 


and, in particular. 


a 7+;+- = 7fe)(fe- m*(h) 


+ {k 2 -llY ^ {k 2 ) v{li). 


(B.16) 
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Photon-induced Cross Sections 


It is possible to introduce a representation for the polarization states of the photons without 
choosing any specific frame. This means that the calculation of the helicity amplitudes in 
(B.15) can be performed in a covariant way. The circularly polarized states of a photon with 
four-momentum k can be expressed as 

e± = T^(ell±ze;^), (B.17) 

where the vectors and e"*" satisfy the conditions 

(ell)" = (e^)" = -l, (B.18) 

A;, ell = fc.e^ = ell-e^ = 0. (B.19) 

If we introduce two arbitrary vectors, p and q, for convenience taken to be light-like {p^ = ?" = 0), 
then we can write 

ejl = 2y/2N[{q-k)pi,-{p-k)q^], (B.20) 

(B.21) 

where the normalization factor, which is the same for both polarization vectors, is fixed by 
(B.18) to be 

N = \[{P-q){p-k){q-k)]~^. (B.22) 

In the amplitude (B.4) the photon polarization vectors only appear in the combination using 
the identity 

= {laldl'i - lagp-t + 7/35a7 “ , (B.23) 

one can write [78]: 

±iV[^|i^(l ± 7 ^) -|5^^(1 T7^) T 2(p- g)|^ 7 ^]. (B.24) 


The choice of the four-vectors p and q is arbitrary. When the photon line is next to an external 
fermion or antifermion line, one usually exploits this freedom of choice and takes either p or 
q equal to the fermion momentum. Use of the equations (A.14), (A.15) for massless fermions 
leads to great simplifications. For the process under study we fix p = I 2 and q = li, so that 


^+{ki) = N +1^) -hhhi^ -1^) - 2 (/i • l2Wii% (b.25) 

with i = 1, 2 and N = (2s tft)” 2 . The factor (1 — 7 ^) in (B.16) only selects the first and the 
third term in the RHS of (B.25), but the first one does not contribute because liv{li) = 0. 
Therefore 


M+ + ; + - 


-^e^iV"s^tt(/2)(l 

-^e^iV"s^tt(/2)(l 

0 , 


5 r ^2(^1- mi 

(A:i-/i)2 


+ 


UH 2 - 

{k2-h? _ 


v{ll) 


l^)i}2+h)v{h) 


(B.26) 


where the last equality follows from u(/ 2)/2 = 0 and liv{li) = 0 . Analogously one can show that 
M++\ _h = 0) which, together with (B.26), gives (B.IO). 
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B.2 Compton Scattering 


The Compton scattering process 

e(/) + 7 (/c) ^ e(/') + 7 (/c') (B.27) 

can be described in terms of the Mandelstam variables 


s = {l + k)\ t={l-l')\ u = {k-l')\ 


(B.28) 


which satisfy (B.3), since the photons are taken to be on shell and the electron mass is 
neglected. The corresponding amplitude to lowest order in QED reads, see Figure B.2, 


M = e^eaC*f^u{l') 




(/ + kf 


(/ - cy 


r 


2(1 + 7^)«(0 


(B.29) 


where e, e* are the polarization vectors of the incoming and outgoing photons; u, u are 
the Dirac spinors of the initial and hnal electrons respectively. We have assumed that 
the incoming electron has right-handed helicity, i.e. its spin is parallel to the direction of 
motion. 



Y Y 



Figure B.2: Feynman diagrams for ey —>• ey. 

The leptonic tensor is given by 

= Yi M'JU = T% + ry, (B.30) 

spins 

and being its symmetric and antisymmetric parts. Using the completeness rela¬ 
tions (A.17), (B.5) one hnds 

T^p{l]l',k') = e^Tr|^^^-^^[(l+7^)f7„^'7^]--^^-^^[(l + 7^)^'7„/'7^] 
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+ (f +mi - k'r ~ 

+ (l+7'')/'Ts(W'-W'+fnTan}. (B.31) 

which is valid also for an incoming virtual photon. 


Unpolarized cross section 

In terms of the variables (B.28), the symmetric part of the leptonic tensor reads, for 
real photons, 




46 

7 : 9ap + 2.S lalf3 + 2 -U I' I'g 

su \ 2 ^ 

+ i {Ictl'fS + h^a) ~ ^ {I'ak'p + l/ska) 

+ u + I'^k'^) 


(B.32) 


From (B.32), one can determine the amplitude squared and averaged over the spins of 
the incoming photon 

spins 

= -2e^(^| + |^, (B.33) 

hence the cross section (B.8) for the process ey —> ey is given by 


dd 

dt 



(B.34) 


Polarized cross section 

From the antisymmetric part of the leptonic tensor, 

Ai 

T^fs{k,l',k') = —^e^^^^{slP + uknk% (B.35) 

su 

with k = I' + k' — I, one gets 

|AM|2 = = 



(B.36) 
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where is the antisymmetric part of the photon polarization density matrix (3.54), 
and the cross section (B.12) for the process ey —ey, where the incoming electron is 
longitudinally polarized and the incoming photon circularly polarized, reads 


dAd 

di 


27ra‘^ 



(B.37) 





Appendix C 

Kinematics of the QED Compton 
Scattering Process 


The kinematics of the QED Compton scattering process in ep e'yp and ep —> e 7 X is 
described in the center-of-mass frame of the outgoing e — 7 system and in the laboratory 
frame, both for a collider and a fixed-target experiment. 


C.l Elastic Channel 

In the elastic channel the process reads 

e{l)+p{P)^e{n+^{k')+p{P'), (C.l) 

with P^ = P'^ = = 0, and f ~ 0 . 

Electron-photon center-of-mass frame 

In the e — 7 center-of-mass frame, we choose the axis to be along the direction of 
the virtual photon exchanged in the reaction, see Figure 5.1. The four-momenta of the 
initial electron and proton are given by: 

/ = {El 0,0, -El), P = {e;, p;sin0;, o, p;cos0;), (c.2) 

where P* = ^— m?. For the outgoing electron and photon we have 

/' = P'*(l, sin6** COS99*, sin6'* sinv?*, cos6'*), (C. 3 ) 

k' = P'*(l, — sin 6^* cos v?*, — sin6'* sin 99*, — cos 6^*). (C. 4 ) 

The four-momentum of the virtual photon is 


k = i' + k' -i = {EiQ, 0 , p :). 


(C.5) 
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with = t and 

El = ^Ef + t . (C.6) 

The overall momentum conservation allows us to write the four-momentum of the hnal 
proton as 


P' = l + P-I' -k'. (C.7) 

We introduce the following Lorentz invariants 

s = (/' + k'f = (C.8) 

i={l- I'f = -2E*E'*{1 + cosr), (C.9) 

u = {l- k'f = -2E*E'*{1 - cosr), (C.IO) 

s = (/ + P)2 = + 2E;e: + 2EIP; cos Of (C.ll) 


T = {P- I'f = -2E'*{E;-p; sine* sine; cos if* - p; cose* cos e*p). (C.12) 

u = {p - k'f = m'^ - 2 e'*{e; + p;sine*sine;cos (p* + p;cose*cose;), (c.13) 

In addition they satisfy: 


s + i + u = t, s + T + U = —t + 3m^. (C.14) 

Using the relations above, it is possible to write the energies of the particles in the labo¬ 
ratory frame in terms of the integration variables s,t,t and the constant s\ 


e: = 


s — t 


El = 


s + t 


and 


s — rrP + t a/(s — 


77'* _ 

2VI ’ 


2 v^ 


Similarly, for the angles we have 



(C.15) 


(C.16) 


* 


t — s — 2i 


cose 


s — t 


(C.17) 
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cos 6* 


2s {s — w?) — {s — t){s — w? + t) 
(s — t)[(s + t — w?Y — Asw?]^ 


(C.18) 


In particular (C.12) and (C.13), through (C.16)-(C.18), express T and U in terms of our 
integration variables s, t, t, (p* and we have used them to relate the laboratory frame 
variables to the integration ones, as shown in the following. 


Laboratory frame 

We choose the laboratory frame such that the initial proton moves along the axis 
and the outgoing electron has zero azimuthal angle. The four-momenta of the particles 
are given by: 


/ = {Ee, 0, 0, -Ee), P = {Ep, 0, 0, Pp), (C.19) 

where Pp = \JEp^ — np, and 

I' = E'^{1, sin6^6, 0, cos6^6), 

k' = sin6*.y COS0..),, sin6*.y sin^..),, cos6*^). (C.20) 

The Lorentz invariants are: 

s = {I' + k'Y = 2E'^E'^{1 — sin 6*e sin 6^.^ cos 0..), — cos 6^6 cos 6^^), (C.21) 


i={l- I'f = -2EeE'^{l + COS^e), 


(C.22) 


u = {I — k'Y = —2EeE'^{l + cos0j), 
s = (I P P)^ = TTi^ P 2Eg{^Ep P Pp), 

T = {P- I'f = m^- 2E'fEp - PpCosOf. 


(C.23) 

(C.24) 

(C.25) 


U = {P- k')^ =m^ - 2E'fEp - PpCosOf, 


(C.26) 


The polar angles in the laboratory frame can be written in terms of the invariants and 
the incident energies: 


Epi- Ee (T - m^) 

cos Oe = - , 

PptpEeiT- m?) 


(C.27) 


Ep it — s — i) — Ee {U — w?) 
Pp {t — s — i) P Ee {U — m?) 


(C.28) 
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In the same way, for the energies of the hnal electron and photon we have: 

^ i Pp + Ee (T - m^) 


(C.29) 




Pp {s — t + i) — Ee {U — m?) 


s — 


The azimuthal angle of the outgoing photon is: 

2E'^E' (1 — cos6e cos 6*^) — s 


COS(p^ = 


2E'^E'^ sin 9e sin 6^ 


(C.30) 


(C.31) 


which is related to the acoplanarity angle by 0 = [tt — 0^ |. Using (C.12) and (C.13) 
together with (C.16)-(C.18), the formulae above for cos6^6, cos6*^, E'^ and cos0^ can 

be expressed in terms of the invariants and the incident energies. Equations (C.27)-(C.31) 
are needed to implement numerically the kinematical region under study, since they relate 
the laboratory variables (energies and angles) to the ones used for the integration. 

For a fixed-target experiment, (C.19) should be replaced by 


I = (Ee, 0, 0, Ee), P = (m, 0, 0, 0), 


(C.32) 


assuming that now the axis is along the incoming electron direction. Keeping (C.20) un¬ 
changed, (C.21)-(C.31) are still valid, provided one performes the following replacements 


Ep — m, 


n — 9p 


(C.33) 


C.2 Inelastic Channel 


For the inelastic channel of the process, 

e{l)+p{P)^e{n+^{k')+X{Px), 

most of the expressions remain the same as given in Appendix C.l. 


(C.34) 


Electron-photon center-of-mass frame 

The relations among the invariants are now: 


s P i u — —, 


s + T + U = 3m^ -I- 


(C.35) 


where = —t and xb is given in (5.38) with = {P + l — I' — k'Y- The only formulae 
which are different from the elastic channel are the ones involving E* and cos 9*. So 
(C.16) will be replaced by 


e; = 


s-Q^-W^ 


P* = 


'{s - Q 2 _ _ 45^2 

2^1 


E'* = ^ , (C.36) 




C.2 Inelastic Channel 


147 


while for the angles: 


and 


cos^* 


+ s + 2t 
s + g2 


2s (s - m2) - (s - g 2 _ ^ g 2 ) 

(s + g2)[(s - g2 - - 4sm2]5 


(C.37) 


(C.38) 


Equations (C.36)-(C.38) reduce to (C.16)-(C.18) of the elastic channel for IE = m and 

g' = -t. 


Laboratory frame 

The definition of the invariants for the inelastic channel is the same as for the elastic 
one: (C.21)-(C.26); (C.35) describes the relation among them. The expressions of cos 6^6, 
cos 6*.y, Eg, E'^ and cos in terms of the integration variables lE^, s, g^, t, (p* are given by 
(C.27)-(C.31) together with (C.13), (C.12) as before, but now (C.36)-(C.38) will replace 
(C.16)-(C.18). 





Appendix D 

Matrix Elements for the Unpolarized 
QEDCS and VCS Processes 


In this Appendix we give the explicit expressions of the matrix elements relative to QED 
Compton scattering (QEDCS) and virtual Compton scattering (VCS) in ep e'jp and 
ep G'yX, with incoming unpolarized electron and proton. 


D.l Elastic Channel 

The amplitude squared of the process 

e(/) + p(P) ^ e(/') + 7 (A;') + p(P') (D.l) 

can be written as in (6.6), i.e. 

I Mel 1^ = I f + I |2 _ 2 gfJe (D.2) 

We will make use of the Lorentz invariants (4.39), (5.2), 
s = (/ + P)^ t = -Q^ = k\ 

s = {l + kf, i = (/-/')^ u = {l-k'f, (D.3) 

with k = I' + k' — I, and 

S={P' + k'f, U ={P- k'f, (D.4) 

which satisfy the relation S = —{s-\-u + U — 2m?). According to the approximation 
explained in Section 6.1, we have 

I l' = ^ P')T“^(/; k'), (D.5) 
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Taf 3 being given by (5.8), and 

£a/3(p. p/) ^ 2 + P'"P^ - ^"^(P ■ P' - m^)}, (D.6) 

with Pi denoting the Dirac form factor of the proton. Furthermore, 

lA/ys = (D.7) 

where L^/s is the leptonic tensor of the non-radiative deep inelastic scattering ep eX 
introduced in (2.17) and 


V^p{P-,P\k') = 

- Tr [(f + m)7^(f - ^'+ m)7„(f’'+ m) 7 / 3 (f - ^'+ m)7^] 

- Tr[(f + m)7/,(f- ^'+ m)7«(f’'+ m)7^(f"+ ^'+ 771)7^3] 

- Tr[(f + m)7“(f'+ ^'+ m)7^(f'+ m)7^(f - ^'+ 777)7'^] |. 

(D.8) 


Explicit calculation gives 


M, 


QEDCS 12 


el 


4e® 
t s h 


71 + 


2m^ 

- ID 

t 


F^t), 


(D.9) 


Mjcs 


4e® r ^ 2m‘^ ^ 

i U' S' ~ i U' S' 


FKi\ 


where 


71 = 2t^ -2 t{s-2 s' -U') + s^ -2 s s' 

+ 4 s'2 + 2 s' 77 + + 4 s' P' + 2 77 P' + 2 P'^, 

B = 2 — 2 t (s + 77 ) + s^ + 77 ^, 

P = (s + 77 )^ [t^ + s^ — 2 t (s — s') — 2 s s' 

+2 s'^ + 2 s' 77 + 77 ^ — 2 777^ (s + 77 — t)] 

+2 (s + 77 ) [P — t s + 77 (—s + 2 s' + 77 )] P' 
+2 [t^ + 77 ^ — t (s + 77 )] P'^. 


(D.IO) 


(D.ll) 

(D.12) 


(D.13) 


In (D.8)-(D.13) we have used the notations s' = s — m^. S' = S — ni^, U' = U — m? for 
compactness and the electron mass has been neglected. 
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The interference term between QEDCS and VCS is given by 


2 _ i; ((; r, k')Mf(P-, P', k'), 


2tt 


(D.14) 


with 


A;') = ^ Tr[^7"(^ + ^ ni/^a(/ - [ (D.15) 


and 


M%{P-,P',k') = e^F,{t)F,{t)\-TT[{r+mh^{r+m)Y{r+¥+mhp] 


+ — Tr[ (f"+ m) 7 „ (f + m) 7 ^(f - ^'+ m)7;3 ] 


The hnal result can be written as 


2SeAC°“A7“- = 4e« F.(«)F.(«), 


t s u t U' S' 


(D.16) 


(D.17) 


with 


D = {(s + h) [t h + s' (s + h)] + [s (s + h) — t (s — h)] U'} 

[2 C + — 2 s s' + 4: s'"^ + 2 s' ii + iS — 2 t {s — 2 s' — U') 

+ 4s' U' + 2u U' + 2U'\ 


(D.18) 


E = -s' -s^ {s' -2 u + U') - u {7 s'+ 2 U') - s {7 s'+ 2 u + 5 U') 

+ 2 [s (h — U') + ii {u + U')] — t {s + u) [s (—2 s' + 3 h — 3 U') 

+ u {—2 s'+ u+ U')]. (D.19) 


D.2 Inelastic Channel 

The amplitude squared of the electron-quark scattering process (6.9), 

e(/) + q{p) e(/') + 7 (A;') + g(p'), (D.20) 

which is needed, according to (6.11), for the calculation of the inelastic reaction 

e{l) + p{P) ^ e(/') + 7(A;') + X{Px) (D.21) 


is given by 


M<? 1^ = I M<?QEDCS |2 ^ I Mqwcs |2 _ 2gfJeM9QEDCS]vf<?vcs7 ^ 22) 
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Within the framework of the parton model, we assume that the initial quark is collinear 
with the parent proton, p = xbP- If we define 

S={p' + k')\ U={p- k')\ (D.23) 

then U = xbU', S = —{s + u + xb U'), with U' = U — w?] U and the other invariants 
useful to describe the process being given in (D.3) and (D.4). Neglecting the masses of 
all the interacting particles in (D.20), the QED Compton scattering amplitude squared 
can be written as 

I M,QED0S i" ^ 1 (D.24) 

where Tai 3 is given in (5.8) and Wafs is the quark tensor, 

WafiiP] p') = 2 Cg (pV^ + P'"/ - ■ P'), (D.25) 


with Cq being the charge of the quark in units of the proton charge e. Analogously, for 
the VCS amplitude squared one has 

I A/.VOS i" = (D,26) 

where 


4:6*6* 
SU 


TlfApip'X) = 


2 fl'o/9 + 2tt) + 2S PaPiS + 2f/ PaPp 

+ {i + t){paP'p + Ppp'a) -{s -i) {p^k'p + Ppk'^ 

+ {U - i) {p^k'p + p'pp\ 


/ ^ 

aJ 


(D.27) 


and Lap is again the leptonic tensor (2.17). The explicit form of (D.24) and (D.26) is 


1 

P 

= 4eV 

Q^su ’ 

(D.28) 

1 f 

= 4 6^6* . . . , 

Uu s 

(D.29) 


with 

F = + 2 {Q* + [s — (2 s' + U') xb] + xb {s' u + u U' 

-s s') + xl (2 s'2 + 2s' U' + U'^)}. (D.30) 

The interference term reads 

2jieAf,QEDCSAf,vcs. = 1 (D.Sl) 
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being the same as in (D.15) and 


^ Tr[/7“;i7^(/+ + i Tr[/7“;i7p(|i - ^')7^)] 

1.0 u 

Therefore we have, from (D.15), (D.31) and (D.32), 

G 


(D.32) 


2 Se = -4e»e’ 


Q‘^ s u t U S 


(D.33) 


with 


G = {—Q^ it {s + it) + Q‘^ {s — it) U' + (<§ + it) [s' ii 

+s (s' + U')] xb} {2 + s^ + -u^ — 2 s s' xb + 2 [s — (2 s' 

+U') xb] +2xb[u (s' + U') + (2 s'2 + 2 s' G' + U'^) xb]]- (D.34) 


The analytic form of the interference term agrees with [87] but differs from [88] in the 
massless case slightly, in particular in (15) of [88], 8 in the hrst line should be replaced 
by 4 and (—8) in the sixth line should be replaced by (—16). However we have checked 
that this does not affect our numerical results for HERA kinematics. 





Appendix E 


Matrix Elements for the Polarized 
QEDCS and VCS Processes 


In this Appendix we give the explicit expressions of the matrix elements relative to QED 
Compton scattering (QEDCS) and virtual Compton scattering (VCS) in ep e'jp and 
ep G'yX, with initial longitudinally polarized electron and proton. 


E.l Elastic Channel 

According to (7.23), 

I AMei 1^ = I 1^ + I AM,y^® f -2^e (E.l) 

is the matrix element squared of the process 

e(/) + p{P) ^ e(/') + 7(A:') + p{P'), (E.2) 

where the incoming electron and proton are longitudinally polarized. In terms of the 
Lorentz invariants (D.3) and (D.4), the spin dependent counterparts of (D.5) and (D.7) 
read 


I t = 1 £y(P; k'), (E.3) 

rW7 = I A('i P', k'), (E,4) 

and being given by (2.64) and (7.5) respectively, while 
CtfAP, S; n = - P’Y , 


(E.5) 
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where S is the spin four-vector of the proton, which satishes the conditions = —1 
and S' ■ P = 0. The tensor can be calculated from (D.8), inserting the polarization 
operator of the proton next to the term m) inside each trace. In (E.3) Pi is the 
Dirac form factor of the proton. Finally one obtains 

I 1^ = -^\-AA + ^ ApI p2(t), (E.6) 

tsu I ts' \ 

I AM^^® l' = [aA+ Apl p2(t), (E.7) 

I el I tU'S' y S's'U' J ^ ^ 

with 

AA = 2t^ + ("S — 2s' — 'u)(s -|- h) — 2t{s — 2s' — U') — 2uU' , (E.8) 

AB = —2t^ -t- s^ — siS -1- 2t^{2s + it) — t(3s^ -1- h^), (E.9) 

AP = (s -|- [—2s'^ -|- sh — 2s'u — -|- 2m^(s + u — t) — t{s' + h)] 

— (s -|- h) [2t^ — 3ts -f s^ -|- s(s' — 2u) + 3u{s' + h)] U' 

— [2t^ + — SU + 2u^ — t(3s + u)] U'‘^. (E.IO) 

We have used the notations s' = s — U' = U — mS, S' = S — mS for compactness. 
The interference term between QEDCS and VCS is given by 

2^e ^(P; P', A;'), (E.ll) 

with 

77/; l',k')=e^i^ + VhJ'Is 1 + 11''! 'r'M - Whe ]} (E.12) 

and 

P', k') = e3F,(/)F,(i)|i Tr|7y(f + m.)Y(r+ V+ m)7«(r+ m)7„l 

+ TTV[77(f + m)Y{r- ¥'+ m)7»(r+ m)7„]|. (E.13) 


4e® 

tsuiU' S' 


. ^ 2mSAE 

AD H- 

s' 


Flit) Frit), 


(E.14) 


Explicit calculation gives 

2 3fJe AMJ^^ 
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where AD and AE read: 

AD = [2t^ + (s — 2s' — u){s + u) — 2t{s — 2s' — U') — 2uU'] 

x{(s + u)\tu + s'{s + u)] + [t(s — h) + s{s + u)]U'}, (E.15) 

AE = [s(s — tY{s — 2t) + (2t^ — t^s — s^)u + (—2t^ — 3ts + + (t + 3s)-u^] U' 

— (s + u){—2Cu — (s + ■u)[s^(s' — 2u) + s'v^ + 2su{s' + u)] 

—t[—7ss'u + 'U^(—2s' + u) + s^(—3s' + Sh)] + t^[2-u(-u — s') + 5(5^ — 2s')]}. 

(E.16) 

E.2 Inelastic Channel 

Eor the corresponding inelastic channel the matrix elements relative to the electron-quark 
scattering process 

e{l) + g(p) ^ e(/') + -f{k') + g(p') (E.17) 

read: 

AM’Qedos i" = ± wig(p-,p')T‘‘l»'(l,l',k') (E,18) 

and 

|2 = (E.19) 

where and are given in (7.5) and in (2.64) respectively. Furthermore, is the 
antisymmetric part of the quark tensor, 

= -2ie^eleai3paP^p'', (E.20) 

Cq being the charge of the quark in units of the proton charge e, and 


An 

Tlp{j>]P,k)= 

{s - i)pp + (u - i)p'p\ {p - p'y. 

(E.21) 

The hnal results can be written as 



1 A aT^QEDCS 

1^ = 4eV 

(E.22) 

and 



1 f 

Usu 

(E.23) 
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where, in terms of the invariants defined in (D.3)-(D.23), 

AF = f -u^+ 2Q^+ 2Q‘^s-2xb[ss' + u{s +U') + Q‘^{2s +U')]. (E.24) 

The spin dependent counterpart of (D.32) reads 


2 VCS* _ _ 

2QH ’ 


(E.25) 


with being the same as in (E.12) and 


rf^ = e=>ediTl'|7''/7“/i7'’(/+ HYu] + i'IVl7"/7”;h7(/t - IfYY)]]- (E.26) 


The final result is given by 


2 


AG AH 
silt sty 


(E.27) 


AG = 2Q'^ + — 2ss'xb — u[u + 2{s' + U')xb\ + 2Q^[s — {2s' + U')xb], (E.28) 


AH = Q‘^[s{u — U)+ u{u + f/)] — xb{s + u) [su + s{s + U') . 


(E.29) 
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